Changeset 198 for trunk/MESOSCALE/LMD_MM_MARS

Timestamp:

Jul 5, 2011, 9:01:02 AM (13 years ago)

Author:

aslmd

Message:

MESOSCALE: Corrected the code so that only option monotonic is called when pd_scalar is true. in WRF v3.1, it is either PD or monotonic. have to put a more advanced flag perhaps, similar to what is done in recent version of WRF. other changes are purely cosmetics.

Location:

trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2

Files:

• dyn_em/module_em.F (modified) (1 diff)
• mars_lmd_new_storm/libf/phymars/aeropacity.F (modified) (1 diff)
• mars_lmd_new_storm/libf/phymars/meso_physiq.F (modified) (1 diff)
• mars_lmd_new_storm/libf/phymars/vdifc.F.modif (modified) (1 diff)

trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/dyn_em/module_em.F

@@ -835 +835 @@
       IF( (rk_step == 3) .and. pd_advection ) THEN
 
-        CALL advect_scalar_pd       ( scalar(ims,kms,jms,im),             &
-                                      scalar_old(ims,kms,jms,im),         &
-                                      advect_tend(ims,kms,jms),           &
-                                      ru, rv, ww, mut, mub, mu_old,       &
-                                      config_flags,                       &
-                                      msfu, msfv, msft, fnm, fnp,         &
-                                      rdx, rdy, rdnw,dt,                  &
-                                      ids, ide, jds, jde, kds, kde,       &
-                                      ims, ime, jms, jme, kms, kme,       &
-                                      its, ite, jts, jte, kts, kte     )
-
+!        CALL advect_scalar_pd       ( scalar(ims,kms,jms,im),             &
+!                                      scalar_old(ims,kms,jms,im),         &
+!                                      advect_tend(ims,kms,jms),           &
+!                                      ru, rv, ww, mut, mub, mu_old,       &
+!                                      config_flags,                       &
+!                                      msfu, msfv, msft, fnm, fnp,         &
+!                                      rdx, rdy, rdnw,dt,                  &
+!                                      ids, ide, jds, jde, kds, kde,       &
+!                                      ims, ime, jms, jme, kms, kme,       &
+!                                      its, ite, jts, jte, kts, kte     )
 
 

trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/mars_lmd_new_storm/libf/phymars/aeropacity.F

@@ -448 +448 @@
 
       iaer=1  !!!! PROVISOIRE !!!!
-      taulocref = 4.25 !10  ! ref optical depth of the local dust storm
+      taulocref = 10.   ! ref optical depth of the local dust storm
       ztoploc = 10      ! target pseudo-altitude of local storm (km)
-      radloc = 0.5       ! radius of dust storm (degree)
+      radloc = 4.0      ! radius of dust storm (degree)
       lonloc = 25       ! center longitude of storm (deg)
       latloc = -3.      ! center latitude of storm (deg)

trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/mars_lmd_new_storm/libf/phymars/meso_physiq.F

@@ -749 +749 @@
 !                write (*,*) "RDUST ", ig, l, rdust(ig-1,l), rdust(ig,l)
 !                rdust(ig,l) = rdust(ig,l-1)
+!               endif         
+!            ENDDO
+!         ENDDO
+!         DO l=2,nlayer
+!            DO ig=1,ngrid
+!               toto = abs(rdust(ig,l-1) - rdust(ig,l))/rdust(ig,l-1)
+!               if ( toto .gt. 0.45 ) then
+!                write (*,*) "RDUST ", ig, l, rdust(ig,l-1), rdust(ig,l)
+!                rdust(ig,l) = rdust(ig,l-1) 
+!     .                      + ( rdust(ig,l) - rdust(ig,l-1) ) / 2.
 !               endif         
 !            ENDDO

trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/mars_lmd_new_storm/libf/phymars/vdifc.F.modif

@@ -444 +444 @@
 
       if(tracer) then
-           PRINT*, 'alphavdifc', alpha_lift(igcm_dust_mass)
+           !PRINT*, 'alphavdifc', alpha_lift(igcm_dust_mass)
 c     Using the wind modified by friction for lifting and  sublimation
 c     ----------------------------------------------------------------