Changeset 2162 for trunk/LMDZ.MARS/libf/aeronomars

Timestamp:

Sep 24, 2019, 7:36:09 PM (5 years ago)

Author:

emillour

Message:

Mars GCM:
Some code cleanup (and preparing next steps):

• Turn calchim into a module and make tendencies module variables in calchim_mod and watercloud_mod
• Externalize in "physiq" the computation of solar zenithal angle (it should be computed at every physics timestep, regardless of iradia)

AB+EM

File:

• trunk/LMDZ.MARS/libf/aeronomars/calchim_mod.F90 (moved) (moved from trunk/LMDZ.MARS/libf/aeronomars/calchim.F90) (4 diffs)

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

@@ +1 @@
+MODULE calchim_mod
+
+  IMPLICIT NONE
+
+  REAL,SAVE,ALLOCATABLE :: zdqchim(:,:,:) ! Tendancy on pq due to photochemistry
+  REAL,SAVE,ALLOCATABLE :: zdqschim(:,:) ! Tendancy on qsurf due to photochemistry
+
+  CONTAINS
+
       subroutine calchim(ngrid,nlayer,nq,                           &
                          ptimestep,pplay,pplev,pt,pdt,dist_sol,mu0,         &
@@ -15 +24 @@
 
       use conc_mod, only: mmean ! mean molecular mass of the atmosphere
-      use comcstfi_h
-      use photolysis_mod
+      use comcstfi_h, only: pi
+      use photolysis_mod, only: jonline, init_photolysis
       use iono_h, only: temp_elect
 
@@ -69 +78 @@
 !=======================================================================
 
-#include "callkeys.h"
+include "callkeys.h"
 
 !     input:
@@ -637 +646 @@
       end if ! of if (output)
 
-      return
-      end
+      end subroutine calchim
+
+
+    subroutine ini_calchim_mod(ngrid,nlayer,nq)
+  
+      implicit none
+  
+      integer,intent(in) :: ngrid ! number of atmospheric columns
+      integer,intent(in) :: nlayer ! number of atmospheric layers
+      integer,intent(in) :: nq ! number of tracers
+
+      allocate(zdqchim(ngrid,nlayer,nq))
+      zdqchim(:,:,:)=0
+      allocate(zdqschim(ngrid,nq))
+      zdqschim(:,:)=0
+
+    end subroutine ini_calchim_mod
+
+
+    subroutine end_calchim_mod
+
+      implicit none
+
+      if (allocated(zdqchim))      deallocate(zdqchim)
+      if (allocated(zdqschim))      deallocate(zdqschim)
+
+    end subroutine end_calchim_mod
+
+END MODULE calchim_mod 
+