Changeset 69 for trunk/mesoscale/LMD_MM_MARS/SRC/PREP_MARS

Timestamp:

Feb 18, 2011, 8:17:40 PM (14 years ago)

Author:

aslmd

Message:

README-000-svn: quelques ajouts

LMD_MM_MARS: serie de changements mineurs

A 0 mesoscale/TESTS/TESTCONVEC.tar.gz
ajout cas test convection bubble [tests nuages CO2]

A 0 mesoscale/LMD_MM_MARS/SIMU/deftank/POLAR_NEST_appere.tar.gz
M 68 mesoscale/LMD_MM_MARS/SRC/WPS/wps_mars/geogrid/GEOGRID.TBL.ARW_Mars
A 0 mesoscale/LMD_MM_MARS/SRC/WPS/wps_mars/geogrid/GEOGRID.TBL.ARW_Mars_POLARTESTS
ajout def avec 3 nests dans les regions polaires
changement de GEOGRID.TBL pour prendre en compte valeurs manquantes polaires pour ALB et TI

M 68 mesoscale/LMD_MM_MARS/SIMU/in_lmdz_mars/compile
M 68 mesoscale/LMD_MM_MARS/SIMU/in_lmdz_mars/myGCM/launch_gcm
M 68 mesoscale/LMDZ.MARS/myGCM/callphys.def
par defaut pointe desormais vers initialisation GCM avec cycle de l'eau

M 68 mesoscale/LMD_MM_MARS/SRC/PREP_MARS/readmeteo.F90
pour passer les champ de traceurs eau vapeur et glace a la surface

M 68 mesoscale/LMD_MM_MARS/SRC/WRFV2/mars_lmd/makegcm_*
M 68 mesoscale/LMD_MM_MARS/SRC/WRFV2/mars_lmd_new/makegcm_*
ajout de l'option csh -f pour importer les options de l'environnement courant

D 68 mesoscale/LMDZ.MARS/myGCM/startfi.nc
D 68 mesoscale/LMDZ.MARS/myGCM/start.nc
ne plus synchroniser (change a chaque run... et cree par launch_gcm de toute facon)

File:

• trunk/mesoscale/LMD_MM_MARS/SRC/PREP_MARS/readmeteo.F90 (modified) (7 diffs)

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

@@ -84 +84 @@
 real, dimension(:,:,:,:), allocatable :: tsoilfile
 real, dimension(:,:,:,:), allocatable :: waterfile, watericefile
+real, dimension(:,:,:), allocatable :: swaterfile, swatericefile
 
 !! Reading the parameter file 
@@ -210 +211 @@
 allocate(waterfile(lonlen,latlen,altlen,timelen))
 allocate(watericefile(lonlen,latlen,altlen,timelen))
+allocate(swaterfile(lonlen,latlen,timelen))
+allocate(swatericefile(lonlen,latlen,timelen))
 allocate(psfile(lonlen,latlen,timelen)) 
 allocate(tsfile(lonlen,latlen,timelen))
@@ -236 +239 @@
 waterfile(:,:,:,:)=0
 watericefile(:,:,:,:)=0
+swaterfile(:,:,:)=0
+swatericefile(:,:,:)=0
 psfile(:,:,:)=0 
 tsfile(:,:,:)=0
@@ -559 +564 @@
     endif
     ierr=NF_GET_VAR_REAL(nid,nvarid,watericefile)
+
+    print *,'Surface Water vapor'
+    ierr=NF_INQ_VARID(nid,"qsurf02",nvarid)
+    if (ierr.ne.NF_NOERR) then
+      write(*,*) "...No qsurf02 - surface Water vapor set to 0"
+      swaterfile(:,:,:)=0.
+    endif
+    ierr=NF_GET_VAR_REAL(nid,nvarid,swaterfile)
+
+    print *,'Surface Water ice'
+    ierr=NF_INQ_VARID(nid,"qsurf01",nvarid)
+    if (ierr.ne.NF_NOERR) then
+      write(*,*) "...No qsurf01 - surface Water ice set to 0"
+      swatericefile(:,:,:)=0.
+    endif
+    ierr=NF_GET_VAR_REAL(nid,nvarid,swatericefile)
 !!------------------------
 !! special water stuff
@@ -983 +1004 @@
 DESC='Water vapor'
 XLVL=200100.
-SLAB=waterfile(:,:,1,time_out(l))
+!SLAB=waterfile(:,:,1,time_out(l))
+SLAB=swaterfile(:,:,time_out(l))
         ! And now put everything in the destination file
         ! ... Header
@@ -1001 +1023 @@
 DESC='Water ice'
 XLVL=200100.
-SLAB=watericefile(:,:,1,time_out(l))
+!SLAB=watericefile(:,:,1,time_out(l))
+SLAB=swatericefile(:,:,time_out(l))
         ! And now put everything in the destination file
         ! ... Header
@@ -1238 +1261 @@
 deallocate(waterfile)
 deallocate(watericefile)
+deallocate(swaterfile)
+deallocate(swatericefile)
 deallocate(psfile) 
 deallocate(tsfile)