Changeset 1529 for trunk/LMDZ.GENERIC/libf/phystd/radcommon_h.F90

Timestamp:

Apr 5, 2016, 10:51:51 AM (9 years ago)

Author:

emillour

Message:

Generic GCM: Towards a cleaner separation between physics and dynamics

• Got rid of references to "dimensions.h" from physics packages: use nbp_lon (=iim), nbp_lat (=jjp1) and nbp_lev from module mod_grid_phy_lmdz (in phy_common) instead.
• Removed module "comhdiff_mod.F90", as it is only used by module surf_heat_transp_mod.F90, moved module variables there.
• Addedin "surf_heat_transp_mod" local versions of some arrays and routines (from dyn3d) required to compute gradient, divergence, etc. on the global dynamics grid. As before, the slab ocean only works in serial.

EM

File:

• trunk/LMDZ.GENERIC/libf/phystd/radcommon_h.F90 (modified) (3 diffs)

trunk/LMDZ.GENERIC/libf/phystd/radcommon_h.F90

@@ -1 @@
-      module radcommon_h
-
-      use radinc_h
+module radcommon_h
+      use radinc_h, only: L_NSPECTI, L_NSPECTV, L_NGAUSS, NTstar, NTstop, &
+                          naerkind, nsizemax
       implicit none
 
@@ -121 +121 @@
       REAL :: omegaREFir(naerkind,nsizemax)
 
-      REAL tstellar ! Stellar brightness temperature (SW)
+      REAL,SAVE :: tstellar ! Stellar brightness temperature (SW)
 
-      real*8 planckir(L_NSPECTI,NTstop-NTstar+1)
+      real*8,save :: planckir(L_NSPECTI,NTstop-NTstar+1)
 
-      real*8 PTOP, TAUREF(L_LEVELS+1)
+      real*8,save :: PTOP
+      real*8,save,allocatable :: TAUREF(:)
 
-      real*8, parameter :: UBARI = 0.5D0
+      real*8,parameter :: UBARI = 0.5D0
 
-      real*8 gweight(L_NGAUSS)
+      real*8,save :: gweight(L_NGAUSS)
 !$OMP THREADPRIVATE(QREFvis,QREFir,omegaREFvis,omegaREFir,&     ! gweight read by master in sugas_corrk
                 !$OMP tstellar,planckir,PTOP,TAUREF)
@@ -144 +145 @@
       real*8, parameter :: grav = 6.672E-11
 
-      real*8 Cmk
-      save Cmk
-      real*8 glat_ig
-      save glat_ig
+      real*8,save :: Cmk
+      real*8,save :: glat_ig
 !$OMP THREADPRIVATE(Cmk,glat_ig)
 
       ! extinction of incoming sunlight (Saturn's rings, eclipses, etc...)
-      REAL, DIMENSION(:), ALLOCATABLE :: eclipse
+      REAL, DIMENSION(:), ALLOCATABLE ,SAVE :: eclipse
 
       !Latitude-dependent gravity
-      REAL, DIMENSION(:), ALLOCATABLE :: glat
+      REAL, DIMENSION(:), ALLOCATABLE , SAVE :: glat
 !$OMP THREADPRIVATE(glat,eclipse)
 
-      end module radcommon_h
+contains
+
+      subroutine ini_radcommon_h
+      use radinc_h, only: L_LEVELS
+      implicit none
+      
+        allocate(TAUREF(L_LEVELS+1))
+      
+      end subroutine ini_radcommon_h
+
+end module radcommon_h