Changeset 73 for trunk/mars

Timestamp:

Feb 20, 2011, 12:20:51 PM (14 years ago)

Author:

aslmd

Message:

LMD_MM_MARS: corrections cycle de l'eau propagees a la nouvelle physique. + corrections readmeteo.F90 [version synchronisee precedemment n etait pas la plus a jour] + corrections api.F90 pour avoir cp, R comme GCM

File:

• trunk/mars/libf/phymars/meso_physiq.F (modified) (3 diffs)

trunk/mars/libf/phymars/meso_physiq.F

@@ -192 +192 @@
       REAL ust(ngridmx)    !! pour LES avec isfflx!=0
       LOGICAL flag_LES     !! pour LES avec isfflx!=0
+      REAL qsurflast(ngridmx) !! pour diagnostics
       real alpha,lay1 ! coefficients for building layers
       integer iloop
@@ -1250 +1251 @@
 
       IF (tracer.AND.water.AND.(ngridmx.NE.1)) THEN
-         if (caps.and.(obliquit.lt.27.)) then
-           ! NB: Updated surface pressure, at grid point 'ngrid', is
-           !     ps(ngrid)=pplev(ngrid,1)+pdpsrf(ngrid)*ptimestep
-           tsurf(ngrid)=1./(1./136.27-r/5.9e+5*alog(0.0095*
-     &                     (pplev(ngrid,1)+pdpsrf(ngrid)*ptimestep)))
-         endif
+!         if (caps.and.(obliquit.lt.27.)) then
+!           ! NB: Updated surface pressure, at grid point 'ngrid', is
+!           !     ps(ngrid)=pplev(ngrid,1)+pdpsrf(ngrid)*ptimestep
+!           tsurf(ngrid)=1./(1./136.27-r/5.9e+5*alog(0.0095*
+!     &                     (pplev(ngrid,1)+pdpsrf(ngrid)*ptimestep)))
+!         endif
+!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
+!!!!! note WRF MESOSCALE AYMERIC -- mot cle "caps"
+!!!!! watercaptag n'est plus utilise que dans vdifc
+!!!!! ... pour que la sublimation ne soit pas stoppee 
+!!!!! ... dans la calotte permanente nord si qsurf=0
+!!!!! on desire garder cet effet regle par caps=T
+!!!!! on a donc commente "if (caps.and.(obliquit.lt.27.))" ci-dessus
+!!!!! --- remplacer ces lignes par qqch de plus approprie
+!!!!!      si on s attaque a la calotte polaire sud
+!!!!! pas d'autre occurrence majeure du mot-cle "caps"
+!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
+
 c       -------------------------------------------------------------
 c       Change of surface albedo (set to 0.4) in case of ground frost
@@ -1563 +1576 @@
          wtsurf(ig) = tsurf(ig)    !! surface temperature
          wco2ice(ig) = co2ice(ig)  !! co2 ice 
+
+         !!! TEMP TEMP TEMP TEMP TEMP TEMP TEMP
+         !!! specific to WRF WRF WRF
+         !!! just to output water ice on surface
+         !!! [it might not be water ice on surface but OK]
+         !!! uncomment the Registry entry
+         qsurflast(ig) = qsurf(ig,nqmx)
+
       enddo
 c