Changeset 2146 for LMDZ5/trunk/libf/phylmd/phys_local_var_mod.F90

Timestamp:

Nov 14, 2014, 9:22:21 PM (10 years ago)

Author:

idelkadi

Message:

Les modifications introduites ont pour but :
1/ d'autoriser le couplage entre INCA-aerosol et les parametrisations de
la nouvelle physique (NP) de LMDZ, en particulier les thermiques et le
transport convectif,
2/ generaliser les routines de calcul de proprietes optiques des
aerosols pour RRTM au cas ou les aerosols sont interactifs
3/ d'inclure les effets LW des aerosols stratospheriques pour RRTM

File:

• LMDZ5/trunk/libf/phylmd/phys_local_var_mod.F90 (modified) (12 diffs)

LMDZ5/trunk/libf/phylmd/phys_local_var_mod.F90

@@ -96 +96 @@
       REAL, SAVE, ALLOCATABLE :: topswcf_aero(:,:),  solswcf_aero(:,:)  ! diag
       !$OMP THREADPRIVATE(topswcf_aero,solswcf_aero)
+! LW radiation diagnostics CK
+      REAL, SAVE, ALLOCATABLE :: toplwad_aero(:),  sollwad_aero(:)      ! diag
+      !$OMP THREADPRIVATE(toplwad_aero,sollwad_aero)
+      REAL, SAVE, ALLOCATABLE :: toplwai_aero(:),  sollwai_aero(:)      ! diag
+      !$OMP THREADPRIVATE(toplwai_aero,sollwai_aero)
+      REAL, SAVE, ALLOCATABLE :: toplwad0_aero(:), sollwad0_aero(:)     ! diag
+      !$OMP THREADPRIVATE(toplwad0_aero,sollwad0_aero)
 ! Special RRTM
       REAL, SAVE, ALLOCATABLE :: ZLWFT0_i(:,:),  ZSWFT0_i(:,:)      ! diag
@@ -140 +147 @@
       REAL, SAVE, ALLOCATABLE :: sconcso4(:) 
       !$OMP THREADPRIVATE(sconcso4) 
+      REAL, SAVE, ALLOCATABLE :: sconcno3(:) 
+      !$OMP THREADPRIVATE(sconcno3) 
       REAL, SAVE, ALLOCATABLE :: sconcoa(:) 
       !$OMP THREADPRIVATE(sconcoa) 
@@ -150 +159 @@
       REAL, SAVE, ALLOCATABLE :: concso4(:,:) 
       !$OMP THREADPRIVATE(concso4) 
+      REAL, SAVE, ALLOCATABLE :: concno3(:,:) 
+      !$OMP THREADPRIVATE(concno3) 
       REAL, SAVE, ALLOCATABLE :: concoa(:,:) 
       !$OMP THREADPRIVATE(concoa) 
@@ -197 +208 @@
 !$OMP THREADPRIVATE(topswcf_aerop, solswcf_aerop)
 
+! additional LW variables CK
+      REAL,ALLOCATABLE,SAVE :: toplwad_aerop(:), sollwad_aerop(:)
+!$OMP THREADPRIVATE(toplwad_aerop, sollwad_aerop)
+      REAL,ALLOCATABLE,SAVE :: toplwai_aerop(:), sollwai_aerop(:)
+!$OMP THREADPRIVATE(toplwai_aerop, sollwai_aerop)
+      REAL,ALLOCATABLE,SAVE :: toplwad0_aerop(:), sollwad0_aerop(:)
+!$OMP THREADPRIVATE(toplwad0_aerop, sollwad0_aerop)
 
 !Ajout de celles nécessaires au phys_output_write_mod
@@ -337 +355 @@
       allocate(topswai_aero(klon), solswai_aero(klon))
       allocate(topswad0_aero(klon), solswad0_aero(klon))
+     ! LW diagnostics CK
+      allocate(toplwad_aero(klon), sollwad_aero(klon))
+      allocate(toplwai_aero(klon), sollwai_aero(klon))
+      allocate(toplwad0_aero(klon), sollwad0_aero(klon))
+      ! end
       allocate(topsw_aero(klon,naero_grp), solsw_aero(klon,naero_grp))
       allocate(topsw0_aero(klon,naero_grp), solsw0_aero(klon,naero_grp))
@@ -363 +386 @@
       allocate(od550lt1aer(klon))                
       allocate(sconcso4(klon))
+      allocate(sconcno3(klon))
       allocate(sconcoa(klon))
       allocate(sconcbc(klon))
@@ -368 +392 @@
       allocate(sconcdust(klon))
       allocate(concso4(klon,klev))
+      allocate(concno3(klon,klev))
       allocate(concoa(klon,klev))
       allocate(concbc(klon,klev))
@@ -392 +417 @@
       ALLOCATE(solsw_aerop(klon,naero_grp), solsw0_aerop(klon,naero_grp))
       ALLOCATE(topswcf_aerop(klon,naero_grp), solswcf_aerop(klon,naero_grp))
+
+! additional LW variables CK
+      ALLOCATE(toplwad_aerop(klon), sollwad_aerop(klon))
+      ALLOCATE(toplwai_aerop(klon), sollwai_aerop(klon))
+      ALLOCATE(toplwad0_aerop(klon), sollwad0_aerop(klon))
 
 ! FH Ajout de celles nécessaires au phys_output_write_mod
@@ -496 +526 @@
       deallocate(topswai_aero,solswai_aero)
       deallocate(topswad0_aero,solswad0_aero)
+      ! LW additional CK
+      deallocate(toplwad_aero,sollwad_aero)
+      deallocate(toplwai_aero,sollwai_aero)
+      deallocate(toplwad0_aero,sollwad0_aero)
+      ! end
       deallocate(topsw_aero,solsw_aero)
       deallocate(topsw0_aero,solsw0_aero)
@@ -517 +552 @@
       deallocate(od550lt1aer)
       deallocate(sconcso4) 
+      deallocate(sconcno3) 
       deallocate(sconcoa) 
       deallocate(sconcbc) 
@@ -522 +558 @@
       deallocate(sconcdust) 
       deallocate(concso4) 
+      deallocate(concno3) 
       deallocate(concoa) 
       deallocate(concbc) 
@@ -549 +586 @@
       deallocate(topswcf_aerop, solswcf_aerop)
 
+!CK LW diagnostics
+      deallocate(toplwad_aerop, sollwad_aerop)
+      deallocate(toplwai_aerop, sollwai_aerop)
+      deallocate(toplwad0_aerop, sollwad0_aerop)
 
 ! FH Ajout de celles nécessaires au phys_output_write_mod