Changeset 1539 for LMDZ5/trunk/libf/phylmd/phys_state_var_mod.F90

Timestamp:

Jun 9, 2011, 12:13:33 AM (13 years ago)

Author:

musat

Message:

Ajouts CFMIP2/CMIP5

• 6eme fichier de sortie "stations" histstn.nc qui necessite 2 fichiers (voir DefLists?): npCFMIP_param.data(_*) contenant le nombre de points (120 pour simulations AMIP, 73 pour aqua) pointlocations.txt(_*) contenant le numero, les coordonnees (lon,lat) et le nom de chaque station
• flag LOGICAL dans tous les appels histwrite_phy pour pouvoir sortir le fichier histstn.nc

NB: 1) les flags de type phys_ que l'on met dans le physiq.def_L* pour ajouter plus de sorties

necessitent dorenavant 6 valeurs, la 6eme correspondant au fichier histstn.nc

2) par defaut le fichier histstn.nc ne sort pas; pour le sortir ajouter les lignes suivantes

dans physiq.def_L*

### Type de fichier : global (n) ou stations (y)
phys_out_filestations = n n n n n y

• introduction de 2 jeux de flags pour les taux des GES; taux actuels avec suffixes _act, taux futurs avec "_per" avec 2 appels au rayonnement si taux "_per" different des taux "_act" (utiles pour diags. CFMIP 4CO2)
• flags "betaCRF" pour calculs CRF pour experiences sensibilite proprietes optiques eau liquide nuageuse avec initialisations par defaut; sinon besoin de fichier beta_crf.data

Ajout flag LOGICAL lCOSP necessaire pour sortir un fichier stations
IM

File:

• LMDZ5/trunk/libf/phylmd/phys_state_var_mod.F90 (modified) (3 diffs)

LMDZ5/trunk/libf/phylmd/phys_state_var_mod.F90

@@ -280 +280 @@
       REAL,ALLOCATABLE,SAVE :: albpla(:)
 !$OMP THREADPRIVATE(albpla)
+
+!IM ajout variables CFMIP2/CMIP5
+      REAL,ALLOCATABLE,SAVE :: heatp(:,:), coolp(:,:)
+!$OMP THREADPRIVATE(heatp, coolp)
+      REAL,ALLOCATABLE,SAVE :: heat0p(:,:), cool0p(:,:)
+!$OMP THREADPRIVATE(heat0p, cool0p)
+      REAL,ALLOCATABLE,SAVE :: radsolp(:), topswp(:), toplwp(:)
+!$OMP THREADPRIVATE(radsolp, topswp, toplwp)
+      REAL,ALLOCATABLE,SAVE :: albplap(:)
+!$OMP THREADPRIVATE(albplap)
+      REAL,ALLOCATABLE,SAVE :: solswp(:), sollwp(:)
+!$OMP THREADPRIVATE(solswp, sollwp)
+      REAL,ALLOCATABLE,SAVE :: sollwdownp(:)
+!$OMP THREADPRIVATE(sollwdownp)
+      REAL,ALLOCATABLE,SAVE :: topsw0p(:),toplw0p(:)
+      REAL,ALLOCATABLE,SAVE :: solsw0p(:),sollw0p(:)
+!$OMP THREADPRIVATE(topsw0p,toplw0p,solsw0p,sollw0p)
+      REAL,ALLOCATABLE,SAVE :: lwdn0p(:,:), lwdnp(:,:)
+      REAL,ALLOCATABLE,SAVE :: lwup0p(:,:), lwupp(:,:)
+!$OMP THREADPRIVATE(lwdn0p, lwdnp, lwup0p, lwupp)
+      REAL,ALLOCATABLE,SAVE :: swdn0p(:,:), swdnp(:,:)
+      REAL,ALLOCATABLE,SAVE :: swup0p(:,:), swupp(:,:)
+!$OMP THREADPRIVATE(swdn0p, swdnp, swup0p, swupp)
+
 ! pbase : cloud base pressure
 ! bbase : cloud base buoyancy
@@ -446 +470 @@
       ALLOCATE(topsw0(klon),toplw0(klon),solsw0(klon),sollw0(klon))
       ALLOCATE(albpla(klon))
+!IM ajout variables CFMIP2/CMIP5
+      ALLOCATE(heatp(klon,klev), coolp(klon,klev))
+      ALLOCATE(heat0p(klon,klev), cool0p(klon,klev))
+      ALLOCATE(radsolp(klon), topswp(klon), toplwp(klon))
+      ALLOCATE(albplap(klon))
+      ALLOCATE(solswp(klon), sollwp(klon))
+      ALLOCATE(sollwdownp(klon))
+      ALLOCATE(topsw0p(klon),toplw0p(klon))
+      ALLOCATE(solsw0p(klon),sollw0p(klon))
+      ALLOCATE(lwdn0p(klon,klevp1), lwdnp(klon,klevp1))
+      ALLOCATE(lwup0p(klon,klevp1), lwupp(klon,klevp1))
+      ALLOCATE(swdn0p(klon,klevp1), swdnp(klon,klevp1))
+      ALLOCATE(swup0p(klon,klevp1), swupp(klon,klevp1))
+
       ALLOCATE(cape(klon))
       ALLOCATE(pbase(klon),bbase(klon))
@@ -539 +577 @@
       deallocate(topsw0,toplw0,solsw0,sollw0)
       deallocate(albpla)
+!IM ajout variables CFMIP2/CMIP5
+      deallocate(heatp, coolp)
+      deallocate(heat0p, cool0p)
+      deallocate(radsolp, topswp, toplwp)
+      deallocate(albplap)
+      deallocate(solswp, sollwp)
+      deallocate(sollwdownp)
+      deallocate(topsw0p,toplw0p)
+      deallocate(solsw0p,sollw0p)
+      deallocate(lwdn0p, lwdnp)
+      deallocate(lwup0p, lwupp)
+      deallocate(swdn0p, swdnp)
+      deallocate(swup0p, swupp)
       deallocate(cape)
       deallocate(pbase,bbase)