Ignore:
Timestamp:
Apr 14, 2017, 4:42:31 PM (7 years ago)
Author:
oboucher
Message:

Introducing dry AOD diagnostics for the total aerosols and specieswise
The calculations are only performed if the diagnostics are requested

File:
1 edited

Legend:

Unmodified
Added
Removed

  • LMDZ5/trunk/libf/phylmd/phys_local_var_mod.F90

    r2849 r2854  
    141141      REAL, SAVE, ALLOCATABLE :: tausum_aero(:,:,:)
    142142      !$OMP THREADPRIVATE(tausum_aero)
     143      REAL, SAVE, ALLOCATABLE :: drytausum_aero(:,:)
     144      !$OMP THREADPRIVATE(drytausum_aero)
    143145      REAL, SAVE, ALLOCATABLE :: tau3d_aero(:,:,:,:)
    144146      !$OMP THREADPRIVATE(tau3d_aero)
    … …  
    167169      REAL, SAVE, ALLOCATABLE :: od550aer(:)
    168170      !$OMP THREADPRIVATE(od550aer)
     171      REAL, SAVE, ALLOCATABLE :: dryod550aer(:)
     172      !$OMP THREADPRIVATE(dryod550aer)
    169173      REAL, SAVE, ALLOCATABLE :: abs550aer(:)
    170174      !$OMP THREADPRIVATE(abs550aer)
    … …  
    500504
    501505IMPLICIT NONE
    502       allocate(t_seri(klon,klev),q_seri(klon,klev),ql_seri(klon,klev),qs_seri(klon,klev))
    503       allocate(u_seri(klon,klev),v_seri(klon,klev))
    504       allocate(l_mixmin(klon,klev,nbsrf), l_mix(klon,klev,nbsrf))
     506      ALLOCATE(t_seri(klon,klev),q_seri(klon,klev),ql_seri(klon,klev),qs_seri(klon,klev))
     507      ALLOCATE(u_seri(klon,klev),v_seri(klon,klev))
     508      ALLOCATE(l_mixmin(klon,klev,nbsrf), l_mix(klon,klev,nbsrf))
    505509      l_mix(:,:,:)=0. ; l_mixmin(:,:,:)=0. ! doit etre initialse car pas toujours remplis
    506510
    507       allocate(tr_seri(klon,klev,nbtr))
    508       allocate(d_t_dyn(klon,klev),d_q_dyn(klon,klev))
    509       allocate(d_ql_dyn(klon,klev),d_qs_dyn(klon,klev))
    510       allocate(d_q_dyn2d(klon),d_ql_dyn2d(klon),d_qs_dyn2d(klon))
    511       allocate(d_u_dyn(klon,klev),d_v_dyn(klon,klev))
    512       allocate(d_tr_dyn(klon,klev,nbtr))                   !RomP
    513       allocate(d_t_con(klon,klev),d_q_con(klon,klev))
    514       allocate(d_u_con(klon,klev),d_v_con(klon,klev))
    515       allocate(d_t_wake(klon,klev),d_q_wake(klon,klev))
    516       allocate(d_t_lsc(klon,klev),d_q_lsc(klon,klev))
    517       allocate(d_t_lwr(klon,klev),d_t_lw0(klon,klev))
    518       allocate(d_t_swr(klon,klev),d_t_sw0(klon,klev))
    519       allocate(d_ql_lsc(klon,klev),d_qi_lsc(klon,klev))
    520       allocate(d_t_ajsb(klon,klev),d_q_ajsb(klon,klev))
    521       allocate(d_t_ajs(klon,klev),d_q_ajs(klon,klev))
     511      ALLOCATE(tr_seri(klon,klev,nbtr))
     512      ALLOCATE(d_t_dyn(klon,klev),d_q_dyn(klon,klev))
     513      ALLOCATE(d_ql_dyn(klon,klev),d_qs_dyn(klon,klev))
     514      ALLOCATE(d_q_dyn2d(klon),d_ql_dyn2d(klon),d_qs_dyn2d(klon))
     515      ALLOCATE(d_u_dyn(klon,klev),d_v_dyn(klon,klev))
     516      ALLOCATE(d_tr_dyn(klon,klev,nbtr))                   !RomP
     517      ALLOCATE(d_t_con(klon,klev),d_q_con(klon,klev))
     518      ALLOCATE(d_u_con(klon,klev),d_v_con(klon,klev))
     519      ALLOCATE(d_t_wake(klon,klev),d_q_wake(klon,klev))
     520      ALLOCATE(d_t_lsc(klon,klev),d_q_lsc(klon,klev))
     521      ALLOCATE(d_t_lwr(klon,klev),d_t_lw0(klon,klev))
     522      ALLOCATE(d_t_swr(klon,klev),d_t_sw0(klon,klev))
     523      ALLOCATE(d_ql_lsc(klon,klev),d_qi_lsc(klon,klev))
     524      ALLOCATE(d_t_ajsb(klon,klev),d_q_ajsb(klon,klev))
     525      ALLOCATE(d_t_ajs(klon,klev),d_q_ajs(klon,klev))
    522526!nrlmd<
    523       allocate(d_t_ajs_w(klon,klev),d_q_ajs_w(klon,klev))
    524       allocate(d_t_ajs_x(klon,klev),d_q_ajs_x(klon,klev))
     527      ALLOCATE(d_t_ajs_w(klon,klev),d_q_ajs_w(klon,klev))
     528      ALLOCATE(d_t_ajs_x(klon,klev),d_q_ajs_x(klon,klev))
    525529!>nrlmd
    526       allocate(d_u_ajs(klon,klev),d_v_ajs(klon,klev))
    527       allocate(d_t_eva(klon,klev),d_q_eva(klon,klev))
    528       allocate(d_ql_eva(klon,klev),d_qi_eva(klon,klev))
    529       allocate(d_t_lscst(klon,klev),d_q_lscst(klon,klev))
    530       allocate(d_t_lscth(klon,klev),d_q_lscth(klon,klev))
    531       allocate(plul_st(klon),plul_th(klon))
    532       allocate(d_t_vdf(klon,klev),d_q_vdf(klon,klev),d_t_diss(klon,klev))
     530      ALLOCATE(d_u_ajs(klon,klev),d_v_ajs(klon,klev))
     531      ALLOCATE(d_t_eva(klon,klev),d_q_eva(klon,klev))
     532      ALLOCATE(d_ql_eva(klon,klev),d_qi_eva(klon,klev))
     533      ALLOCATE(d_t_lscst(klon,klev),d_q_lscst(klon,klev))
     534      ALLOCATE(d_t_lscth(klon,klev),d_q_lscth(klon,klev))
     535      ALLOCATE(plul_st(klon),plul_th(klon))
     536      ALLOCATE(d_t_vdf(klon,klev),d_q_vdf(klon,klev),d_t_diss(klon,klev))
    533537!nrlmd+jyg<
    534       allocate(d_t_vdf_w(klon,klev),d_q_vdf_w(klon,klev))
    535       allocate(d_t_vdf_x(klon,klev),d_q_vdf_x(klon,klev))
     538      ALLOCATE(d_t_vdf_w(klon,klev),d_q_vdf_w(klon,klev))
     539      ALLOCATE(d_t_vdf_x(klon,klev),d_q_vdf_x(klon,klev))
    536540!>nrlmd+jyg
    537       allocate(d_u_vdf(klon,klev),d_v_vdf(klon,klev))
    538       allocate(d_t_oli(klon,klev),d_t_oro(klon,klev))
    539       allocate(d_u_oli(klon,klev),d_v_oli(klon,klev))
    540       allocate(d_u_oro(klon,klev),d_v_oro(klon,klev))
    541       allocate(d_t_lif(klon,klev),d_t_ec(klon,klev))
    542       allocate(d_u_lif(klon,klev),d_v_lif(klon,klev))
    543       allocate(d_ts(klon,nbsrf), d_tr(klon,klev,nbtr))
     541      ALLOCATE(d_u_vdf(klon,klev),d_v_vdf(klon,klev))
     542      ALLOCATE(d_t_oli(klon,klev),d_t_oro(klon,klev))
     543      ALLOCATE(d_u_oli(klon,klev),d_v_oli(klon,klev))
     544      ALLOCATE(d_u_oro(klon,klev),d_v_oro(klon,klev))
     545      ALLOCATE(d_t_lif(klon,klev),d_t_ec(klon,klev))
     546      ALLOCATE(d_u_lif(klon,klev),d_v_lif(klon,klev))
     547      ALLOCATE(d_ts(klon,nbsrf), d_tr(klon,klev,nbtr))
    544548! Special RRTM
    545       allocate(ZLWFT0_i(klon,klev+1),ZSWFT0_i(klon,klev+1),ZFLDN0(klon,klev+1))
    546       allocate(ZFLUP0(klon,klev+1),ZFSDN0(klon,klev+1),ZFSUP0(klon,klev+1))
    547 !
    548       allocate(topswad_aero(klon), solswad_aero(klon))
    549       allocate(topswai_aero(klon), solswai_aero(klon))
    550       allocate(topswad0_aero(klon), solswad0_aero(klon))
    551      ! LW diagnostics CK
    552       allocate(toplwad_aero(klon), sollwad_aero(klon))
    553       allocate(toplwai_aero(klon), sollwai_aero(klon))
    554       allocate(toplwad0_aero(klon), sollwad0_aero(klon))
    555       ! end
    556       allocate(topsw_aero(klon,naero_grp), solsw_aero(klon,naero_grp))
    557       allocate(topsw0_aero(klon,naero_grp), solsw0_aero(klon,naero_grp))
    558       allocate(topswcf_aero(klon,3), solswcf_aero(klon,3))
    559       allocate(du_gwd_hines(klon,klev),dv_gwd_hines(klon,klev))
    560       allocate(dv_gwd_rando(klon,klev),dv_gwd_front(klon,klev))
    561       allocate(east_gwstress(klon,klev),west_gwstress(klon,klev))
    562       allocate(d_t_hin(klon,klev))
    563       allocate(d_q_ch4(klon,klev))
    564 !      allocate(tausum_aero(klon,nwave,naero_spc))
    565 !      allocate(tau3d_aero(klon,klev,nwave,naero_spc))
    566       allocate(stratomask(klon,klev))
    567 !--correction mini bug OB
    568       allocate(tausum_aero(klon,nwave,naero_tot))
    569       allocate(tau3d_aero(klon,klev,nwave,naero_tot))
    570       allocate(scdnc(klon, klev))
    571       allocate(cldncl(klon))
    572       allocate(reffclwtop(klon))
    573       allocate(lcc(klon))
    574       allocate(reffclws(klon, klev))
    575       allocate(reffclwc(klon, klev))
    576       allocate(cldnvi(klon))
    577       allocate(lcc3d(klon, klev))
    578       allocate(lcc3dcon(klon, klev))
    579       allocate(lcc3dstra(klon, klev))
    580       allocate(od443aer(klon))
    581       allocate(od550aer(klon))
    582       allocate(od865aer(klon))
    583       allocate(abs550aer(klon))
    584       allocate(ec550aer(klon,klev))
    585       allocate(od550lt1aer(klon))
    586       allocate(sconcso4(klon))
    587       allocate(sconcno3(klon))
    588       allocate(sconcoa(klon))
    589       allocate(sconcbc(klon))
    590       allocate(sconcss(klon))
    591       allocate(sconcdust(klon))
    592       allocate(concso4(klon,klev))
    593       allocate(concno3(klon,klev))
    594       allocate(concoa(klon,klev))
    595       allocate(concbc(klon,klev))
    596       allocate(concss(klon,klev))
    597       allocate(concdust(klon,klev))
    598       allocate(loadso4(klon))
    599       allocate(loadoa(klon))
    600       allocate(loadbc(klon))
    601       allocate(loadss(klon))
    602       allocate(loaddust(klon))
    603       allocate(loadno3(klon))
    604       allocate(load_tmp1(klon))
    605       allocate(load_tmp2(klon))
    606       allocate(load_tmp3(klon))
     549      ALLOCATE(ZLWFT0_i(klon,klev+1),ZSWFT0_i(klon,klev+1),ZFLDN0(klon,klev+1))
     550      ALLOCATE(ZFLUP0(klon,klev+1),ZFSDN0(klon,klev+1),ZFSUP0(klon,klev+1))
     551!
     552      ALLOCATE(topswad_aero(klon), solswad_aero(klon))
     553      ALLOCATE(topswai_aero(klon), solswai_aero(klon))
     554      ALLOCATE(topswad0_aero(klon), solswad0_aero(klon))
     555      ALLOCATE(toplwad_aero(klon), sollwad_aero(klon))
     556      ALLOCATE(toplwai_aero(klon), sollwai_aero(klon))
     557      ALLOCATE(toplwad0_aero(klon), sollwad0_aero(klon))
     558      ALLOCATE(topsw_aero(klon,naero_grp), solsw_aero(klon,naero_grp))
     559      ALLOCATE(topsw0_aero(klon,naero_grp), solsw0_aero(klon,naero_grp))
     560      ALLOCATE(topswcf_aero(klon,3), solswcf_aero(klon,3))
     561      ALLOCATE(du_gwd_hines(klon,klev),dv_gwd_hines(klon,klev))
     562      ALLOCATE(dv_gwd_rando(klon,klev),dv_gwd_front(klon,klev))
     563      ALLOCATE(east_gwstress(klon,klev),west_gwstress(klon,klev))
     564      ALLOCATE(d_t_hin(klon,klev))
     565      ALLOCATE(d_q_ch4(klon,klev))
     566      ALLOCATE(stratomask(klon,klev))
     567      ALLOCATE(tausum_aero(klon,nwave,naero_tot))
     568      ALLOCATE(drytausum_aero(klon,naero_tot))
     569      ALLOCATE(tau3d_aero(klon,klev,nwave,naero_tot))
     570      ALLOCATE(scdnc(klon, klev))
     571      ALLOCATE(cldncl(klon))
     572      ALLOCATE(reffclwtop(klon))
     573      ALLOCATE(lcc(klon))
     574      ALLOCATE(reffclws(klon, klev))
     575      ALLOCATE(reffclwc(klon, klev))
     576      ALLOCATE(cldnvi(klon))
     577      ALLOCATE(lcc3d(klon, klev))
     578      ALLOCATE(lcc3dcon(klon, klev))
     579      ALLOCATE(lcc3dstra(klon, klev))
     580      ALLOCATE(od443aer(klon))
     581      ALLOCATE(od550aer(klon))
     582      ALLOCATE(od865aer(klon))
     583      ALLOCATE(dryod550aer(klon))
     584      ALLOCATE(abs550aer(klon))
     585      ALLOCATE(ec550aer(klon,klev))
     586      ALLOCATE(od550lt1aer(klon))
     587      ALLOCATE(sconcso4(klon))
     588      ALLOCATE(sconcno3(klon))
     589      ALLOCATE(sconcoa(klon))
     590      ALLOCATE(sconcbc(klon))
     591      ALLOCATE(sconcss(klon))
     592      ALLOCATE(sconcdust(klon))
     593      ALLOCATE(concso4(klon,klev))
     594      ALLOCATE(concno3(klon,klev))
     595      ALLOCATE(concoa(klon,klev))
     596      ALLOCATE(concbc(klon,klev))
     597      ALLOCATE(concss(klon,klev))
     598      ALLOCATE(concdust(klon,klev))
     599      ALLOCATE(loadso4(klon))
     600      ALLOCATE(loadoa(klon))
     601      ALLOCATE(loadbc(klon))
     602      ALLOCATE(loadss(klon))
     603      ALLOCATE(loaddust(klon))
     604      ALLOCATE(loadno3(klon))
     605      ALLOCATE(load_tmp1(klon))
     606      ALLOCATE(load_tmp2(klon))
     607      ALLOCATE(load_tmp3(klon))
    607608
    608609!IM ajout variables CFMIP2/CMIP5
    … …  
    772773USE indice_sol_mod
    773774IMPLICIT NONE
    774       deallocate(t_seri,q_seri,ql_seri,qs_seri)
    775       deallocate(u_seri,v_seri)
    776       deallocate(l_mixmin,l_mix)
    777 
    778       deallocate(tr_seri)
    779       deallocate(d_t_dyn,d_q_dyn)
    780       deallocate(d_ql_dyn,d_qs_dyn)
    781       deallocate(d_q_dyn2d,d_ql_dyn2d,d_qs_dyn2d)
    782       deallocate(d_u_dyn,d_v_dyn)
    783       deallocate(d_tr_dyn)                      !RomP
    784       deallocate(d_t_con,d_q_con)
    785       deallocate(d_u_con,d_v_con)
    786       deallocate(d_t_wake,d_q_wake)
    787       deallocate(d_t_lsc,d_q_lsc)
    788       deallocate(d_t_lwr,d_t_lw0)
    789       deallocate(d_t_swr,d_t_sw0)
    790       deallocate(d_ql_lsc,d_qi_lsc)
    791       deallocate(d_t_ajsb,d_q_ajsb)
    792       deallocate(d_t_ajs,d_q_ajs)
     775      DEALLOCATE(t_seri,q_seri,ql_seri,qs_seri)
     776      DEALLOCATE(u_seri,v_seri)
     777      DEALLOCATE(l_mixmin,l_mix)
     778
     779      DEALLOCATE(tr_seri)
     780      DEALLOCATE(d_t_dyn,d_q_dyn)
     781      DEALLOCATE(d_ql_dyn,d_qs_dyn)
     782      DEALLOCATE(d_q_dyn2d,d_ql_dyn2d,d_qs_dyn2d)
     783      DEALLOCATE(d_u_dyn,d_v_dyn)
     784      DEALLOCATE(d_tr_dyn)                      !RomP
     785      DEALLOCATE(d_t_con,d_q_con)
     786      DEALLOCATE(d_u_con,d_v_con)
     787      DEALLOCATE(d_t_wake,d_q_wake)
     788      DEALLOCATE(d_t_lsc,d_q_lsc)
     789      DEALLOCATE(d_t_lwr,d_t_lw0)
     790      DEALLOCATE(d_t_swr,d_t_sw0)
     791      DEALLOCATE(d_ql_lsc,d_qi_lsc)
     792      DEALLOCATE(d_t_ajsb,d_q_ajsb)
     793      DEALLOCATE(d_t_ajs,d_q_ajs)
    793794!nrlmd<
    794       deallocate(d_t_ajs_w,d_q_ajs_w)
    795       deallocate(d_t_ajs_x,d_q_ajs_x)
     795      DEALLOCATE(d_t_ajs_w,d_q_ajs_w)
     796      DEALLOCATE(d_t_ajs_x,d_q_ajs_x)
    796797!>nrlmd
    797       deallocate(d_u_ajs,d_v_ajs)
    798       deallocate(d_t_eva,d_q_eva)
    799       deallocate(d_ql_eva,d_qi_eva)
    800       deallocate(d_t_lscst,d_q_lscst)
    801       deallocate(d_t_lscth,d_q_lscth)
    802       deallocate(plul_st,plul_th)
    803       deallocate(d_t_vdf,d_q_vdf,d_t_diss)
     798      DEALLOCATE(d_u_ajs,d_v_ajs)
     799      DEALLOCATE(d_t_eva,d_q_eva)
     800      DEALLOCATE(d_ql_eva,d_qi_eva)
     801      DEALLOCATE(d_t_lscst,d_q_lscst)
     802      DEALLOCATE(d_t_lscth,d_q_lscth)
     803      DEALLOCATE(plul_st,plul_th)
     804      DEALLOCATE(d_t_vdf,d_q_vdf,d_t_diss)
    804805!nrlmd+jyg<
    805       deallocate(d_t_vdf_w,d_q_vdf_w)
    806       deallocate(d_t_vdf_x,d_q_vdf_x)
     806      DEALLOCATE(d_t_vdf_w,d_q_vdf_w)
     807      DEALLOCATE(d_t_vdf_x,d_q_vdf_x)
    807808!>nrlmd+jyg
    808       deallocate(d_u_vdf,d_v_vdf)
    809       deallocate(d_t_oli,d_t_oro)
    810       deallocate(d_u_oli,d_v_oli)
    811       deallocate(d_u_oro,d_v_oro)
    812       deallocate(d_t_lif,d_t_ec)
    813       deallocate(d_u_lif,d_v_lif)
    814       deallocate(d_ts, d_tr)
    815       deallocate(topswad_aero,solswad_aero)
    816       deallocate(topswai_aero,solswai_aero)
    817       deallocate(topswad0_aero,solswad0_aero)
    818       ! LW additional CK
    819       deallocate(toplwad_aero,sollwad_aero)
    820       deallocate(toplwai_aero,sollwai_aero)
    821       deallocate(toplwad0_aero,sollwad0_aero)
    822       ! end
    823       deallocate(topsw_aero,solsw_aero)
    824       deallocate(topsw0_aero,solsw0_aero)
    825       deallocate(topswcf_aero,solswcf_aero)
    826       deallocate(stratomask)
    827       deallocate(tausum_aero)
    828       deallocate(tau3d_aero)
    829       deallocate(scdnc)
    830       deallocate(cldncl)
    831       deallocate(reffclwtop)
    832       deallocate(lcc)
    833       deallocate(reffclws)
    834       deallocate(reffclwc)
    835       deallocate(cldnvi)
    836       deallocate(lcc3d)
    837       deallocate(lcc3dcon)
    838       deallocate(lcc3dstra)
    839       deallocate(od443aer)
    840       deallocate(od550aer)
    841       deallocate(od865aer)
    842       deallocate(abs550aer)
    843       deallocate(ec550aer)
    844       deallocate(od550lt1aer)
    845       deallocate(sconcso4)
    846       deallocate(sconcno3)
    847       deallocate(sconcoa)
    848       deallocate(sconcbc)
    849       deallocate(sconcss)
    850       deallocate(sconcdust)
    851       deallocate(concso4)
    852       deallocate(concno3)
    853       deallocate(concoa)
    854       deallocate(concbc)
    855       deallocate(concss)
    856       deallocate(concdust)
    857       deallocate(loadso4)
    858       deallocate(loadoa)
    859       deallocate(loadbc)
    860       deallocate(loadss)
    861       deallocate(loaddust)
    862       deallocate(loadno3)
    863       deallocate(load_tmp1)
    864       deallocate(load_tmp2)
    865       deallocate(load_tmp3)
    866       deallocate(du_gwd_hines,dv_gwd_hines,d_t_hin)
    867       deallocate(d_q_ch4)
    868       deallocate(dv_gwd_rando,dv_gwd_front)
    869       deallocate(east_gwstress,west_gwstress)
     809      DEALLOCATE(d_u_vdf,d_v_vdf)
     810      DEALLOCATE(d_t_oli,d_t_oro)
     811      DEALLOCATE(d_u_oli,d_v_oli)
     812      DEALLOCATE(d_u_oro,d_v_oro)
     813      DEALLOCATE(d_t_lif,d_t_ec)
     814      DEALLOCATE(d_u_lif,d_v_lif)
     815      DEALLOCATE(d_ts, d_tr)
     816      DEALLOCATE(topswad_aero,solswad_aero)
     817      DEALLOCATE(topswai_aero,solswai_aero)
     818      DEALLOCATE(topswad0_aero,solswad0_aero)
     819      DEALLOCATE(toplwad_aero,sollwad_aero)
     820      DEALLOCATE(toplwai_aero,sollwai_aero)
     821      DEALLOCATE(toplwad0_aero,sollwad0_aero)
     822      DEALLOCATE(topsw_aero,solsw_aero)
     823      DEALLOCATE(topsw0_aero,solsw0_aero)
     824      DEALLOCATE(topswcf_aero,solswcf_aero)
     825      DEALLOCATE(stratomask)
     826      DEALLOCATE(tausum_aero)
     827      DEALLOCATE(drytausum_aero)
     828      DEALLOCATE(tau3d_aero)
     829      DEALLOCATE(scdnc)
     830      DEALLOCATE(cldncl)
     831      DEALLOCATE(reffclwtop)
     832      DEALLOCATE(lcc)
     833      DEALLOCATE(reffclws)
     834      DEALLOCATE(reffclwc)
     835      DEALLOCATE(cldnvi)
     836      DEALLOCATE(lcc3d)
     837      DEALLOCATE(lcc3dcon)
     838      DEALLOCATE(lcc3dstra)
     839      DEALLOCATE(od443aer)
     840      DEALLOCATE(od550aer)
     841      DEALLOCATE(od865aer)
     842      DEALLOCATE(dryod550aer)
     843      DEALLOCATE(abs550aer)
     844      DEALLOCATE(ec550aer)
     845      DEALLOCATE(od550lt1aer)
     846      DEALLOCATE(sconcso4)
     847      DEALLOCATE(sconcno3)
     848      DEALLOCATE(sconcoa)
     849      DEALLOCATE(sconcbc)
     850      DEALLOCATE(sconcss)
     851      DEALLOCATE(sconcdust)
     852      DEALLOCATE(concso4)
     853      DEALLOCATE(concno3)
     854      DEALLOCATE(concoa)
     855      DEALLOCATE(concbc)
     856      DEALLOCATE(concss)
     857      DEALLOCATE(concdust)
     858      DEALLOCATE(loadso4)
     859      DEALLOCATE(loadoa)
     860      DEALLOCATE(loadbc)
     861      DEALLOCATE(loadss)
     862      DEALLOCATE(loaddust)
     863      DEALLOCATE(loadno3)
     864      DEALLOCATE(load_tmp1)
     865      DEALLOCATE(load_tmp2)
     866      DEALLOCATE(load_tmp3)
     867      DEALLOCATE(du_gwd_hines,dv_gwd_hines,d_t_hin)
     868      DEALLOCATE(d_q_ch4)
     869      DEALLOCATE(dv_gwd_rando,dv_gwd_front)
     870      DEALLOCATE(east_gwstress,west_gwstress)
    870871
    871872!IM ajout variables CFMIP2/CMIP5
    872       deallocate(topswad_aerop, solswad_aerop)
    873       deallocate(topswai_aerop, solswai_aerop)
    874       deallocate(topswad0_aerop, solswad0_aerop)
    875       deallocate(topsw_aerop, topsw0_aerop)
    876       deallocate(solsw_aerop, solsw0_aerop)
    877       deallocate(topswcf_aerop, solswcf_aerop)
     873      DEALLOCATE(topswad_aerop, solswad_aerop)
     874      DEALLOCATE(topswai_aerop, solswai_aerop)
     875      DEALLOCATE(topswad0_aerop, solswad0_aerop)
     876      DEALLOCATE(topsw_aerop, topsw0_aerop)
     877      DEALLOCATE(solsw_aerop, solsw0_aerop)
     878      DEALLOCATE(topswcf_aerop, solswcf_aerop)
    878879
    879880!CK LW diagnostics
    880       deallocate(toplwad_aerop, sollwad_aerop)
    881       deallocate(toplwai_aerop, sollwai_aerop)
    882       deallocate(toplwad0_aerop, sollwad0_aerop)
     881      DEALLOCATE(toplwad_aerop, sollwad_aerop)
     882      DEALLOCATE(toplwai_aerop, sollwai_aerop)
     883      DEALLOCATE(toplwad0_aerop, sollwad0_aerop)
    883884
    884885! FH Ajout de celles nécessaires au phys_output_write_mod
Note: See TracChangeset for help on using the changeset viewer.