Changeset 967 for LMDZ4/trunk/libf/phylmd/phys_state_var_mod.F90

Timestamp:

Jun 11, 2008, 1:14:41 PM (16 years ago)

Author:

Laurent Fairhead

Message:

Pour des raisons pratiques (besoin de tester facilement les parametrisations
physiques en 1D et 3D), les arguments des routines de lecture et d'ecriture
de l'etat initial de la physique ont disparu des appels à ces routines et sont
maintenant passés par le module phys_state_var_mod
LF

File:

• LMDZ4/trunk/libf/phylmd/phys_state_var_mod.F90 (modified) (7 diffs)

LMDZ4/trunk/libf/phylmd/phys_state_var_mod.F90

@@ -10 +10 @@
       INTEGER, PARAMETER :: nout=3
       INTEGER, PARAMETER :: napisccp=1
+      INTEGER, SAVE :: radpas
       REAL, SAVE :: dtime, co2_ppm_etat0, solaire_etat0
+!$OMP THREADPRIVATE(radpas)
+!$OMP THREADPRIVATE(dtime, co2_ppm_etat0, solaire_etat0)
+
       REAL, ALLOCATABLE, SAVE :: rlat(:), rlon(:), pctsrf(:,:)
+!$OMP THREADPRIVATE(rlat, rlon, pctsrf)
       REAL, ALLOCATABLE, SAVE :: ftsol(:,:)
-      character(len=6), SAVE :: ocean
-      logical, SAVE :: ok_veget 
+!$OMP THREADPRIVATE(ftsol)
+!      character(len=6), SAVE :: ocean
+!!!!!!$OMP THREADPRIVATE(ocean)
+!      logical, SAVE :: ok_veget 
+!!!!!!$OMP THREADPRIVATE(ok_veget)
       REAL, ALLOCATABLE, SAVE :: falb1(:,:), falb2(:,:)
+!$OMP THREADPRIVATE(falb1, falb2)
       REAL, ALLOCATABLE, SAVE :: rain_fall(:), snow_fall(:)
+!$OMP THREADPRIVATE( rain_fall, snow_fall)
       REAL, ALLOCATABLE, SAVE :: solsw(:), sollw(:)
+!$OMP THREADPRIVATE(solsw, sollw)
       REAL, ALLOCATABLE, SAVE :: radsol(:)
+!$OMP THREADPRIVATE(radsol)
+
 !clesphy0 param physiq
 !
@@ -24 +37 @@
 !
       REAL, ALLOCATABLE, SAVE :: zmea(:), zstd(:), zsig(:), zgam(:)
+!$OMP THREADPRIVATE(zmea, zstd, zsig, zgam)
       REAL, ALLOCATABLE, SAVE :: zthe(:), zpic(:), zval(:)
+!$OMP THREADPRIVATE(zthe, zpic, zval)
 !     REAL tabcntr0(100)
       REAL, ALLOCATABLE, SAVE :: rugoro(:)
+!$OMP THREADPRIVATE(rugoro)
       REAL, ALLOCATABLE, SAVE :: t_ancien(:,:), q_ancien(:,:)
+!$OMP THREADPRIVATE(t_ancien, q_ancien)
       LOGICAL, SAVE :: ancien_ok
+!$OMP THREADPRIVATE(ancien_ok)
       REAL, ALLOCATABLE, SAVE :: clwcon(:,:),rnebcon(:,:)
+!$OMP THREADPRIVATE(clwcon,rnebcon)
       REAL, ALLOCATABLE, SAVE :: ratqs(:,:)
+!$OMP THREADPRIVATE(ratqs)
       REAL, ALLOCATABLE, SAVE :: pbl_tke(:,:,:) ! turb kinetic energy
+!$OMP THREADPRIVATE(pbl_tke)
       REAL, ALLOCATABLE, SAVE :: zmax0(:), f0(:) ! 
+!$OMP THREADPRIVATE(zmax0,f0)
       REAL, ALLOCATABLE, SAVE :: ema_work1(:,:), ema_work2(:,:)
+!$OMP THREADPRIVATE(ema_work1,ema_work2)
       REAL, ALLOCATABLE, SAVE :: entr_therm(:,:), fm_therm(:,:)
+!$OMP THREADPRIVATE(entr_therm,fm_therm)
       REAL, ALLOCATABLE, SAVE :: detr_therm(:,:)
+!$OMP THREADPRIVATE(detr_therm)
 !IM 150408
 !     pour phsystoke avec thermiques
       REAL,ALLOCATABLE,SAVE :: clwcon0th(:,:),rnebcon0th(:,:)
+!$OMP THREADPRIVATE(clwcon0th,rnebcon0th)
 ! radiation outputs
       REAL,ALLOCATABLE,SAVE :: swdn0(:,:), swdn(:,:)
+!$OMP THREADPRIVATE(swdn0,swdn)
       REAL,ALLOCATABLE,SAVE :: swup0(:,:), swup(:,:)
+!$OMP THREADPRIVATE(swup0,swup)
       REAL,ALLOCATABLE,SAVE :: SWdn200clr(:), SWdn200(:)
+!$OMP THREADPRIVATE(SWdn200clr,SWdn200)
       REAL,ALLOCATABLE,SAVE :: SWup200clr(:), SWup200(:)
+!$OMP THREADPRIVATE(SWup200clr,SWup200)
       REAL,ALLOCATABLE,SAVE :: lwdn0(:,:), lwdn(:,:)
+!$OMP THREADPRIVATE(lwdn0,lwdn)
       REAL,ALLOCATABLE,SAVE :: lwup0(:,:), lwup(:,:)
+!$OMP THREADPRIVATE(lwup0,lwup)
       REAL,ALLOCATABLE,SAVE :: LWdn200clr(:), LWdn200(:)
+!$OMP THREADPRIVATE(LWdn200clr,LWdn200)
       REAL,ALLOCATABLE,SAVE :: LWup200clr(:), LWup200(:)
+!$OMP THREADPRIVATE(LWup200clr,LWup200)
       REAL,ALLOCATABLE,SAVE :: LWdnTOA(:), LWdnTOAclr(:)
+!$OMP THREADPRIVATE(LWdnTOA,LWdnTOAclr)
 ! pressure level
       REAL,ALLOCATABLE,SAVE :: tsumSTD(:,:,:)
+!$OMP THREADPRIVATE(tsumSTD)
       REAL,ALLOCATABLE,SAVE :: usumSTD(:,:,:), vsumSTD(:,:,:)
+!$OMP THREADPRIVATE(usumSTD,vsumSTD)
       REAL,ALLOCATABLE,SAVE :: wsumSTD(:,:,:), phisumSTD(:,:,:)
+!$OMP THREADPRIVATE(wsumSTD,phisumSTD)
       REAL,ALLOCATABLE,SAVE :: qsumSTD(:,:,:), rhsumSTD(:,:,:)
+!$OMP THREADPRIVATE(qsumSTD,rhsumSTD)
       REAL,ALLOCATABLE,SAVE :: tnondef(:,:,:) 
+!$OMP THREADPRIVATE(tnondef)
       REAL,ALLOCATABLE,SAVE :: uvsumSTD(:,:,:)
+!$OMP THREADPRIVATE(uvsumSTD)
       REAL,ALLOCATABLE,SAVE :: vqsumSTD(:,:,:)
+!$OMP THREADPRIVATE(vqsumSTD)
       REAL,ALLOCATABLE,SAVE :: vTsumSTD(:,:,:)
+!$OMP THREADPRIVATE(vTsumSTD)
       REAL,ALLOCATABLE,SAVE :: wqsumSTD(:,:,:)
+!$OMP THREADPRIVATE(wqsumSTD)
       REAL,ALLOCATABLE,SAVE :: vphisumSTD(:,:,:)
+!$OMP THREADPRIVATE(vphisumSTD)
       REAL,ALLOCATABLE,SAVE :: wTsumSTD(:,:,:)
+!$OMP THREADPRIVATE(wTsumSTD)
       REAL,ALLOCATABLE,SAVE :: u2sumSTD(:,:,:)
+!$OMP THREADPRIVATE(u2sumSTD)
       REAL,ALLOCATABLE,SAVE :: v2sumSTD(:,:,:)
+!$OMP THREADPRIVATE(v2sumSTD)
       REAL,ALLOCATABLE,SAVE :: T2sumSTD(:,:,:)
+!$OMP THREADPRIVATE(T2sumSTD)
       INTEGER,ALLOCATABLE,SAVE :: seed_old(:,:)
+!$OMP THREADPRIVATE(seed_old)
       REAL,ALLOCATABLE,SAVE :: zuthe(:),zvthe(:)
+!$OMP THREADPRIVATE(zuthe,zvthe)
       REAL,ALLOCATABLE,SAVE :: alb_neig(:)
+!$OMP THREADPRIVATE(alb_neig)
 !cloud base mass flux
       REAL,ALLOCATABLE,SAVE :: ema_workcbmf(:), ema_cbmf(:)
+!$OMP THREADPRIVATE(ema_workcbmf,ema_cbmf)
 !cloud base pressure & cloud top pressure
       REAL,ALLOCATABLE,SAVE :: ema_pcb(:), ema_pct(:)
+!$OMP THREADPRIVATE(ema_pcb,ema_pct)
       REAL,ALLOCATABLE,SAVE :: Ma(:,:)        ! undilute upward mass flux
+!$OMP THREADPRIVATE(Ma)
       REAL,ALLOCATABLE,SAVE :: qcondc(:,:)    ! in-cld water content from convect
+!$OMP THREADPRIVATE(qcondc)
       REAL,ALLOCATABLE,SAVE :: wd(:) ! sb
+!$OMP THREADPRIVATE(wd)
       REAL,ALLOCATABLE,SAVE :: sigd(:)
+!$OMP THREADPRIVATE(sigd)
 !
       REAL,ALLOCATABLE,SAVE :: cin(:)
+!$OMP THREADPRIVATE(cin)
 ! ftd : differential heating between wake and environment
       REAL,ALLOCATABLE,SAVE :: ftd(:,:)
+!$OMP THREADPRIVATE(ftd)
 ! fqd : differential moistening between wake and environment
       REAL,ALLOCATABLE,SAVE :: fqd(:,:)     
+!$OMP THREADPRIVATE(fqd)
 !34EK
 ! -- Variables de controle de ALE et ALP
@@ -86 +147 @@
 !      convection d'Emanuel pour le declenchement et la regulation
       REAL,ALLOCATABLE,SAVE :: ALE(:)
+!$OMP THREADPRIVATE(ALE)
 !ALP : Puissance  disponible pour soulevement
       REAL,ALLOCATABLE,SAVE :: ALP(:)
+!$OMP THREADPRIVATE(ALP)
 !
 ! nouvelles variables pour le couplage convection-couche limite
       REAL,ALLOCATABLE,SAVE :: Ale_bl(:)
+!$OMP THREADPRIVATE(Ale_bl)
       REAL,ALLOCATABLE,SAVE :: Alp_bl(:)
+!$OMP THREADPRIVATE(Alp_bl)
       INTEGER,ALLOCATABLE,SAVE :: lalim_conv(:)
+!$OMP THREADPRIVATE(lalim_conv)
       REAL,ALLOCATABLE,SAVE :: wght_th(:,:)
+!$OMP THREADPRIVATE(wght_th)
 !
 ! variables de la wake
@@ -103 +170 @@
 ! dt_wake, dq_wake: LS tendencies due to wake
       REAL,ALLOCATABLE,SAVE :: wake_deltat(:,:)
+!$OMP THREADPRIVATE(wake_deltat)
       REAL,ALLOCATABLE,SAVE :: wake_deltaq(:,:)
+!$OMP THREADPRIVATE(wake_deltaq)
       REAL,ALLOCATABLE,SAVE :: wake_Cstar(:)
+!$OMP THREADPRIVATE(wake_Cstar)
       REAL,ALLOCATABLE,SAVE :: wake_s(:)
+!$OMP THREADPRIVATE(wake_s)
       REAL,ALLOCATABLE,SAVE :: wake_fip(:)
+!$OMP THREADPRIVATE(wake_fip)
       REAL,ALLOCATABLE,SAVE :: dt_wake(:,:)
+!$OMP THREADPRIVATE(dt_wake)
       REAL,ALLOCATABLE,SAVE :: dq_wake(:,:)
+!$OMP THREADPRIVATE(dq_wake)
 !
 ! pfrac_impa : Produits des coefs lessivage impaction
@@ -114 +188 @@
 ! pfrac_1nucl: Produits des coefs lessi nucl (alpha = 1) 
       REAL,ALLOCATABLE,SAVE :: pfrac_impa(:,:), pfrac_nucl(:,:)
+!$OMP THREADPRIVATE(pfrac_impa,pfrac_nucl)
       REAL,ALLOCATABLE,SAVE :: pfrac_1nucl(:,:)
+!$OMP THREADPRIVATE(pfrac_1nucl)
 !
       REAL,ALLOCATABLE,SAVE :: total_rain(:), nday_rain(:)  
+!$OMP THREADPRIVATE(total_rain,nday_rain)
       REAL,ALLOCATABLE,SAVE :: paire_ter(:)
+!$OMP THREADPRIVATE(paire_ter)
 ! albsol1: albedo du sol total pour SW visible
 ! albsol2: albedo du sol total pour SW proche IR
       REAL,ALLOCATABLE,SAVE :: albsol1(:), albsol2(:)
+!$OMP THREADPRIVATE(albsol1,albsol2)
       REAL,ALLOCATABLE,SAVE :: wo(:,:)
+!$OMP THREADPRIVATE(wo)
 !
 ! heat : chauffage solaire
@@ -132 +212 @@
 ! toplwdownclr : downward CS LW flux at TOA
       REAL,ALLOCATABLE,SAVE :: clwcon0(:,:),rnebcon0(:,:)
+!$OMP THREADPRIVATE(clwcon0,rnebcon0)
       REAL,ALLOCATABLE,SAVE :: heat(:,:)   
+!$OMP THREADPRIVATE(heat)
       REAL,ALLOCATABLE,SAVE :: heat0(:,:)
+!$OMP THREADPRIVATE(heat0)
       REAL,ALLOCATABLE,SAVE :: cool(:,:)
+!$OMP THREADPRIVATE(cool)
       REAL,ALLOCATABLE,SAVE :: cool0(:,:)
+!$OMP THREADPRIVATE(cool0)
       REAL,ALLOCATABLE,SAVE :: topsw(:), toplw(:)
+!$OMP THREADPRIVATE(topsw,toplw)
       REAL,ALLOCATABLE,SAVE :: sollwdown(:)
+!$OMP THREADPRIVATE(sollwdown)
       REAL,ALLOCATABLE,SAVE :: sollwdownclr(:)
+!$OMP THREADPRIVATE(sollwdownclr)
       REAL,ALLOCATABLE,SAVE :: toplwdown(:)
+!$OMP THREADPRIVATE(toplwdown)
       REAL,ALLOCATABLE,SAVE :: toplwdownclr(:)
+!$OMP THREADPRIVATE(toplwdownclr)
       REAL,ALLOCATABLE,SAVE :: topsw0(:),toplw0(:),solsw0(:),sollw0(:)
+!$OMP THREADPRIVATE(topsw0,toplw0,solsw0,sollw0)
       REAL,ALLOCATABLE,SAVE :: albpla(:)
+!$OMP THREADPRIVATE(albpla)
 ! pbase : cloud base pressure
 ! bbase : cloud base buoyancy
       REAL,ALLOCATABLE,SAVE :: cape(:)
+!$OMP THREADPRIVATE(cape)
       REAL,ALLOCATABLE,SAVE :: pbase(:)
+!$OMP THREADPRIVATE(pbase)
       REAL,ALLOCATABLE,SAVE :: bbase(:)
+!$OMP THREADPRIVATE(bbase)
 !
       REAL,SAVE,ALLOCATABLE :: zqasc(:,:)
+!$OMP THREADPRIVATE( zqasc)
       INTEGER,ALLOCATABLE,SAVE :: ibas_con(:), itop_con(:)
+!$OMP THREADPRIVATE(ibas_con,itop_con)
       REAL,SAVE,ALLOCATABLE :: rain_con(:)
+!$OMP THREADPRIVATE(rain_con)
       REAL,SAVE,ALLOCATABLE :: snow_con(:)
+!$OMP THREADPRIVATE(snow_con)
 !
 ! sulfate_pi : SO4 aerosol concentration [ug/m3] (pre-industrial value)
 
       REAL,SAVE,ALLOCATABLE :: sulfate_pi(:, :)
+!$OMP THREADPRIVATE(sulfate_pi)
       REAL,SAVE,ALLOCATABLE :: rlonPOS(:)
+!$OMP THREADPRIVATE(rlonPOS)
       REAL,SAVE,ALLOCATABLE :: newsst(:)
+!$OMP THREADPRIVATE(newsst)
       REAL,SAVE,ALLOCATABLE :: u10m(:,:), v10m(:,:)
+!$OMP THREADPRIVATE(u10m,v10m)
 !
 ! ok_ade=T -ADE=topswad-topsw
@@ -168 +271 @@
 !topswad, solswad : Aerosol direct effect
       REAL,SAVE,ALLOCATABLE :: topswad(:), solswad(:)
+!$OMP THREADPRIVATE(topswad,solswad)
 !topswai, solswai : Aerosol indirect effect
       REAL,SAVE,ALLOCATABLE :: topswai(:), solswai(:)
+!$OMP THREADPRIVATE(topswai,solswai)
       REAL,SAVE,ALLOCATABLE :: tau_ae(:,:,:), piz_ae(:,:,:)
+!$OMP THREADPRIVATE(tau_ae,piz_ae)
       REAL,SAVE,ALLOCATABLE :: cg_ae(:,:,:)
+!$OMP THREADPRIVATE(cg_ae)
 
 ! Les variables suivants uniquement pour un configuration avec INCA
 ! topswad_inca, solswad_inca : Aerosol direct effect
       REAL,SAVE,ALLOCATABLE :: topswad_inca(:), solswad_inca(:)
+!$OMP THREADPRIVATE(topswad_inca,solswad_inca)
 ! topswad0_inca, solswad0_inca : Aerosol direct effect
       REAL,SAVE,ALLOCATABLE :: topswad0_inca(:), solswad0_inca(:)
+!$OMP THREADPRIVATE(topswad0_inca,solswad0_inca)
 ! topswai_inca, solswai_inca : Aerosol indirect effect
       REAL,SAVE,ALLOCATABLE :: topswai_inca(:), solswai_inca(:)
+!$OMP THREADPRIVATE(topswai_inca,solswai_inca)
       REAL,SAVE,ALLOCATABLE :: topsw_inca(:,:), solsw_inca(:,:)
+!$OMP THREADPRIVATE(topsw_inca,solsw_inca)
       REAL,SAVE,ALLOCATABLE :: topsw0_inca(:,:), solsw0_inca(:,:)
+!$OMP THREADPRIVATE(topsw0_inca,solsw0_inca)
       REAL,SAVE,ALLOCATABLE :: tau_inca(:,:,:,:)
+!$OMP THREADPRIVATE(tau_inca)
       REAL,SAVE,ALLOCATABLE :: piz_inca(:,:,:,:)
+!$OMP THREADPRIVATE(piz_inca)
       REAL,SAVE,ALLOCATABLE :: cg_inca(:,:,:,:)
+!$OMP THREADPRIVATE(cg_inca)
       REAL,SAVE,ALLOCATABLE :: ccm(:,:,:)
+!$OMP THREADPRIVATE(ccm)
 
 CONTAINS