Changeset 159 for trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2

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)

File:

• trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/mars_lmd_new_storm/libf/phymars/meso_physiq.F (modified) (2 diffs)

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