Changeset 81 for trunk/mars/libf/phymars

Timestamp:

Feb 28, 2011, 6:53:54 PM (14 years ago)

Author:

aslmd

Message:

LMD_MM_MARS: changements mineurs cycle de l'eau [sorties plus physiques, du type vmr, mic_precip, ...]. quelques modifications cosmetiques des scripts d execution.

File:

• trunk/mars/libf/phymars/meso_physiq.F (modified) (1 diff)

trunk/mars/libf/phymars/meso_physiq.F

@@ -1613 +1613 @@
 ccc**************** WRF OUTPUT **************************
 ccc**************** WRF OUTPUT **************************
-      do ig=1,ngrid
-         wtsurf(ig) = tsurf(ig)    !! surface temperature
-         wco2ice(ig) = co2ice(ig)  !! co2 ice 
-
-         !!! specific to WRF WRF WRF
-         !!! just to output water ice on surface
-         !!! uncomment the Registry entry
-         IF (igcm_h2o_ice .ne. 0) qsurfice(ig) = qsurf(ig,igcm_h2o_ice)
-
-      enddo
+      !do ig=1,ngrid
+      !   wtsurf(ig) = tsurf(ig)    !! surface temperature
+      !   wco2ice(ig) = co2ice(ig)  !! co2 ice 
+      !
+      !   !!! specific to WRF WRF WRF
+      !   !!! just to output water ice on surface
+      !   !!! uncomment the Registry entry
+      !   IF (igcm_h2o_ice .ne. 0) qsurfice(ig) = qsurf(ig,igcm_h2o_ice)
+      !
+      !   !!! "VMR_ICE"   "VOL. MIXING RATIO ICE"           "ppm"
+      !   IF (igcm_h2o_ice .ne. 0) THEN
+      !     vmr=zq(1:ngridmx,1:nlayermx,igcm_h2o_ice)*mugaz/mmol(igcm_h2o_ice)
+      !   ENDIF
+      !
+      !enddo
+      wtsurf(1:ngrid) = tsurf(1:ngrid)    !! surface temperature
+      wco2ice(1:ngrid) = co2ice(1:ngrid)  !! co2 ice
+      mtot(1:ngrid) = mtot(1:ngrid) * 1.e6 / 917.
+      IF (igcm_h2o_ice .ne. 0) THEN      
+        qsurfice(1:ngrid) = qsurf(1:ngrid,igcm_h2o_ice)
+        vmr=1.e6 * zq(1:ngrid,1:nlayer,igcm_h2o_ice)
+     .           * mugaz / mmol(igcm_h2o_ice)
+      ENDIF
+
 c
 c THIS INCLUDE IS AUTOMATICALLY GENERATED FROM REGISTRY