Changeset 73 for trunk/mesoscale/LMD_MM_MARS/SRC

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

Location:

trunk/mesoscale/LMD_MM_MARS/SRC

Files:

• POSTPROC/api.F90 (modified) (1 diff)
• PREP_MARS/readmeteo.F90 (modified) (8 diffs)
• WRFV2/mars_lmd/libf/phymars/meso_physiq.F (modified) (1 diff)
• WRFV2/phys/module_lmd_driver.F (modified) (4 diffs)

trunk/mesoscale/LMD_MM_MARS/SRC/POSTPROC/api.F90

@@ -48 +48 @@
       ! MARS CONSTANTS
       ! 
-      REAL, PARAMETER :: Rd  = 192.     ! gas constant m2 s-2 K-1
-      REAL, PARAMETER :: Cp  = 844.6    ! r= 8.314511E+0*1000.E+0/mugaz
-!REAL, PARAMETER :: Rd  = 191.0
-!REAL, PARAMETER :: Cp  = 744.5
+      !REAL, PARAMETER :: Rd  = 192.     ! gas constant m2 s-2 K-1
+      !REAL, PARAMETER :: Cp  = 844.6    ! r= 8.314511E+0*1000.E+0/mugaz
+      REAL, PARAMETER :: Rd  = 191.0
+      REAL, PARAMETER :: Cp  = 744.5
       REAL, PARAMETER :: RCP = Rd/Cp
       REAL, PARAMETER :: p0  = 610.

trunk/mesoscale/LMD_MM_MARS/SRC/PREP_MARS/readmeteo.F90

@@ -74 +74 @@
 !! Intermediate data arrays
 integer :: k,l,m,n,p
-real, dimension(:), allocatable :: lat,lon,time,alt,aps,bps,levels,dsoilvert
+real, dimension(:), allocatable :: lat,lon,time,alt,aps,bps,levels,vertdsoil
 real, dimension(:,:), allocatable :: vide,ones,ghtsfile
 real, dimension(:,:), allocatable :: interm
@@ -207 +207 @@
 allocate(dsoilfile(lonlen,latlen,altlen,timelen))
 allocate(isoilfile(lonlen,latlen,altlen,timelen))
-allocate(dsoilvert(altlen))
+allocate(vertdsoil(altlen))
 !allocate(tfileorig(lonlen,latlen,altlen,timelen))
 allocate(ufile(lonlen,latlen,altlen,timelen))
@@ -240 +240 @@
 isoilfile(:,:,:,:)=0 
 dsoilfile(:,:,:,:)=0 
-dsoilvert(:)=0.
+vertdsoil(:)=0.
 !tfileorig(:,:,:,:)=0
 !ufileorig(:,:,:,:)=0
@@ -591 +591 @@
 
     print *,'Surface Water ice'
-    ierr=NF_INQ_VARID(nid,"qsurf01",nvarid)
+!!!!!! ATTENTION ATTENTION 
+!!!!!! water ice a la surface est qsurf(ig,nqmx)
+    ierr=NF_INQ_VARID(nid,"qsurf02",nvarid)
     if (ierr.ne.NF_NOERR) then
-      write(*,*) "...No qsurf01 - surface Water ice set to 0"
+      write(*,*) "...No qsurf02 - surface Water ice set to 0"
       swatericefile(:,:,:)=0.
     else
@@ -648 +650 @@
         write(*,*) "...No soildepth - Set to -999"  !!! see soil_settings in LMD physics
         DO l=1,altlen
-                dsoilvert(l)=-999.
+                vertdsoil(l)=-999.
         ENDDO
     else
-        ierr=NF_GET_VAR_REAL(nid,nvarid,dsoilvert)
+        ierr=NF_GET_VAR_REAL(nid,nvarid,vertdsoil)
     endif
     print *, 'wait a minute' !! AS: I know this could be better
@@ -657 +659 @@
      DO n=1,latlen
       DO p=1,timelen
-       dsoilfile(m,n,:,p) = dsoilvert(:) 
+       dsoilfile(m,n,:,p) = vertdsoil(:) 
       ENDDO
      ENDDO
@@ -1050 +1052 @@
 !------------------------! 
 FIELD='SM100200'
-UNITS='kg/kg'
+UNITS='kg/m2'
 DESC='Surf water ice'
 XLVL=200100.
@@ -1470 +1472 @@
 deallocate(isoilfile)
 deallocate(dsoilfile)
-!deallocate(dsoilvert)
 !deallocate(tfileorig)
 deallocate(ufile)

trunk/mesoscale/LMD_MM_MARS/SRC/WRFV2/mars_lmd/libf/phymars/meso_physiq.F

@@ -1321 +1321 @@
        
       IF (tracer.AND.water.AND.ngridmx.NE.1) THEN
-         if (caps.and.(obliquit.lt.27.)) then
-           tsurf(ngrid)=1/(1/136.27-r/5.9e+5*alog(0.0095*ps(ngrid)))
-         endif
+         !if (caps.and.(obliquit.lt.27.)) then
+         !  tsurf(ngrid)=1/(1/136.27-r/5.9e+5*alog(0.0095*ps(ngrid)))
+         !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

trunk/mesoscale/LMD_MM_MARS/SRC/WRFV2/phys/module_lmd_driver.F

@@ -154 +154 @@
 !-------------------------------------------
    INTEGER ::    i,j,k,its,ite,jts,jte,ij
-   INTEGER ::    subs
+   INTEGER ::    subs,iii
 
    ! *** for LMD physics
@@ -623 +623 @@
     CASE(0)  !! NO TRACERS (mars=0)
       wtnom(nq) = 'co2'
+    CASE(1)
+      wtnom(1) = 'h2o_vap'  
+      wtnom(2) = 'h2o_ice'
+    CASE(2)
+      wtnom(1) = 'dust'      
     CASE(11) !! newwater mars==11 scalar:qh2o,qh2o_ice,qdust,qdust_number,qco2
       wtnom(nq) = 'co2'
@@ -665 +670 @@
 #ifdef NEWPHYS
 q_prof(:,1:nq) = SCALAR(i,kps:kpe,j,2:nq+1)  !! the names were set above !! one dummy tracer in WRF
-q_prof(:,nq) = 0.95 !! cas du CO2 !! TEMPORAIRE
+  !!! CAS DU CO2
+  DO iii=1,nq 
+   IF ( wtnom(iii) .eq. 'co2' ) q_prof(:,iii) = 0.95
+  ENDDO
 #else
 SELECT CASE (MARS_MODE)
@@ -790 +798 @@
 ! Tracer at surface !
 !-------------------!
-#ifdef NEWPHYS
-  !!! a faire !!!
-  PRINT *, 'WARNING WARNING no tracer at surface'
-  qsurf_val(:)=0.
-#else
 SELECT CASE (MARS_MODE)
     CASE(0)  !! NO TRACERS (mars=0)
     qsurf_val(:)=0.
     CASE(1)  !! WATER CYCLE (mars=1)
+#ifdef NEWPHYS    
+    qsurf_val(2)=MARS_WICE(i,j)    !! logique avec wtnom(2) = 'h2o_ice' defini ci-dessus
+    qsurf_val(1)=0.
+#else
     qsurf_val(2)=MARS_WICE(i,j)    !! attention... H2O ice is tracer nqmx in qsurf in LMD physics
-    qsurf_val(1)=0.                
+    qsurf_val(1)=0.
+#endif    
     CASE(2)  !! DUST CYCLE (mars=2)
     qsurf_val(:)=0.
+#ifdef NEWPHYS
+    CASE(11)
+    qsurf_val(:)=0. !! provisoire        
+#endif
 END SELECT
-#endif
 
 !-------------------!