Changeset 802 for trunk/MESOSCALE/LMD_MM_MARS/SRC/PREP_MARS

Timestamp:

Oct 5, 2012, 2:52:43 AM (12 years ago)

Author:

aslmd

Message:

MESOSCALE : adaptation for polar runs with new physics. co2 mr is now taken from GCM and given as an input (previously was initialized at 0.95). added mars=10 which is a case with CO2 being the only tracer.

File:

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

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

@@ -88 +88 @@
 real, dimension(:,:,:), allocatable :: swatericefile!, swaterfile
 real, dimension(:,:,:,:), allocatable :: dustfile,dustnfile
+real, dimension(:,:,:,:), allocatable :: co2file
 real, dimension(:,:,:,:), allocatable :: ccnfile,ccnnfile
 
@@ -230 +231 @@
 allocate(dustfile(lonlen,latlen,altlen,timelen))
 allocate(dustnfile(lonlen,latlen,altlen,timelen))
+allocate(co2file(lonlen,latlen,altlen,timelen))
 allocate(ccnfile(lonlen,latlen,altlen,timelen))
 allocate(ccnnfile(lonlen,latlen,altlen,timelen))
@@ -285 +287 @@
 dustfile(:,:,:,:)=0
 dustnfile(:,:,:,:)=0
+co2file(:,:,:,:)=0
 ccnfile(:,:,:,:)=0
 ccnnfile(:,:,:,:)=0
@@ -636 +639 @@
 !! special water stuff
 !!------------------------
+
+    print *,'CO2 mass mixing ratio'
+    ierr=NF_INQ_VARID(nid,"co2",nvarid)
+    if (ierr.ne.NF_NOERR) then
+      write(*,*) "...No co2 - co2 mixing ratio set to 0.95"
+      co2file(:,:,:,:)=0.95
+    else
+      ierr=NF_GET_VAR_REAL(nid,nvarid,co2file)
+    endif
 
 !!------------------------
@@ -1230 +1242 @@
 !    ... Copy&Paste part !
 !------------------------!
+FIELD='CO2'
+UNITS='kg/kg'
+DESC='CO2 mixing ratio'
+XLVL=200100.
+SLAB=co2file(:,:,1,time_out(l))
+        ! And now put everything in the destination file
+        ! ... Header
+        write(1) IFV
+        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
+        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
+        ! ... Data
+        write(1) SLAB
+!print *,'The field '//DESC//' was written to '//output
+
+!------------------------!
+! >>> Write a variable   !
+!    ... Copy&Paste part !
+!------------------------!
 FIELD='DUSTQ'
 UNITS='kg/kg'
@@ -1533 +1563 @@
         XLVL=levels(k)
         SLAB=isoilfile(:,:,k,time_out(l))
+                ! And now put everything in the destination file
+                ! ... Header
+        write(1) IFV
+        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
+        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
+                ! ... Data
+                write(1) SLAB
+END DO
+!print *,'The field '//DESC//' was written to '//output
+
+!------------------------!
+! >>> Write a variable   !
+!    ... Copy&Paste part !
+!------------------------!
+FIELD='CO2'
+UNITS='kg/kg'
+DESC='CO2 mixing ratio'
+DO k = 1,altlen
+        XLVL=levels(k)
+        SLAB=co2file(:,:,k,time_out(l))
                 ! And now put everything in the destination file
                 ! ... Header
@@ -1676 +1726 @@
 deallocate(dustfile)
 deallocate(dustnfile)
+deallocate(co2file)
 deallocate(ccnfile)
 deallocate(ccnnfile)