Changeset 1194 for trunk/LMDZ.GENERIC/libf/phystd/callcorrk.F90

Timestamp:

Mar 4, 2014, 11:32:02 AM (11 years ago)

Author:

sglmd

Message:

Latitude-dependent gravity field added. Option oblate = .true. in callphys.def, and three additional variables needed: J2, Rmean and MassPlanet?.

File:

• trunk/LMDZ.GENERIC/libf/phystd/callcorrk.F90 (modified) (4 diffs)

trunk/LMDZ.GENERIC/libf/phystd/callcorrk.F90

@@ -180 +180 @@
       REAL*8 muvarrad(L_LEVELS)
 
-         
 !===============================================================
 !     Initialization on first call
@@ -227 +226 @@
          call suaer_corrk       ! set up aerosol optical properties
 
-         Cmk= 0.01 * 1.0 / (g * mugaz * 1.672621e-27) ! q_main=1.0 assumed
 
          if((igcm_h2o_vap.eq.0) .and. varactive)then
@@ -714 +712 @@
 
 
+         Cmk= 0.01 * 1.0 / (glat(ig) * mugaz * 1.672621e-27) ! q_main=1.0 assumed
+         glat_ig=glat(ig)
+
 !-----------------------------------------------------------------------
 !     Shortwave
 
          if(fract(ig) .ge. 1.0e-4) then ! only during daylight!
-!            if(sum(eclipse) .eq. 0.) then
-!                   write(*,*) 'but eclipse is still 0 in physiq !'
-!               endif
             if((ngrid.eq.1).and.(global1d))then
                do nw=1,L_NSPECTV
@@ -801 +799 @@
          DO l=2,L_NLAYRAD
             dtsw(ig,L_NLAYRAD+1-l)=(fmnetv(l)-fmnetv(l-1))  &
-                *g/(cpp*scalep*(plevrad(2*l+1)-plevrad(2*l-1)))
+                *glat(ig)/(cpp*scalep*(plevrad(2*l+1)-plevrad(2*l-1)))
             dtlw(ig,L_NLAYRAD+1-l)=(fmneti(l)-fmneti(l-1))  &
-                *g/(cpp*scalep*(plevrad(2*l+1)-plevrad(2*l-1)))
+                *glat(ig)/(cpp*scalep*(plevrad(2*l+1)-plevrad(2*l-1)))
          END DO      
 
 !     These are values at top of atmosphere
          dtsw(ig,L_NLAYRAD)=(fmnetv(1)-nfluxtopv)           &
-             *g/(cpp*scalep*(plevrad(3)-plevrad(1)))
+             *glat(ig)/(cpp*scalep*(plevrad(3)-plevrad(1)))
          dtlw(ig,L_NLAYRAD)=(fmneti(1)-nfluxtopi)           &
-             *g/(cpp*scalep*(plevrad(3)-plevrad(1)))
+             *glat(ig)/(cpp*scalep*(plevrad(3)-plevrad(1)))
 
 !     ---------------------------------------------------------------