Changeset 4088 for LMDZ6/trunk


Ignore:
Timestamp:
Mar 10, 2022, 8:03:20 AM (2 years ago)
Author:
fhourdin
Message:

Correction isotopes + version identique sur phylmd

Location:
LMDZ6/trunk/libf
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • LMDZ6/trunk/libf/phylmd/phys_state_var_mod.F90

    r4059 r4088  
    6262      REAL, ALLOCATABLE, SAVE :: swradcorr(:)
    6363!$OMP THREADPRIVATE(swradcorr)
     64#ifdef ISO
     65      REAL,ALLOCATABLE,SAVE :: xtsol(:,:),fxtevap(:,:,:)
     66!$OMP THREADPRIVATE(xtsol,fxtevap)
     67      REAL, ALLOCATABLE, SAVE :: xtrain_fall(:,:), xtsnow_fall(:,:)
     68!$OMP THREADPRIVATE(xtrain_fall,xtsnow_fall)
     69#endif
    6470
    6571!clesphy0 param physiq
     
    8086      REAL, ALLOCATABLE, SAVE :: prw_ancien(:), prlw_ancien(:), prsw_ancien(:)
    8187!$OMP THREADPRIVATE(prw_ancien, prlw_ancien, prsw_ancien)
     88#ifdef ISO
     89      REAL, ALLOCATABLE, SAVE :: xt_ancien(:,:,:),xtl_ancien(:,:,:),xts_ancien(:,:,:)
     90!$OMP THREADPRIVATE(xt_ancien,xtl_ancien,xts_ancien)
     91#endif
    8292      REAL, ALLOCATABLE, SAVE :: u_ancien(:,:), v_ancien(:,:)
    8393!$OMP THREADPRIVATE(u_ancien, v_ancien)
     
    221231      REAL,ALLOCATABLE,SAVE :: fqd(:,:)     
    222232!$OMP THREADPRIVATE(fqd)
     233#ifdef ISO
     234      REAL, ALLOCATABLE, SAVE :: fxtd(:,:,:)
     235!$OMP THREADPRIVATE(fxtd)
     236#endif
    223237!34EK
    224238! -- Variables de controle de ALE et ALP
     
    260274      REAL,ALLOCATABLE,SAVE :: wake_deltaq(:,:)
    261275!$OMP THREADPRIVATE(wake_deltaq)
     276#ifdef ISO
     277      REAL, ALLOCATABLE, SAVE :: wake_deltaxt(:,:,:)
     278!$OMP THREADPRIVATE(wake_deltaxt)
     279#endif
    262280      REAL,ALLOCATABLE,SAVE :: wake_s(:)
    263281!$OMP THREADPRIVATE(wake_s)
     
    397415!$OMP THREADPRIVATE(snow_con)
    398416!
     417#ifdef ISO
     418      REAL,SAVE,ALLOCATABLE :: xtrain_con(:,:)
     419!$OMP THREADPRIVATE(xtrain_con)
     420      REAL,SAVE,ALLOCATABLE :: xtsnow_con(:,:)
     421!$OMP THREADPRIVATE(xtsnow_con)
     422#endif
    399423      REAL,SAVE,ALLOCATABLE :: rlonPOS(:)
    400424!$OMP THREADPRIVATE(rlonPOS)
     
    431455!$OMP THREADPRIVATE(ratqs_inter)
    432456
     457#ifdef ISO
     458#ifdef ISOTRAC
     459      INTEGER,SAVE,ALLOCATABLE :: bassin_map(:)
     460!$OMP THREADPRIVATE(bassin_map)
     461      INTEGER,SAVE,ALLOCATABLE :: boite_map(:,:)
     462!$OMP THREADPRIVATE(boite_map)
     463#endif   
     464#endif
    433465      REAL, ALLOCATABLE, SAVE:: du_gwd_rando(:, :), du_gwd_front(:, :)
    434466      !$OMP THREADPRIVATE(du_gwd_rando, du_gwd_front)
     
    466498USE aero_mod
    467499USE infotrac_phy, ONLY : nbtr
     500#ifdef ISO
     501USE infotrac_phy, ONLY : ntraciso,niso
     502#endif
    468503USE indice_sol_mod
    469504use config_ocean_skin_m, only: activate_ocean_skin
     
    659694      ALLOCATE(cg_aero_lw_rrtm(klon,klev,2,nbands_lw_rrtm))
    660695      ALLOCATE(ccm(klon,klev,nbands))
     696
     697#ifdef ISO
     698      ALLOCATE(xtsol(niso,klon),fxtevap(ntraciso,klon,nbsrf))
     699      ALLOCATE(fxtd(ntraciso,klon,klev))
     700      ALLOCATE(wake_deltaxt(ntraciso,klon,klev))
     701      ALLOCATE(xt_ancien(ntraciso,klon,klev))
     702      ALLOCATE(xtl_ancien(ntraciso,klon,klev))
     703      ALLOCATE(xts_ancien(ntraciso,klon,klev))
     704      ALLOCATE(xtrain_fall(ntraciso,klon))
     705      ALLOCATE(xtsnow_fall(ntraciso,klon))
     706      ALLOCATE(xtrain_con(ntraciso,klon))
     707      ALLOCATE(xtsnow_con(ntraciso,klon))
     708#ifdef ISOTRAC
     709      ALLOCATE(bassin_map(klon))
     710      ALLOCATE(boite_map(klon,klev)) 
     711#endif     
     712#endif
     713
    661714      ALLOCATE(ale_bl_trig(klon))
    662715      ALLOCATE(ratqs_inter(klon,klev))
     
    811864           delta_sst)
    812865
     866#ifdef ISO   
     867      DEALLOCATE(xtsol,fxtevap) 
     868      DEALLOCATE(xt_ancien,xtl_ancien,xts_ancien, fxtd, wake_deltaxt)
     869      DEALLOCATE(xtrain_fall, xtsnow_fall, xtrain_con, xtsnow_con)
     870#ifdef ISOTRAC
     871      DEALLOCATE(bassin_map,boite_map)
     872#endif       
     873#endif
    813874      is_initialized=.FALSE.
    814875     
  • LMDZ6/trunk/libf/phylmdiso/phys_state_var_mod.F90

    r4040 r4088  
    100100      REAL, ALLOCATABLE, SAVE :: clwcon(:,:),rnebcon(:,:)
    101101!$OMP THREADPRIVATE(clwcon,rnebcon)
     102      REAL, ALLOCATABLE, SAVE :: rneb_ancien(:,:)
     103!$OMP THREADPRIVATE(rneb_ancien)
    102104      REAL, ALLOCATABLE, SAVE :: qtc_cv(:,:),sigt_cv(:,:)
    103105!$OMP THREADPRIVATE(qtc_cv,sigt_cv)
     
    549551!!! Rom P <<<
    550552      ALLOCATE(clwcon(klon,klev),rnebcon(klon,klev))
     553      ALLOCATE(rneb_ancien(klon,klev))
    551554      ALLOCATE(qtc_cv(klon,klev),sigt_cv(klon,klev))
    552555      ALLOCATE(ratqs(klon,klev))
     
    742745      DEALLOCATE(zthe, zpic, zval)
    743746      DEALLOCATE(rugoro, t_ancien, q_ancien, clwcon, rnebcon)
    744       DEALLOCATE(qs_ancien, ql_ancien)
     747      DEALLOCATE(qs_ancien, ql_ancien, rneb_ancien)
    745748      DEALLOCATE(prw_ancien, prlw_ancien, prsw_ancien)
    746749      DEALLOCATE(qtc_cv,sigt_cv)
Note: See TracChangeset for help on using the changeset viewer.