Changeset 1179 for LMDZ4/branches/LMDZ4-dev/libf/phylmd/phys_local_var_mod.F90

Timestamp:

Jun 11, 2009, 4:18:47 PM (15 years ago)

Author:

jghattas

Message:

• Ajout de l'interpolation vertical pour les nouveaux fichiers de forcage des aerosols. Utilisant les anciennes fichiers de SO4 pas d'interpolation possible. Convergence numerique avec la version precedente en utilisant les anciens fichiers des SO4. aerosol_optic.F90 change du nom pour readaerosol_optic.F90 (lecture d'aerosol + optic) Les fichiers de forcage aerosol doit maintenant avoir le suffix .nc.

• Correction des bugs pour inca et certain diagnostiques optionelles de radlwsw.
• Ajout de test pour le choix advection schema.

File:

• LMDZ4/branches/LMDZ4-dev/libf/phylmd/phys_local_var_mod.F90 (modified) (3 diffs)

LMDZ4/branches/LMDZ4-dev/libf/phylmd/phys_local_var_mod.F90

@@ -69 +69 @@
       REAL, SAVE, ALLOCATABLE :: d_ts(:,:), d_tr(:,:,:)
       !$OMP THREADPRIVATE(d_ts, d_tr)
+
+! diagnostique pour le rayonnement
+      REAL, SAVE, ALLOCATABLE :: topswad_aero(:),  solswad_aero(:)  ! diag
+      !$OMP THREADPRIVATE(topswad_aero,solswad_aero)
+      REAL, SAVE, ALLOCATABLE :: topswai_aero(:),  solswai_aero(:)  ! diag
+      !$OMP THREADPRIVATE(topswai_aero,solswai_aero)
+      REAL, SAVE, ALLOCATABLE :: topswad0_aero(:), solswad0_aero(:) ! pas utilise, eventuellment pour diag
+      !$OMP THREADPRIVATE(topswad0_aero,solswad0_aero)
+      REAL, SAVE, ALLOCATABLE :: topsw_aero(:,:),  solsw_aero(:,:)  ! pas utilise, eventuellment pour diag
+      !$OMP THREADPRIVATE(topsw_aero,solsw_aero)
+      REAL, SAVE, ALLOCATABLE :: topsw0_aero(:,:), solsw0_aero(:,:) ! pas utilise, eventuellment pour diag
+      !$OMP THREADPRIVATE(topsw0_aero,solsw0_aero)
 CONTAINS
 
@@ -100 +112 @@
       allocate(d_u_lif(klon,klev),d_v_lif(klon,klev))
       allocate(d_ts(klon,klev), d_tr(klon,klev,nbtr))
+      allocate(topswad_aero(klon), solswad_aero(klon))
+      allocate(topswai_aero(klon), solswai_aero(klon))
+      allocate(topswad0_aero(klon), solswad0_aero(klon))
+      allocate(topsw_aero(klon,9), solsw_aero(klon,9))
+      allocate(topsw0_aero(klon,9), solsw0_aero(klon,9))
       allocate(d_u_hin(klon,klev),d_v_hin(klon,klev),d_t_hin(klon,klev))
 
@@ -132 +149 @@
       deallocate(d_u_lif,d_v_lif)
       deallocate(d_ts, d_tr)
+      deallocate(topswad_aero,solswad_aero)
+      deallocate(topswai_aero,solswai_aero)
+      deallocate(topswad0_aero,solswad0_aero)
+      deallocate(topsw_aero,solsw_aero)
+      deallocate(topsw0_aero,solsw0_aero)
       deallocate(d_u_hin,d_v_hin,d_t_hin)