Changeset 1877 for LMDZ5/trunk


Ignore:
Timestamp:
Sep 25, 2013, 3:35:04 PM (11 years ago)
Author:
Laurent Fairhead
Message:

Ménage sur le code pour éliminer les calculs spécifiques ISCCP hors COSP
Un certain nombre de variables non utilisées dans physiq.F90 ont aussi
été supprimée. Environ 840 lignes supprimées du code physiq.F90

Code cleanup to eliminate specific references to ISCCP outside the COSP
library. Unused variables have been cleaned up from the physiq.F90 routine
as well.

Location:
LMDZ5/trunk/libf/phylmd
Files:
5 deleted
6 edited

Legend:

Unmodified
Added
Removed

  • LMDZ5/trunk/libf/phylmd/clesphys.h

    r1849 r1877  
    6666! aer_type: pour utiliser un fichier constant dans readaerosol
    6767       CHARACTER*8 :: aer_type
    68        LOGICAL ok_isccp, ok_regdyn
     68       LOGICAL ok_regdyn
    6969       REAL lonmin_ins, lonmax_ins, latmin_ins, latmax_ins
    7070       REAL ecrit_ins, ecrit_hf, ecrit_day
    … …  
    9494     &     , lev_histins, lev_histLES, lev_histdayNMC, levout_histNMC   &
    9595     &     , pasphys, ok_histNMC, freq_outNMC, freq_calNMC              &
    96      &     , type_run, ok_isccp, ok_regdyn, ok_cosp                     &
     96     &     , type_run, ok_regdyn, ok_cosp                               &
    9797     &     , ok_mensuelCOSP,ok_journeCOSP,ok_hfCOSP                     &
    9898     &     , lonmin_ins, lonmax_ins, latmin_ins, latmax_ins             &

  • LMDZ5/trunk/libf/phylmd/conema3.h

    r879 r1877  
    66      logical ok_adj_ema      ! F
    77      integer iflag_clw      ! 0
    8           integer iflag_cvl_sigd
     8      integer iflag_cvl_sigd
    99      real sig1feed      ! 1.
    1010      real sig2feed      ! 0.95

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

    r1865 r1877  
    161161  REAL, SAVE :: freq_outNMC_omp(3), freq_calNMC_omp(3)
    162162  CHARACTER*4, SAVE :: type_run_omp
    163   LOGICAL,SAVE :: ok_isccp_omp
    164163  LOGICAL,SAVE :: ok_cosp_omp
    165164  LOGICAL,SAVE :: ok_mensuelCOSP_omp,ok_journeCOSP_omp,ok_hfCOSP_omp
    … …  
    13511350  type_run_omp = 'AMIP'
    13521351  call getin('type_run',type_run_omp)
    1353 
    1354 !
    1355 !Config Key  = ok_isccp
    1356 !Config Desc =
    1357 !Config Def  = .true.
    1358 !Config Help =
    1359 !
    1360 ! ok_isccp = .true.
    1361   ok_isccp_omp = .false.
    1362   call getin('ok_isccp',ok_isccp_omp)
    13631352
    13641353!
    … …  
    17991788    iflag_cvl_sigd = iflag_cvl_sigd_omp
    18001789    type_run = type_run_omp
    1801     ok_isccp = ok_isccp_omp
    18021790    ok_cosp = ok_cosp_omp
    18031791    ok_mensuelCOSP = ok_mensuelCOSP_omp
    … …  
    19721960  write(lunout,*)' iflag_clos = ', iflag_clos
    19731961  write(lunout,*)' type_run = ',type_run
    1974   write(lunout,*)' ok_isccp = ',ok_isccp
    19751962  write(lunout,*)' ok_cosp = ',ok_cosp
    19761963  write(lunout,*)' ok_mensuelCOSP = ',ok_mensuelCOSP

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

    r1865 r1877  
    6969  !!   =====================
    7070#define histNMC
    71   !#define histISCCP
    7271  !!======================================================================
    7372  !!    modif   ( P. Le Van ,  12/10/98 )
    … …  
    184183  real facteur
    185184
    186   REAL zz,znum,zden
    187185  REAL wmax_th(klon)
    188186  REAL tau_overturning_th(klon)
    … …  
    260258  !
    261259  REAL PVteta(klon,nbteta)
    262   REAL zx_tmp_3dte(iim,jjmp1,nbteta)
    263260  !
    264261  !MI Amip2 PV a theta constante
    … …  
    274271  include "declare_STDlev.h"
    275272  !
    276   CHARACTER*4 bb2
    277   CHARACTER*2 bb3
    278273  !
    279274  include "radopt.h"
    … …  
    281276  !
    282277
    283   REAL convliq(klon,klev)  ! eau liquide nuageuse convective
    284   REAL convfra(klon,klev)  ! fraction nuageuse convective
    285 
    286   REAL cldl_c(klon),cldm_c(klon),cldh_c(klon) !nuages bas, moyen et haut
    287   REAL cldt_c(klon),cldq_c(klon) !nuage total, eau liquide integree
    288   REAL cldl_s(klon),cldm_s(klon),cldh_s(klon) !nuages bas, moyen et haut
    289   REAL cldt_s(klon),cldq_s(klon) !nuage total, eau liquide integree
    290 
    291   INTEGER linv, kp1
    292   ! flwp, fiwp = Liquid Water Path & Ice Water Path (kg/m2)
    293   ! flwc, fiwc = Liquid Water Content & Ice Water Content (kg/kg)
    294   REAL flwp_c(klon), fiwp_c(klon)
    295   REAL flwc_c(klon,klev), fiwc_c(klon,klev)
    296   REAL flwp_s(klon), fiwp_s(klon)
    297   REAL flwc_s(klon,klev), fiwc_s(klon,klev)
    298 
    299 
    300   !IM ISCCP simulator v3.4
    301   ! dans clesphys.h top_height, overlap
    302   !v3.4
    303   INTEGER debug, debugcol
     278
     279  INTEGER debug
     280  INTEGER n
    304281  !ym      INTEGER npoints
    305282  !ym      PARAMETER(npoints=klon)
    306283  !
    307   INTEGER sunlit(klon) !sunlit=1 if day; sunlit=0 if night
    308284  INTEGER nregISCtot
    309285  PARAMETER(nregISCtot=1)
    … …  
    315291  INTEGER imin_debut, nbpti
    316292  INTEGER jmin_debut, nbptj
    317   !IM parametres ISCCP BEG
    318   INTEGER nbapp_isccp
    319   !     INTEGER nbapp_isccp,isccppas
    320   !     PARAMETER(isccppas=6) !appel du simulateurs tous les 6pas de temps de la physique
    321   !                           !i.e. toutes les 3 heures
    322   INTEGER n
    323   INTEGER ifreq_isccp(napisccp), freqin_pdt(napisccp)
    324   DATA ifreq_isccp/3/
    325   SAVE ifreq_isccp
    326   !$OMP THREADPRIVATE(ifreq_isccp)
    327   CHARACTER*5 typinout(napisccp)
    328   DATA typinout/'i3od'/
    329   SAVE typinout
    330   !$OMP THREADPRIVATE(typinout)
    331   !IM verif boxptop BEG
    332   CHARACTER*1 verticaxe(napisccp)
    333   DATA verticaxe/'1'/
    334   SAVE verticaxe
    335   !$OMP THREADPRIVATE(verticaxe)
    336   !IM verif boxptop END
    337   INTEGER nvlev(napisccp)
    338   !     INTEGER nvlev
    339   REAL t1, aa
    340   REAL seed_re(klon,napisccp)
    341   !ym !!!! A voir plus tard
    342   !ym      INTEGER iphy(iim,jjmp1)
    343   !IM parametres ISCCP END
    344   !
    345   ! ncol = nb. de sous-colonnes pour chaque maille du GCM
    346   ! ncolmx = No. max. de sous-colonnes pour chaque maille du GCM
    347   !      INTEGER ncol(napisccp), ncolmx, seed(klon,napisccp)
    348   INTEGER,SAVE :: ncol(napisccp)
    349   !$OMP THREADPRIVATE(ncol)
    350   INTEGER ncolmx, seed(klon,napisccp)
    351   REAL nbsunlit(nregISCtot,klon,napisccp)  !nbsunlit : moyenne de sunlit
    352   !     PARAMETER(ncolmx=1500)
    353   PARAMETER(ncolmx=300)
    354   !
    355   !IM verif boxptop BEG
    356   REAL vertlev(ncolmx,napisccp)
    357   !IM verif boxptop END
    358   !
    359   REAL,SAVE :: tautab_omp(0:255),tautab(0:255)
    360   INTEGER,SAVE :: invtau_omp(-20:45000),invtau(-20:45000)
    361   !$OMP THREADPRIVATE(tautab,invtau)
    362   REAL emsfc_lw
    363   PARAMETER(emsfc_lw=0.99)
    364   !     REAL    ran0                      ! type for random number fuction
    365   !
    366   REAL cldtot(klon,klev)
    367   ! variables de haut en bas pour le simulateur ISCCP
    368   REAL dtau_s(klon,klev) !tau nuages startiformes
    369   REAL dtau_c(klon,klev) !tau nuages convectifs
    370   REAL dem_s(klon,klev)  !emissivite nuages startiformes
    371   REAL dem_c(klon,klev)  !emissivite nuages convectifs
    372   !
    373   ! variables de haut en bas pour le simulateur ISCCP
    374   REAL pfull(klon,klev)
    375   REAL phalf(klon,klev+1)
    376   REAL qv(klon,klev)
    377   REAL cc(klon,klev)
    378   REAL conv(klon,klev)
    379   REAL dtau_sH2B(klon,klev)
    380   REAL dtau_cH2B(klon,klev)
    381   REAL at(klon,klev)
    382   REAL dem_sH2B(klon,klev)
    383   REAL dem_cH2B(klon,klev)
    384   !
    385   INTEGER kmax, lmax, lmax3
    386   PARAMETER(kmax=8, lmax=8, lmax3=3)
    387   INTEGER kmaxm1, lmaxm1
    388   PARAMETER(kmaxm1=kmax-1, lmaxm1=lmax-1)
    389   INTEGER iimx7, jjmx7, jjmp1x7
    390   PARAMETER(iimx7=iim*kmaxm1, jjmx7=jjm*lmaxm1,  &
    391        jjmp1x7=jjmp1*lmaxm1)
    392   !
    393   ! output from ISCCP simulator
    394   REAL fq_isccp(klon,kmaxm1,lmaxm1,napisccp)
    395   REAL fq_is_true(klon,kmaxm1,lmaxm1,napisccp)
    396   REAL totalcldarea(klon,napisccp)
    397   REAL meanptop(klon,napisccp)
    398   REAL meantaucld(klon,napisccp)
    399   REAL boxtau(klon,ncolmx,napisccp)
    400   REAL boxptop(klon,ncolmx,napisccp)
    401   REAL zx_tmp_fi3d_bx(klon,ncolmx)
    402   REAL zx_tmp_3d_bx(iim,jjmp1,ncolmx)
    403   !
    404   REAL cld_fi3d(klon,lmax3)
    405   REAL cld_3d(iim,jjmp1,lmax3)
    406   !
    407   INTEGER iw, iwmax
    408   REAL wmin, pas_w
    409   !     PARAMETER(wmin=-100.,pas_w=10.,iwmax=30)
    410   !IM 051005     PARAMETER(wmin=-200.,pas_w=10.,iwmax=40)
    411   PARAMETER(wmin=-100.,pas_w=10.,iwmax=20)
    412   REAL o500(klon)
    413   !
    414 
    415   ! sorties ISCCP
    416 
    417   integer nid_isccp
    418   save nid_isccp       
    419   !$OMP THREADPRIVATE(nid_isccp)
    420 
    421   REAL zx_tau(kmaxm1), zx_pc(lmaxm1), zx_o500(iwmax)
    422   DATA zx_tau/0.0, 0.3, 1.3, 3.6, 9.4, 23., 60./
    423   SAVE zx_tau
    424   DATA zx_pc/180., 310., 440., 560., 680., 800., 1000./
    425   SAVE zx_pc
    426   !$OMP THREADPRIVATE(zx_tau,zx_pc)
    427   ! cldtopres pression au sommet des nuages
    428   REAL cldtopres(lmaxm1), cldtopres3(lmax3)
    429   DATA cldtopres/180., 310., 440., 560., 680., 800., 1000./
    430   DATA cldtopres3/440., 680., 1000./
    431   SAVE cldtopres,cldtopres3
    432   !$OMP THREADPRIVATE(cldtopres,cldtopres3)
    433   !IM 051005 BEG
    434   INTEGER komega, nhoriRD
    435 
    436   ! taulev: numero du niveau de tau dans les sorties ISCCP
    437   CHARACTER *4 taulev(kmaxm1)
    438   !     DATA taulev/'tau1','tau2','tau3','tau4','tau5','tau6','tau7'/
    439   DATA taulev/'tau0','tau1','tau2','tau3','tau4','tau5','tau6'/
    440   CHARACTER *3 pclev(lmaxm1)
    441   DATA pclev/'pc1','pc2','pc3','pc4','pc5','pc6','pc7'/
    442   SAVE taulev,pclev
    443   !$OMP THREADPRIVATE(taulev,pclev)
    444   !
    445   ! cnameisccp
    446   CHARACTER *29 cnameisccp(lmaxm1,kmaxm1)
    447   !IM bad 151205     DATA cnameisccp/'pc< 50hPa, tau< 0.3',
    448   DATA cnameisccp/'pc= 50-180hPa, tau< 0.3', &
    449        'pc= 180-310hPa, tau< 0.3', &
    450        'pc= 310-440hPa, tau< 0.3', &
    451        'pc= 440-560hPa, tau< 0.3', &
    452        'pc= 560-680hPa, tau< 0.3', &
    453        'pc= 680-800hPa, tau< 0.3', &
    454        'pc= 800-1000hPa, tau< 0.3', &
    455        'pc= 50-180hPa, tau= 0.3-1.3', &
    456        'pc= 180-310hPa, tau= 0.3-1.3', &
    457        'pc= 310-440hPa, tau= 0.3-1.3', &
    458        'pc= 440-560hPa, tau= 0.3-1.3', &
    459        'pc= 560-680hPa, tau= 0.3-1.3', &
    460        'pc= 680-800hPa, tau= 0.3-1.3', &
    461        'pc= 800-1000hPa, tau= 0.3-1.3', &
    462        'pc= 50-180hPa, tau= 1.3-3.6', &
    463        'pc= 180-310hPa, tau= 1.3-3.6', &
    464        'pc= 310-440hPa, tau= 1.3-3.6', &
    465        'pc= 440-560hPa, tau= 1.3-3.6', &
    466        'pc= 560-680hPa, tau= 1.3-3.6', &
    467        'pc= 680-800hPa, tau= 1.3-3.6', &
    468        'pc= 800-1000hPa, tau= 1.3-3.6', &
    469        'pc= 50-180hPa, tau= 3.6-9.4', &
    470        'pc= 180-310hPa, tau= 3.6-9.4', &
    471        'pc= 310-440hPa, tau= 3.6-9.4', &
    472        'pc= 440-560hPa, tau= 3.6-9.4', &
    473        'pc= 560-680hPa, tau= 3.6-9.4', &
    474        'pc= 680-800hPa, tau= 3.6-9.4', &
    475        'pc= 800-1000hPa, tau= 3.6-9.4', &
    476        'pc= 50-180hPa, tau= 9.4-23', &
    477        'pc= 180-310hPa, tau= 9.4-23', &
    478        'pc= 310-440hPa, tau= 9.4-23', &
    479        'pc= 440-560hPa, tau= 9.4-23', &
    480        'pc= 560-680hPa, tau= 9.4-23', &
    481        'pc= 680-800hPa, tau= 9.4-23', &
    482        'pc= 800-1000hPa, tau= 9.4-23', &
    483        'pc= 50-180hPa, tau= 23-60', &
    484        'pc= 180-310hPa, tau= 23-60', &
    485        'pc= 310-440hPa, tau= 23-60', &
    486        'pc= 440-560hPa, tau= 23-60', &
    487        'pc= 560-680hPa, tau= 23-60', &
    488        'pc= 680-800hPa, tau= 23-60', &
    489        'pc= 800-1000hPa, tau= 23-60', &
    490        'pc= 50-180hPa, tau> 60.', &
    491        'pc= 180-310hPa, tau> 60.', &
    492        'pc= 310-440hPa, tau> 60.', &
    493        'pc= 440-560hPa, tau> 60.', &
    494        'pc= 560-680hPa, tau> 60.', &
    495        'pc= 680-800hPa, tau> 60.', &
    496        'pc= 800-1000hPa, tau> 60.'/
    497   SAVE cnameisccp
    498   !$OMP THREADPRIVATE(cnameisccp)
    499   !
    500   !     REAL zx_lonx7(iimx7), zx_latx7(jjmp1x7)
    501   !     INTEGER nhorix7
    502293  !IM: region='3d' <==> sorties en global
    503294  CHARACTER*3 region
    504295  PARAMETER(region='3d')
    505   !
    506   !IM ISCCP simulator v3.4
    507   !
    508296  logical ok_hf
    509297  !
    510   integer nid_hf, nid_hf3d
    511   save ok_hf, nid_hf, nid_hf3d
    512   !$OMP THREADPRIVATE(ok_hf, nid_hf, nid_hf3d)
    513   !  QUESTION : noms de variables ?
     298  save ok_hf
     299  !$OMP THREADPRIVATE(ok_hf)
    514300
    515301  INTEGER        longcles
    … …  
    544330  !$OMP THREADPRIVATE(bas, top)
    545331
    546   REAL wdn(klon), tdn(klon), qdn(klon)
    547332  !
    548333  !=================================================================================================
    … …  
    559344  ! variables supplementaires de concvl
    560345  REAL Tconv(klon,klev)
    561   REAL ment(klon,klev,klev),sij(klon,klev,klev)
    562   REAL dd_t(klon,klev),dd_q(klon,klev)
     346  REAL sij(klon,klev,klev)
    563347
    564348  real, save :: alp_bl_prescr=0.
    … …  
    677461  REAL beta_prec_fisrt(klon,klev) ! taux de conv de l'eau cond (fisrt)
    678462  ! RomP <<<
    679   INTEGER       :: iii
    680   REAL          :: calday
    681463
    682464  !IM cf FH pour Tiedtke 080604
    … …  
    779561  REAL zdtime, zlongi
    780562  !
    781   CHARACTER*2 str2
    782   CHARACTER*2 iqn
    783   !
    784563  REAL qcheck
    785564  REAL z_avant(klon), z_apres(klon), z_factor(klon)
    … …  
    787566  !
    788567  REAL za, zb
    789   REAL zx_t, zx_qs, zdelta, zcor, zfra, zlvdcp, zlsdcp
     568  REAL zx_t, zx_qs, zdelta, zcor, zlvdcp, zlsdcp
    790569  real zqsat(klon,klev)
    791   INTEGER i, k, iq, ig, j, nsrf, ll, l, iiq, iff
     570  INTEGER i, k, iq, ig, j, nsrf, ll, l, iiq
    792571  REAL t_coup
    793572  PARAMETER (t_coup=234.0)
    … …  
    813592  ! -- convect43:
    814593  INTEGER ntra              ! nb traceurs pour convect4.3
    815   REAL pori_con(klon)    ! pressure at the origin level of lifted parcel
    816   REAL dtma_con(klon),dtlcl_con(klon)
    817594  REAL dtvpdt1(klon,klev), dtvpdq1(klon,klev)
    818595  REAL dplcldt(klon), dplcldr(klon)
    … …  
    849626  save ratqsbas,ratqshaut,tau_ratqs
    850627  !$OMP THREADPRIVATE(ratqsbas,ratqshaut,tau_ratqs)
    851   real zpt_conv(klon,klev)
    852628
    853629  ! Parametres lies au nouveau schema de nuages (SB, PDF)
    … …  
    872648  !======================================================================
    873649  !
    874   !IM cf. AM 081204 BEG
    875   !   declarations pour sortir sur une sous-region
    876   integer imin_ins,imax_ins,jmin_ins,jmax_ins
    877   save imin_ins,imax_ins,jmin_ins,jmax_ins
    878   !$OMP THREADPRIVATE(imin_ins,imax_ins,jmin_ins,jmax_ins)
    879   !      real lonmin_ins,lonmax_ins,latmin_ins
    880   !     s  ,latmax_ins
    881   !     data lonmin_ins,lonmax_ins,latmin_ins
    882   !    s  ,latmax_ins/
    883   ! valeurs initiales     s   -5.,20.,41.,55./   
    884   !    s   100.,130.,-20.,20./
    885   !    s   -180.,180.,-90.,90./
    886   !======================================================================
    887   !IM cf. AM 081204 END
    888 
     650 
    889651  !
    890652  integer itau_w   ! pas de temps ecriture = itap + itau_phy
    … …  
    900662  REAL tabcntr0( length       )
    901663  !
    902   INTEGER ndex2d(iim*jjmp1),ndex3d(iim*jjmp1*klev)
     664  INTEGER ndex2d(iim*jjmp1)
    903665  !IM
    904   INTEGER ndex2d1(iwmax)
    905666  !
    906667  !IM AMIP2 BEG
    … …  
    932693  REAL zx_tmp_fi2d(klon)      ! variable temporaire grille physique
    933694  REAL zx_tmp_fi3d(klon,klev) ! variable temporaire pour champs 3D
    934   REAL zx_tmp_fi3d1(klon,klev+1) !variable temporaire pour champs 3D (kelvp1)
    935   REAL(KIND=8) zx_tmp2_fi3d(klon,klev) ! variable temporaire pour champs 3D
    936   REAL zx_tmp_2d(iim,jjmp1), zx_tmp_3d(iim,jjmp1,klev)
     695  REAL zx_tmp_2d(iim,jjmp1)
    937696  REAL zx_lon(iim,jjmp1), zx_lat(iim,jjmp1)
    938697  !
    939   INTEGER nid_day, nid_mth, nid_ins, nid_mthnmc, nid_daynmc
    940   INTEGER nid_hfnmc, nid_day_seri, nid_ctesGCM
    941   SAVE nid_day, nid_mth, nid_ins, nid_mthnmc, nid_daynmc
    942   SAVE nid_hfnmc, nid_day_seri, nid_ctesGCM
    943   !$OMP THREADPRIVATE(nid_day, nid_mth, nid_ins)
    944   !$OMP THREADPRIVATE(nid_mthnmc, nid_daynmc, nid_hfnmc)
     698  INTEGER nid_day_seri, nid_ctesGCM
     699  SAVE nid_day_seri, nid_ctesGCM
    945700  !$OMP THREADPRIVATE(nid_day_seri,nid_ctesGCM)
    946701  !
    947702  !IM 280405 BEG
    948   INTEGER nid_bilKPins, nid_bilKPave
    949   SAVE nid_bilKPins, nid_bilKPave
    950   !$OMP THREADPRIVATE(nid_bilKPins, nid_bilKPave)
     703!  INTEGER nid_bilKPins, nid_bilKPave
     704!  SAVE nid_bilKPins, nid_bilKPave
     705!  !$OMP THREADPRIVATE(nid_bilKPins, nid_bilKPave)
    951706  !
    952707  REAL ve_lay(klon,klev) ! transport meri. de l'energie a chaque niveau vert.
    … …  
    955710  REAL uq_lay(klon,klev) ! transport zonal de l'eau a chaque niveau vert.
    956711  !
    957   INTEGER nhori, nvert, nvert1, nvert3
    958   REAL zsto, zsto1, zsto2
    959   REAL zstophy, zstorad, zstohf, zstoday, zstomth, zout
    960   REAL zcals(napisccp), zcalh(napisccp), zoutj(napisccp)
    961   REAL zout_isccp(napisccp)
    962   SAVE zcals, zcalh, zoutj, zout_isccp
    963   !$OMP THREADPRIVATE(zcals, zcalh, zoutj, zout_isccp)
    964 
     712  INTEGER nhori, nvert
     713  REAL zsto
     714  REAL zstophy, zout
     715 
    965716  real zjulian
    966717  save zjulian
    … …  
    979730  parameter (prof2d_on = 1, prof3d_on = 2, &
    980731       prof2d_av = 3, prof3d_av = 4)
    981   character*30 nom_fichier
    982   character*40 varname
    983   character*40 vartitle
    984   character*20 varunits
    985732  !     Variables liees au bilan d'energie et d'enthalpi
    986733  REAL ztsol(klon)
    987   REAL      h_vcol_tot, h_dair_tot, h_qw_tot, h_ql_tot &
    988        , h_qs_tot, qw_tot, ql_tot, qs_tot , ec_tot
    989   SAVE      h_vcol_tot, h_dair_tot, h_qw_tot, h_ql_tot &
    990        , h_qs_tot, qw_tot, ql_tot, qs_tot , ec_tot
    991   !$OMP THREADPRIVATE(h_vcol_tot, h_dair_tot, h_qw_tot, h_ql_tot)
    992   !$OMP THREADPRIVATE(h_qs_tot, qw_tot, ql_tot, qs_tot , ec_tot)
    993   REAL      d_h_vcol, d_h_dair, d_qt, d_qw, d_ql, d_qs, d_ec
     734  REAL      d_h_vcol, d_qt, d_qw, d_ql, d_qs, d_ec
    994735  REAL      d_h_vcol_phy
    995736  REAL      fs_bound, fq_bound
    … …  
    997738  !$OMP THREADPRIVATE(d_h_vcol_phy)
    998739  REAL      zero_v(klon)
    999   CHARACTER*15 ztit
     740  CHARACTER*40 ztit
    1000741  INTEGER   ip_ebil  ! PRINT level for energy conserv. diag.
    1001742  SAVE      ip_ebil
    … …  
    1005746  SAVE      if_ebil
    1006747  !$OMP THREADPRIVATE(if_ebil)
    1007   !+jld ec_conser
    1008   REAL ZRCPD
    1009   !-jld ec_conser
    1010748  REAL q2m(klon,nbsrf)  ! humidite a 2m
    1011749
    … …  
    1016754
    1017755
    1018   REAL re_top(klon), fl_top(klon) ! CDR at top of liquid water clouds
    1019 
    1020756  ! Aerosol optical properties
    1021757  CHARACTER*4, DIMENSION(naero_grp) :: rfname
    1022   REAL, DIMENSION(klon)          :: aerindex     ! POLDER aerosol index
    1023758  REAL, DIMENSION(klon,klev)     :: mass_solu_aero    ! total mass concentration for all soluble aerosols[ug/m3]
    1024759  REAL, DIMENSION(klon,klev)     :: mass_solu_aero_pi ! - " - (pre-industrial value)
    1025   INTEGER :: naero ! aerosol species
    1026760
    1027761  ! Parameters
    … …  
    1048782  LOGICAL,SAVE :: first=.true.
    1049783  !$OMP THREADPRIVATE(first)
    1050 
    1051   integer iunit
    1052784
    1053785  integer, save::  read_climoz ! read ozone climatology
    … …  
    12851017     ENDIF
    12861018     !
    1287      IF(ok_isccp.AND.iflag_con.LE.2) THEN
    1288         abort_message = 'ISCCP-like outputs may be available for KE' // &
    1289              '(iflag_con >= 3); for Tiedtke (iflag_con=-2) put ok_isccp=n'
    1290         CALL abort_gcm (modname,abort_message,1)
    1291      ENDIF
    12921019     !
    12931020     ! Initialiser les compteurs:
    … …  
    15251252     !$OMP BARRIER
    15261253
    1527 #undef histISCCP
    1528 #ifdef histISCCP
    1529      include "ini_histISCCP.h"
    1530 #endif
    15311254
    15321255     include "ini_histday_seri.h"
    … …  
    15571280        CALL VTe(VTphysiq)
    15581281        CALL VTb(VTinca)
    1559         !         iii = MOD(NINT(xjour),360)
    1560         !         calday = REAL(iii) + jH_cur
    15611282        calday = REAL(days_elapsed) + jH_cur
    15621283        WRITE(lunout,*) 'initial time chemini', days_elapsed, calday
    … …  
    29132634     !--fin STRAT AEROSOL
    29142635
    2915      !IM calcul nuages par le simulateur ISCCP
    2916      !
    2917 #ifdef histISCCP
    2918      IF (ok_isccp) THEN
    2919         !
    2920         !IM lecture invtau, tautab des fichiers formattes
    2921         !
    2922         IF (debut) THEN
    2923            !$OMP MASTER
    2924            !
    2925            open(99,file='tautab.formatted', FORM='FORMATTED')
    2926            read(99,'(f30.20)') tautab_omp
    2927            close(99)
    2928            !
    2929            open(99,file='invtau.formatted',form='FORMATTED')
    2930            read(99,'(i10)') invtau_omp
    2931 
    2932            !     print*,'calcul_simulISCCP invtau_omp',invtau_omp
    2933            !     write(6,'(a,8i10)') 'invtau_omp',(invtau_omp(i),i=1,100)
    2934 
    2935            close(99)
    2936            !$OMP END MASTER
    2937            !$OMP BARRIER
    2938            tautab=tautab_omp
    2939            invtau=invtau_omp
    2940            !
    2941         ENDIF !debut
    2942         !
    2943         !IM appel simulateur toutes les  NINT(freq_ISCCP/dtime) heures
    2944         IF (MOD(itap,NINT(freq_ISCCP/dtime)).EQ.0) THEN
    2945            include "calcul_simulISCCP.h"
    2946         ENDIF !(MOD(itap,NINT(freq_ISCCP/dtime))
    2947      ENDIF !ok_isccp
    2948 #endif
    29492636
    29502637     !   On prend la somme des fractions nuageuses et des contenus en eau
    … …  
    39673654
    39683655
    3969 #ifdef histISCCP
    3970   include "write_histISCCP.h"
    3971 #endif
    3972 
    39733656  include "write_histday_seri.h"
    39743657

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

    r1836 r1877  
    1 !$Id $
    2 !$Id $
     1!$Id$
    32MODULE phytrac_mod
    43!=================================================================================

  • LMDZ5/trunk/libf/phylmd/write_paramLMDZ_phy.h

    r1862 r1877  
    2727!
    2828      ndex2d = 0
    29       itau_w = itau_phy + itap + start_time * day_step / iphysiq
     29      itau_w = itau_phy + itap + int(start_time * day_step / iphysiq)
    3030!
    3131! Variables globales
Note: See TracChangeset for help on using the changeset viewer.