Changeset 648 for trunk/LMDZ.MARS/libf

Timestamp:

May 4, 2012, 10:24:20 AM (13 years ago)

Author:

emillour

Message:

Mars GCM:

• some syntax corrections in thermcall_main_mars, vdif_cd, pbl_parameters which cause problems when compiling with some strict compilers (g95, gfortran)
• added an initialisation of 'varian' in initracer for cases when using conrath dust; because that value can be is used elsewhere (e.g. surfacearea)

JYC+EM

Location:

trunk/LMDZ.MARS/libf/phymars

Files:

• initracer.F (modified) (1 diff)
• pbl_parameters.F (modified) (1 diff)
• thermcell_main_mars.F90 (modified) (1 diff)
• vdif_cd.F (modified) (1 diff)

trunk/LMDZ.MARS/libf/phymars/initracer.F

@@ -440 +440 @@
       else
 
+       ! initialize varian, which may be used (e.g. by surfacearea)
+       ! even with conrath dust
+       nueff_lift = 0.5
+       varian=sqrt(log(1.+nueff_lift))
+
        if (dustbin.gt.1) then
         print*,'initracer: STOP!',

trunk/LMDZ.MARS/libf/phymars/pbl_parameters.F

@@ -149 +149 @@
 !       zu2(ig)=MAX(pu(ig,1)*pu(ig,1)+pv(ig,1)*pv(ig,1)
 !     &                                ,(0.3*wstar_in(ig))**2)
-       zu2(ig)=pu(ig,1)*pu(ig,1)+pv(ig,1)*pv(ig,1)+
+       zu2(ig)=pu(ig,1)*pu(ig,1)+pv(ig,1)*pv(ig,1)
      &     + (log(1.+0.7*wstar_in(ig) + 2.3*wstar_in(ig)**2))**2
 

trunk/LMDZ.MARS/libf/phymars/thermcell_main_mars.F90

@@ -1241 +1241 @@
 
       if (1 .eq. 0) then
-      call thermcell_dqup(ngridmx,nlayermx,ptimestep     &
-     &     ,fm,entr,  &
-     &    masse,ztv,zdthladj)
+!      call thermcell_dqup(ngridmx,nlayermx,ptimestep     &
+!     &     ,fm,entr,  &
+!     &    masse,ztv,zdthladj)
       else
 

trunk/LMDZ.MARS/libf/phymars/vdif_cd.F

@@ -150 +150 @@
 !         zu2(ig)=MAX(pu(ig,1)*pu(ig,1) + pv(ig,1)*pv(ig,1),             &
 !     &      (0.3*wstar(ig))**2)
-          zu2(ig)=pu(ig,1)*pu(ig,1) + pv(ig,1)*pv(ig,1) +
+          zu2(ig)=pu(ig,1)*pu(ig,1) + pv(ig,1)*pv(ig,1) 
      &     + (log(1.+0.7*wstar(ig) + 2.3*wstar(ig)**2))**2