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

Timestamp:

Mar 1, 2012, 12:57:47 AM (13 years ago)

Author:

aslmd

Message:

MESOSCALE: for improved watercycle runs: ccn_mass and ccn_number are now passed through boundaries.

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 :: ccnfile,ccnnfile
 
 !! Reading the parameter file 
@@ -229 +230 @@
 allocate(dustfile(lonlen,latlen,altlen,timelen))
 allocate(dustnfile(lonlen,latlen,altlen,timelen))
+allocate(ccnfile(lonlen,latlen,altlen,timelen))
+allocate(ccnnfile(lonlen,latlen,altlen,timelen))
 allocate(psfile(lonlen,latlen,timelen)) 
 allocate(tsfile(lonlen,latlen,timelen))
@@ -282 +285 @@
 dustfile(:,:,:,:)=0
 dustnfile(:,:,:,:)=0
+ccnfile(:,:,:,:)=0
+ccnnfile(:,:,:,:)=0
 psfile(:,:,:)=0 
 tsfile(:,:,:)=0
@@ -651 +656 @@
     else
       ierr=NF_GET_VAR_REAL(nid,nvarid,dustnfile)
+    endif
+
+    print *,'CCN mass'
+    ierr=NF_INQ_VARID(nid,"ccn",nvarid)
+    if (ierr.ne.NF_NOERR) then
+      write(*,*) "...No ccn - CCN mass set to 0"
+      dustfile(:,:,:,:)=0.
+    else
+      ierr=NF_GET_VAR_REAL(nid,nvarid,ccnfile)
+    endif
+
+    print *,'CCN number'
+    ierr=NF_INQ_VARID(nid,"ccnN",nvarid)
+    if (ierr.ne.NF_NOERR) then
+      write(*,*) "...No ccnN - CCN number set to 0"
+      dustnfile(:,:,:,:)=0.
+    else
+      ierr=NF_GET_VAR_REAL(nid,nvarid,ccnnfile)
     endif
 !!------------------------
@@ -1241 +1264 @@
 !------------------------!
 ! >>> Write a variable   !
+!    ... Copy&Paste part !
+!------------------------!
+FIELD='CCNQ'
+UNITS='kg/kg'
+DESC='CCN mixing ratio'
+XLVL=200100.
+SLAB=ccnfile(:,:,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='CCNN'
+UNITS='part/kg'
+DESC='CCN number density'
+XLVL=200100.
+SLAB=ccnnfile(:,:,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   !
 !     PHOTOCHEMISTRY     !
 !------------------------!
@@ -1513 +1573 @@
         XLVL=levels(k)
         SLAB=dustnfile(:,:,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='CCNQ'
+UNITS='kg/kg'
+DESC='CCN mixing ratio'
+DO k = 1,altlen
+        XLVL=levels(k)
+        SLAB=ccnfile(:,:,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='CCNN'
+UNITS='part/kg'
+DESC='CCN number density'
+DO k = 1,altlen
+        XLVL=levels(k)
+        SLAB=ccnnfile(:,:,k,time_out(l))
                 ! And now put everything in the destination file
                 ! ... Header
@@ -1576 +1676 @@
 deallocate(dustfile)
 deallocate(dustnfile)
+deallocate(ccnfile)
+deallocate(ccnnfile)
 deallocate(psfile) 
 deallocate(tsfile)