Changeset 1888 for trunk/LMDZ.MARS/libf/aeronomars/calchim.F90

Timestamp:

Jan 9, 2018, 11:53:29 AM (7 years ago)

Author:

emillour

Message:

Mars GCM:
Updated the calculation of the dissociation and ionization
branching ratios, using the values in the Schunk and Nagy book.
New datafiles *branchingratio_schunkandnagy2000_param.dat must be loaded
for the "EUVDAT" subdirectory of the standard "datadir" directory.
Main changes are:

• param_v4_h.F90 -> New definition for the O2 ionization branching ratio
• param_read_e107.F -> Read the new files containing the S&N branching ratios
• paramfoto_compact.F -> Mainly cleaning of the code and the comments.

Also correction of a bug affecting the calculation of CO losses

• chemthermos.F90 -> Small modification to add the possibility of including

NO and O2 nightglow rates to the outputs

• calchim.F90 and calchim_asis.F90 -> account for change in arguments in

calls to chemthermos
FGG

File:

• trunk/LMDZ.MARS/libf/aeronomars/calchim.F90 (modified) (4 diffs)

trunk/LMDZ.MARS/libf/aeronomars/calchim.F90

@@ -166 +166 @@
       real :: surfdust1d(nlayer)   ! Dust surface area (cm2/cm3)
       real :: jo3(nlayer)          ! Photodissociation rate O3->O1D (s-1)
+      real :: em_no(nlayer)        !  NO nightglow emission rate
+      real :: em_o2(nlayer)        !  O2 nightglow emission rate      
 
 !     for output:
@@ -172 +174 @@
       parameter (output = .true.)
       real :: jo3_3d(ngrid,nlayer)  ! Photodissociation rate O3->O1D (s-1)
+      real :: emission_no(ngrid,nlayer) !NO emission rate
+      real :: emission_o2(ngrid,nlayer) !O2 emission rate
 
 !=======================================================================
@@ -665 +669 @@
          if (thermochem) then
             call chemthermos(ig,nlayer,lswitch,chemthermod,zycol,ztemp,&
-                             zdens,zpress,zlocal,szacol,ptimestep,zday)
+                             zdens,zpress,zlocal,szacol,ptimestep,zday,&
+                             em_no,em_o2)
+            do l=1,nlayer
+               emission_no(ig,l)=em_no(l)
+               emission_o2(ig,l)=em_o2(l)
+            enddo
          end if
 
@@ -715 +724 @@
             call writediagfi(ngrid,'jo3','j o3->o1d',    &
                              's-1',3,jo3_3d(1,1))
+            call writediagfi(ngrid,'emission_no',        &
+                 'NO nightglow emission rate','cm-3 s-1',3,emission_no)
+            call writediagfi(ngrid,'emission_o2',        &
+                 'O2 nightglow emission rate','cm-3 s-1',3,emission_o2)
+
            if (callstats) then
               call wstats(ngrid,'jo3','j o3->o1d',       &
                           's-1',3,jo3_3d(1,1))
+              call wstats(ngrid,'emission_no',           &
+                   'NO nightglow emission rate','cm-3 s-1',3,emission_no)
+              call wstats(ngrid,'emission_o2',           &
+                   'O2 nightglow emission rate','cm-3 s-1',3,emission_o2)
            endif
          end if ! of if (ngrid.gt.1)