/[MITgcm]/MITgcm/pkg/ebm/ebm_atmosphere.F
ViewVC logotype

Diff of /MITgcm/pkg/ebm/ebm_atmosphere.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.3 by heimbach, Fri May 21 21:45:35 2004 UTC revision 1.4 by heimbach, Wed Jul 28 19:54:37 2004 UTC
# Line 165  CADJ STORE DTDy35(:) = comlev1_bibj, key Line 165  CADJ STORE DTDy35(:) = comlev1_bibj, key
165  #endif  #endif
166    
167          DO no_so = 1, 2          DO no_so = 1, 2
168           gamma = -T35(no_so)*beta*Hw*Nw*Nw/           IF ( DTDy35(no_so).NE.0. .AND. T35(no_so).NE.0. ) THEN
169              gamma = -T35(no_so)*beta*Hw*Nw*Nw/
170       <        (gravity*f0*DTDy35(no_so))       <        (gravity*f0*DTDy35(no_so))
171           kappa = Hw/(1 + gamma)            kappa = Hw/(1 + gamma)
172           De = Hw/(0.48 + 1.48*gamma)            De = Hw/(0.48 + 1.48*gamma)
173           C = 0.6*gravity*kappa*kappa*Nw/            C = 0.6*gravity*kappa*kappa*Nw/
174       <        (Tw*f0*f0)       <        (Tw*f0*f0)
175           Cs = rho_air*cp*C*            Cs = rho_air*cp*C*
176       <        (1/(1/Hw+1/De) - 1/(1/Hw+1/De+1/dz))       <        (1/(1/Hw+1/De) - 1/(1/Hw+1/De+1/dz))
177           Cf = htil*2.97e12*C/(T35(no_so)**3)*(            Cf = htil*2.97e12*C/(T35(no_so)**3)*(
178       <        1/(1/De + (5420*tau /(T35(no_so)**2)))       <        1/(1/De + (5420*tau /(T35(no_so)**2)))
179       <     - 1/(1/De+5420*tau/(T35(no_so)**2)+1/dz))       <        - 1/(1/De+5420*tau/(T35(no_so)**2)+1/dz))
180           Cl = Cf*lv            Cl = Cf*lv
181           Hd35(no_so) = 2.*PI*rSphere*COS(lat(2)*deg2rad)            Hd35(no_so) = 2.*PI*rSphere*COS(lat(2)*deg2rad)
182       <        *(Cs + Cl*exp(-5420./T35(no_so)))       <        *(Cs + Cl*exp(-5420./T35(no_so)))
183       <        *(abs(DTDy35(no_so))**trans_eff)       <        *(abs(DTDy35(no_so))**trans_eff)
184           Fw35(no_so) = 2.*PI*rSphere*COS(lat(2)*deg2rad)            Fw35(no_so) = 2.*PI*rSphere*COS(lat(2)*deg2rad)
185       <        *(abs(DTDy35(no_so))**trans_eff)       <        *(abs(DTDy35(no_so))**trans_eff)
186       <        *Cf*exp(-5420./T35(no_so))       <        *Cf*exp(-5420./T35(no_so))
187  c        write(0,*) no_so, Hd35(no_so), Fw35(no_so)           ELSE
188              Hd35(no_so) = 0.
189              Fw35(no_so) = 0.
190             ENDIF
191          ENDDO          ENDDO
192    c
193          Fw35(1) = 929944128.          Fw35(1) = 929944128.
194          Fw35(2) = 678148032.          Fw35(2) = 678148032.
195    c
196  #ifdef EBM_VERSION_1BASIN  #ifdef EBM_VERSION_1BASIN
197  c      Fw35(2) = 0.7*Fw35(2)  c      Fw35(2) = 0.7*Fw35(2)
198  #else  #else
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22