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

Timestamp:

Oct 21, 2011, 11:31:28 AM (14 years ago)

Author:

aslmd

Message:

MESOSCALE : Preparatory commit for the ultimate option mars=42 which

would allow mesoscale modeling with photochemistry.

[see r76 method to add 'mars' options]
Modified module_lmd_driver.F and Registry.EM and runmeso
Modified readmeteo.F90 and introduced an option -DPHOTOCHEM
...
Transparent to the casual user
Option mars=42 not yet finished though -- so do not use!

Location:

trunk/MESOSCALE/LMD_MM_MARS/SRC/PREP_MARS

Files:

• compile_pgf (modified) (1 diff)
• readmeteo.F90 (modified) (8 diffs)

trunk/MESOSCALE/LMD_MM_MARS/SRC/PREP_MARS/compile_pgf

@@ -8 +8 @@
 -L$NETCDF/lib -lnetcdf \
 -I$NETCDF/include \
--o readmeteo.exe
+-o readmeteo.exe #-DPHOTOCHEM
 
 pgf90 create_readmeteo.F90 \
 -L$NETCDF/lib -lnetcdf \
 -I$NETCDF/include \
--o create_readmeteo.exe
+-o create_readmeteo.exe #-DPHOTOCHEM \
 
 \rm fix_no_info.inc 2> /dev/null

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

@@ -98 +98 @@
 character*13, dimension(:), allocatable :: date_out
 character*19, dimension(:), allocatable :: date_out2
+
+#ifdef PHOTOCHEM
+real, dimension(:,:,:,:,:), allocatable :: chemtrac
+integer :: nchemtrac,i
+CHARACTER*20,DIMENSION(:),ALLOCATABLE :: wtnom
+#endif
+
+
 !***************************************************************************
 !***************************************************************************
@@ -230 +238 @@
 allocate(interm(lonlen,latlen))
 allocate(gwparam(lonlen,latlen,5))
-allocate(ghtsfile(lonlen,latlen))       !! no scan axis
+allocate(ghtsfile(lonlen,latlen))    !! no scan axis
 allocate(vide(lonlen,latlen))
 allocate(ones(lonlen,latlen))
 allocate(lat(latlen), lon(lonlen), alt(altlen), time(timelen))
 allocate(aps(altlen),bps(altlen),levels(altlen))
+#ifdef PHOTOCHEM
+nchemtrac = 14
+allocate(wtnom(nchemtrac))
+wtnom(1)  = "c_co2"
+wtnom(2)  = "c_co"
+wtnom(3)  = "c_o"
+wtnom(4)  = "c_o1d"
+wtnom(5)  = "c_o2"
+wtnom(6)  = "c_o3"
+wtnom(7)  = "c_h"
+wtnom(8)  = "c_h2"
+wtnom(9)  = "c_oh"
+wtnom(10) = "c_ho2"
+wtnom(11) = "c_h2o2"
+wtnom(12) = "c_ch4"
+wtnom(13) = "c_n2"
+wtnom(14) = "c_ar"
+allocate(chemtrac(lonlen,latlen,altlen,timelen,nchemtrac))
+chemtrac(:,:,:,:,:)=0
+#endif
 
 tfile(:,:,:,:)=0
@@ -523 +551 @@
    IF (ierr .NE. NF_NOERR) THEN
       ierr = NF_INQ_VARID (nid,"t",nvarid)
-        IF (ierr .NE. NF_NOERR) THEN
+        IF (ierr .NE. NF_NOERR) THEN
           PRINT *, "Error: Readmeteo <t> not found"
           stop
-        ENDIF
+        ENDIF
    ENDIF
 #ifdef NC_DOUBLE
@@ -627 +655 @@
 !! special dust stuff
 !!------------------------
-
 
 !SELECT CASE(ident)
@@ -677 +704 @@
 !!!!!!!!
 !!!!!!!! new physics
+
+
+!!!!!!!!!!!!!!!!!!!!!!!!NEW PHYSICS + PHOTOCHEM
+!!!!!!!!!!!!!!!!!!!!!!!!NEW PHYSICS + PHOTOCHEM
+#ifdef PHOTOCHEM
+    print *,'photochem'
+    DO i=1,nchemtrac
+     print *,wtnom(i)
+     ierr=NF_INQ_VARID(nid,wtnom(i),nvarid)
+     if (ierr.ne.NF_NOERR) then
+       write(*,*) "...No ",wtnom(i), " - set to 0"
+       chemtrac(:,:,:,:,i)=0.
+     else
+       ierr=NF_GET_VAR_REAL(nid,nvarid,chemtrac(:,:,:,:,i))
+     endif
+    ENDDO
+#endif
+!!!!!!!!!!!!!!!!!!!!!!!!NEW PHYSICS + PHOTOCHEM
+!!!!!!!!!!!!!!!!!!!!!!!!NEW PHYSICS + PHOTOCHEM
+
 
 
@@ -1192 +1239 @@
 !print *,'The field '//DESC//' was written to '//output
 
+!------------------------!
+! >>> Write a variable   !
+!     PHOTOCHEMISTRY     !
+!------------------------!
+#ifdef PHOTOCHEM
+    DO i=1,nchemtrac
+       FIELD=wtnom(i)
+       UNITS='units'
+       DESC='desc'
+       XLVL=200100.
+       SLAB=chemtrac(:,:,1,time_out(l),i)
+       ! 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
+    ENDDO
+#endif
 
 !!----------------------------------------------------
@@ -1455 +1522 @@
 END DO
 !print *,'The field '//DESC//' was written to '//output
+
+
+!------------------------!
+! >>> Write a variable   !
+!     PHOTOCHEMISTRY     !
+!------------------------!
+#ifdef PHOTOCHEM
+    DO i=1,nchemtrac
+       FIELD=wtnom(i)
+       UNITS='units'
+       DESC='desc'
+       DO k = 1,altlen
+         XLVL=levels(k)
+         SLAB=chemtrac(:,:,k,time_out(l),i)
+         ! 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
+    ENDDO
+#endif
 
 print *,'****done file '//output, int(100.*float(l)/float(FILES)), ' % '
@@ -1498 +1589 @@
 deallocate(aps,bps,levels)
 
+#ifdef PHOTOCHEM
+deallocate(chemtrac)
+deallocate(wtnom)
+#endif
 
 print *, '------------------------'