Changeset 159 for trunk/MESOSCALE/LMD_MM_MARS

Timestamp:

Jun 16, 2011, 11:43:20 AM (13 years ago)

Author:

aslmd

Message:

MESOSCALE: modifications mineures LESnophys (stretch_grid) et mars_lmd_new_storm (essais pour stabiliser RDUST)

Location:

trunk/MESOSCALE/LMD_MM_MARS/SRC

Files:

• LESnophys_/modif/module_initialize_les.F (modified) (2 diffs)
• WRFV2/mars_lmd_new_storm/libf/phymars/meso_physiq.F (modified) (2 diffs)

trunk/MESOSCALE/LMD_MM_MARS/SRC/LESnophys_/modif/module_initialize_les.F

@@ -164 +164 @@
 !  stretch_grid = .true.
 !  FOR LES, set stretch to false
+#ifdef NOPHYS
+   stretch_grid = .true.
+#else
    stretch_grid = .false.
+#endif
    delt = 3.
 !   z_scale = .50
@@ -233 +237 @@
 
    IF (stretch_grid) THEN ! exponential stretch for eta (nearly constant dz)
+     print *, '*** STRETCH GRID ***' 
      DO k=1, kde
       grid%znw(k) = (exp(-(k-1)/float(kde-1)/z_scale) - exp(-1./z_scale))/ &

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

@@ -211 +211 @@
       logical tracerdyn
       CHARACTER (len=20) :: wtnom(nqmx) ! tracer name
+      REAL toto
 
 c   outputs:
@@ -725 +726 @@
      $     zdtlw,zdtsw,fluxsurf_lw,fluxsurf_sw,fluxtop_lw,fluxtop_sw,
      &     tauref,tau,aerosol,ccn,rdust,rice,nuice,zdqnorm,dsodust)
+
+!!!!!!!!!!!!!!!!!!!!!! TEST TEST
+!         DO l=2,nlayer
+!            DO ig=1,ngrid
+!               toto = abs(rdust(ig,l) - rdust(ig,l-1))
+!               if ( toto .gt. 0.5*rdust(ig,l) ) then
+!                    write (*,*) "LIMITING RADIUS", ig, l
+!                    write (*,*) rdust(ig,l-1), rdust(ig,l)
+!                   rdust(ig,l) = rdust(ig,l-1)
+!               endif         
+!            ENDDO
+!         ENDDO
+!
+!         DO l=2,nlayer
+!            DO ig=2,ngrid
+!               toto = abs(rdust(ig-1,l) - rdust(ig,l))/rdust(ig,l)
+!               if ( toto .gt. 0.45 ) then
+!                write (*,*) "RDUST ", ig, l, rdust(ig-1,l), rdust(ig,l)
+!                rdust(ig,l) = rdust(ig,l-1)
+!               endif         
+!            ENDDO
+!         ENDDO