Changeset 3461 for LMDZ6/trunk/libf


Ignore:
Timestamp:
Feb 1, 2019, 8:29:43 PM (6 years ago)
Author:
oboucher
Message:

Writing out some parameters that were not output before
As usual cleaning out the routine but a lot more needed

File:
1 edited

Legend:

Unmodified
Added
Removed

  • LMDZ6/trunk/libf/phylmd/physiq_mod.F90

    r3460 r3461  
    335335    LOGICAL ok_gust ! pour activer l'effet des gust sur flux surface
    336336    PARAMETER (ok_gust=.FALSE.)
    337     integer iflag_radia     ! active ou non le rayonnement (MPL)
    338     save iflag_radia
     337    INTEGER iflag_radia     ! active ou non le rayonnement (MPL)
     338    SAVE iflag_radia
    339339    !$OMP THREADPRIVATE(iflag_radia)
    340340    !======================================================================
    … …  
    371371    !======================================================================
    372372    LOGICAL ok_journe ! sortir le fichier journalier
    373     save ok_journe
     373    SAVE ok_journe
    374374    !$OMP THREADPRIVATE(ok_journe)
    375375    !
    376376    LOGICAL ok_mensuel ! sortir le fichier mensuel
    377     save ok_mensuel
     377    SAVE ok_mensuel
    378378    !$OMP THREADPRIVATE(ok_mensuel)
    379379    !
    380380    LOGICAL ok_instan ! sortir le fichier instantane
    381     save ok_instan
     381    SAVE ok_instan
    382382    !$OMP THREADPRIVATE(ok_instan)
    383383    !
    384384    LOGICAL ok_LES ! sortir le fichier LES
    385     save ok_LES                           
     385    SAVE ok_LES                           
    386386    !$OMP THREADPRIVATE(ok_LES)                 
    387387    !
    388388    LOGICAL callstats ! sortir le fichier stats
    389     save callstats                           
     389    SAVE callstats                           
    390390    !$OMP THREADPRIVATE(callstats)                 
    391391    !
    … …  
    393393    PARAMETER (ok_region=.FALSE.)
    394394    !======================================================================
    395     real seuil_inversion
    396     save seuil_inversion
     395    REAL seuil_inversion
     396    SAVE seuil_inversion
    397397    !$OMP THREADPRIVATE(seuil_inversion)
    398     integer iflag_ratqs
    399     save iflag_ratqs
     398    INTEGER iflag_ratqs
     399    SAVE iflag_ratqs
    400400    !$OMP THREADPRIVATE(iflag_ratqs)
    401401    real facteur
    … …  
    404404    REAL tau_overturning_th(klon)
    405405
    406     integer lmax_th(klon)
    407     integer limbas(klon)
    408     real ratqscth(klon,klev)
    409     real ratqsdiff(klon,klev)
    410     real zqsatth(klon,klev)
     406    INTEGER lmax_th(klon)
     407    INTEGER limbas(klon)
     408    REAL ratqscth(klon,klev)
     409    REAL ratqsdiff(klon,klev)
     410    REAL zqsatth(klon,klev)
    411411
    412412    !======================================================================
    … …  
    505505    CHARACTER*3 region
    506506    PARAMETER(region='3d')
    507     logical ok_hf
    508     !
    509     save ok_hf
     507    LOGICAL ok_hf
     508    !
     509    SAVE ok_hf
    510510    !$OMP THREADPRIVATE(ok_hf)
    511511
    512     INTEGER,PARAMETER :: longcles=20
    513     REAL,SAVE :: clesphy0(longcles)
     512    INTEGER, PARAMETER :: longcles=20
     513    REAL, SAVE :: clesphy0(longcles)
    514514    !$OMP THREADPRIVATE(clesphy0)
    515515    !
    516516    ! Variables propres a la physique
    517     INTEGER itap
    518     SAVE itap                   ! compteur pour la physique
     517    INTEGER, SAVE :: itap         ! compteur pour la physique
    519518    !$OMP THREADPRIVATE(itap)
    520519
    … …  
    522521    !$OMP THREADPRIVATE(abortphy)
    523522    !
    524     REAL,save ::  solarlong0
     523    REAL,SAVE ::  solarlong0
    525524    !$OMP THREADPRIVATE(solarlong0)
    526525
    … …  
    539538    ! Variables liees a la convection de K. Emanuel (sb):
    540539    !
    541     REAL bas, top             ! cloud base and top levels
    542     SAVE bas
    543     SAVE top
     540    REAL, SAVE :: bas, top             ! cloud base and top levels
    544541    !$OMP THREADPRIVATE(bas, top)
    545542    !------------------------------------------------------------------
    … …  
    576573!!    REAL, DIMENSION(klon,klev)     :: dql_sat
    577574
    578     real, save :: alp_bl_prescr=0.
    579     real, save :: ale_bl_prescr=0.
    580 
    581     real, save :: wake_s_min_lsp=0.1
    582 
     575    REAL, SAVE :: alp_bl_prescr=0.
     576    REAL, SAVE :: ale_bl_prescr=0.
     577    REAL, SAVE :: wake_s_min_lsp=0.1
    583578    !$OMP THREADPRIVATE(alp_bl_prescr,ale_bl_prescr)
    584579    !$OMP THREADPRIVATE(wake_s_min_lsp)
    585580
    586 
    587     real ok_wk_lsp(klon)
     581    REAL ok_wk_lsp(klon)
    588582
    589583    !RC
    … …  
    735729    REAL :: jD_eq
    736730
    737     LOGICAL, parameter :: new_orbit = .true.
     731    LOGICAL, parameter :: new_orbit = .TRUE.
    738732
    739733    !
    … …  
    921915    INTEGER kcbot(klon), kctop(klon), kdtop(klon)
    922916    !
    923     real ratqsbas,ratqshaut,tau_ratqs
    924     save ratqsbas,ratqshaut,tau_ratqs
     917    REAL ratqsbas,ratqshaut,tau_ratqs
     918    SAVE ratqsbas,ratqshaut,tau_ratqs
    925919    !$OMP THREADPRIVATE(ratqsbas,ratqshaut,tau_ratqs)
    926920    REAL, SAVE :: ratqsp0=50000., ratqsdp=20000.
    … …  
    928922
    929923    ! Parametres lies au nouveau schema de nuages (SB, PDF)
    930     real fact_cldcon
    931     real facttemps
    932     logical ok_newmicro
    933     save ok_newmicro
     924    REAL, SAVE :: fact_cldcon
     925    REAL, SAVE :: facttemps
     926    !$OMP THREADPRIVATE(fact_cldcon,facttemps)
     927    LOGICAL, SAVE :: ok_newmicro
    934928    !$OMP THREADPRIVATE(ok_newmicro)
    935     !real ref_liq_pi(klon,klev), ref_ice_pi(klon,klev)
    936     save fact_cldcon,facttemps
    937     !$OMP THREADPRIVATE(fact_cldcon,facttemps)
    938 
    939     integer iflag_cld_th
    940     save iflag_cld_th
     929
     930    INTEGER, SAVE :: iflag_cld_th
    941931    !$OMP THREADPRIVATE(iflag_cld_th)
    942932!IM logical ptconv(klon,klev)  !passe dans phys_local_var_mod
    943933    !IM cf. AM 081204 BEG
    944     logical ptconvth(klon,klev)
     934    LOGICAL ptconvth(klon,klev)
    945935    !IM cf. AM 081204 END
    946936    !
    … …  
    949939    !======================================================================
    950940    !
    951 
    952941    !
    953942!JLD    integer itau_w   ! pas de temps ecriture = itap + itau_phy
    … …  
    10151004!JLD    REAL zstophy, zout
    10161005
    1017     character*20 modname
    1018     character*80 abort_message
    1019     logical, save ::  ok_sync, ok_sync_omp
     1006    CHARACTER*20 modname
     1007    CHARACTER*80 abort_message
     1008    LOGICAL, SAVE ::  ok_sync, ok_sync_omp
    10201009    !$OMP THREADPRIVATE(ok_sync)
    1021     real date0
     1010    REAL date0
    10221011
    10231012    ! essai writephys
    1024     integer fid_day, fid_mth, fid_ins
    1025     parameter (fid_ins = 1, fid_day = 2, fid_mth = 3)
    1026     integer prof2d_on, prof3d_on, prof2d_av, prof3d_av
    1027     parameter (prof2d_on = 1, prof3d_on = 2, &
    1028          prof2d_av = 3, prof3d_av = 4)
     1013    INTEGER fid_day, fid_mth, fid_ins
     1014    PARAMETER (fid_ins = 1, fid_day = 2, fid_mth = 3)
     1015    INTEGER prof2d_on, prof3d_on, prof2d_av, prof3d_av
     1016    PARAMETER (prof2d_on = 1, prof3d_on = 2, prof2d_av = 3, prof3d_av = 4)
    10291017    REAL ztsol(klon)
    10301018    REAL q2m(klon,nbsrf)  ! humidite a 2m
    … …  
    10781066    ! Declaration des constantes et des fonctions thermodynamiques
    10791067    !
    1080     LOGICAL,SAVE :: first=.true.
     1068    LOGICAL,SAVE :: first=.TRUE.
    10811069    !$OMP THREADPRIVATE(first)
    10821070
    … …  
    11721160    REAL zzz
    11731161    !albedo SB >>>
    1174     real,dimension(6),save :: SFRWL
     1162    REAL,DIMENSION(6), SAVE :: SFRWL
    11751163!$OMP THREADPRIVATE(SFRWL)
    11761164    !albedo SB <<<
    11771165
    11781166    !--OB variables for mass fixer (hard coded for now)
    1179     logical, parameter :: mass_fixer=.false.
    1180     real qql1(klon),qql2(klon),corrqql
     1167    LOGICAL, PARAMETER :: mass_fixer=.FALSE.
     1168    REAL qql1(klon),qql2(klon),corrqql
    11811169
    11821170    REAL pi
    … …  
    12241212
    12251213    ! Quick check on pressure levels:
    1226     call assert(paprs(:, nbp_lev + 1) < paprs(:, nbp_lev), &
     1214    CALL assert(paprs(:, nbp_lev + 1) < paprs(:, nbp_lev), &
    12271215            "physiq_mod paprs bad order")
    12281216
    … …  
    12321220       !CR:nvelles variables convection/poches froides
    12331221
    1234        print*, '================================================='
    1235        print*, 'Allocation des variables locales et sauvegardees'
     1222       WRITE(lunout,*) '================================================='
     1223       WRITE(lunout,*) 'Allocation des variables locales et sauvegardees'
     1224       WRITE(lunout,*) '================================================='
    12361225       CALL phys_local_var_init
    12371226       !
    … …  
    12741263
    12751264       itau_con=0
    1276        first=.false.
     1265       first=.FALSE.
    12771266
    12781267    ENDIF  ! first
    … …  
    13031292! secondes
    13041293       tau_gl=86400.*tau_gl
    1305        print*,'debut physiq_mod tau_gl=',tau_gl
     1294       WRITE(lunout,*) 'debut physiq_mod tau_gl=',tau_gl
     1295
    13061296       CALL getin_p('iflag_alp_wk_cond', iflag_alp_wk_cond)
    13071297       CALL getin_p('random_notrig_max',random_notrig_max)
    … …  
    13341324       CALL getin_p('NVM',nvm_lmdz)
    13351325
     1326       WRITE(lunout,*) 'iflag_alp_wk_cond=',  iflag_alp_wk_cond
     1327       WRITE(lunout,*) 'random_ntrig_max=',   random_notrig_max
     1328       WRITE(lunout,*) 'ok_adjwk=',           ok_adjwk
     1329       WRITE(lunout,*) 'iflag_adjwk=',        iflag_adjwk
     1330       WRITE(lunout,*) 'qtcon_multistep_max=',dtcon_multistep_max
     1331       WRITE(lunout,*) 'qdcon_multistep_max=',dqcon_multistep_max
     1332       WRITE(lunout,*) 'ratqsp0=',            ratqsp0
     1333       WRITE(lunout,*) 'ratqsdp=',            ratqsdp
     1334       WRITE(lunout,*) 'iflag_wake_tend=',    iflag_wake_tend
     1335       WRITE(lunout,*) 'ok_bad_ecmwf_thermo=',ok_bad_ecmwf_thermo
     1336       WRITE(lunout,*) 'ok_bug_cv_trac=',     ok_bug_cv_trac
     1337       WRITE(lunout,*) 'ok_bug_split_th=',    ok_bug_split_th
     1338       WRITE(lunout,*) 'fl_ebil=',            fl_ebil
     1339       WRITE(lunout,*) 'fl_cor_ebil=',        fl_cor_ebil
     1340       WRITE(lunout,*) 'iflag_phytrac=',      iflag_phytrac
     1341       WRITE(lunout,*) 'NVM=',                nvm_lmdz
     1342
    13361343       !--PC: defining fields to be exchanged between LMDz, ORCHIDEE and NEMO
    13371344       WRITE(lunout,*) 'Call to infocfields from physiq'
    … …  
    13971404       ENDIF
    13981405
    1399 
    1400 
    14011406       rnebcon0(:,:) = 0.0
    14021407       clwcon0(:,:) = 0.0
    … …  
    14681473          abort_message='nbre de pas de temps physique n est pas multiple ' &
    14691474               // 'de nbapp_cv'
    1470           call abort_physic(modname,abort_message,1)
     1475          CALL abort_physic(modname,abort_message,1)
    14711476       ENDIF
    14721477       IF (MOD(NINT(86400./phys_tstep),nbapp_wk).EQ.0) THEN
    … …  
    14801485          abort_message='nbre de pas de temps physique n est pas multiple ' &
    14811486               // 'de nbapp_wk'
    1482           call abort_physic(modname,abort_message,1)
     1487          CALL abort_physic(modname,abort_message,1)
    14831488       ENDIF
    14841489       ! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    … …  
    15361541       ! FH : if ok_sync=.true. , the time axis is written at each time step
    15371542       ! in the output files. Only at the end in the opposite case
    1538        ok_sync_omp=.false.
     1543       ok_sync_omp=.FALSE.
    15391544       CALL getin('ok_sync',ok_sync_omp)
    15401545       CALL phys_output_open(longitude_deg,latitude_deg,nCFMIP,tabijGCM, &
    … …  
    18841889
    18851890       !albedo SB >>>
    1886        select case(nsw)
    1887        case(2)
     1891       SELECT CASE(nsw)
     1892       CASE(2)
    18881893          SFRWL(1)=0.45538747
    18891894          SFRWL(2)=0.54461211
    1890        case(4)
     1895       CASE(4)
    18911896          SFRWL(1)=0.45538747
    18921897          SFRWL(2)=0.32870591
    18931898          SFRWL(3)=0.18568763
    18941899          SFRWL(4)=3.02191470E-02
    1895        case(6)
     1900       CASE(6)
    18961901          SFRWL(1)=1.28432794E-03
    18971902          SFRWL(2)=0.12304168
    … …  
    19001905          SFRWL(5)=0.18568763
    19011906          SFRWL(6)=3.02191470E-02
    1902        end select
     1907       END SELECT
    19031908
    19041909
    … …  
    19731978    ENDIF
    19741979
    1975 
    19761980    ! Tendances bidons pour les processus qui n'affectent pas certaines
    19771981    ! variables.
    … …  
    22072211                      ro3i, 'C', press_cen_climoz, pplay, wo, paprs(:,1),    &
    22082212                      time_climoz )
    2209           END IF
     2213          ENDIF
    22102214          ! Convert from mole fraction of ozone to column density of ozone in a
    22112215          ! cell, in kDU:
    … …  
    22272231            (du0,dv0,d_t_eva,d_q_eva,d_ql_eva,d_qi_eva,paprs,&
    22282232               'eva',abortphy,flag_inhib_tend,itap,0)
    2229     call prt_enerbil('eva',itap)
     2233    CALL prt_enerbil('eva',itap)
    22302234
    22312235    !=========================================================================
    … …  
    23842388       ENDIF
    23852389
    2386 
    2387 
    23882390       CALL pbl_surface(  &
    23892391            phys_tstep,     date0,     itap,    days_elapsed+1, &
    … …  
    24552457
    24562458
    2457 
    2458 
    2459 
    24602459       !---------------------------------------------------------------------
    24612460       ! ajout des tendances de la diffusion turbulente
    … …  
    24692468               'vdf',abortphy,flag_inhib_tend,itap,0)
    24702469       ENDIF
    2471        call prt_enerbil('vdf',itap)
     2470       CALL prt_enerbil('vdf',itap)
    24722471       !--------------------------------------------------------------------
    24732472
    … …  
    29372936    CALL add_phys_tend(d_u_con, d_v_con, d_t_con, d_q_con, dql0, dqi0, paprs, &
    29382937         'convection',abortphy,flag_inhib_tend,itap,0)
    2939     call prt_enerbil('convection',itap)
     2938    CALL prt_enerbil('convection',itap)
    29402939
    29412940    !-------------------------------------------------------------------------
    … …  
    30683067       CALL add_phys_tend(du0,dv0,d_t_wake,d_q_wake,dql0,dqi0,paprs,'wake', &
    30693068            abortphy,flag_inhib_tend,itap,0)
    3070        call prt_enerbil('wake',itap)
     3069       CALL prt_enerbil('wake',itap)
    30713070       !------------------------------------------------------------------------
    30723071
    … …  
    30773076            (d_deltat_wk, d_deltaq_wk, d_s_wk, d_dens_a_wk, d_dens_wk, wake_k, &
    30783077             'wake', abortphy)
    3079           call prt_enerbil('wake',itap)
     3078          CALL prt_enerbil('wake',itap)
    30803079       ENDIF   ! (iflag_wake_tend .GT. 0.)
    30813080       !
    … …  
    32193218                   (d_deltat_the, d_deltaq_the, dsig0, ddens0, ddens0, wake_k, 'the', abortphy)
    32203219             ENDIF
    3221              call prt_enerbil('the',itap)
     3220             CALL prt_enerbil('the',itap)
    32223221          !
    32233222          ENDIF  ! (mod(iflag_pbl_split/10,10) .GE. 1)
    … …  
    32253224          CALL add_phys_tend(d_u_ajs,d_v_ajs,d_t_ajs,d_q_ajs,  &
    32263225                             dql0,dqi0,paprs,'thermals', abortphy,flag_inhib_tend,itap,0)
    3227           call prt_enerbil('thermals',itap)
     3226          CALL prt_enerbil('thermals',itap)
    32283227          !
    32293228!
    … …  
    32883287          CALL add_phys_tend(du0,dv0,d_t_ajsb,d_q_ajsb,dql0,dqi0,paprs, &
    32893288               'ajsb',abortphy,flag_inhib_tend,itap,0)
    3290           call prt_enerbil('ajsb',itap)
     3289          CALL prt_enerbil('ajsb',itap)
    32913290          d_t_ajs(:,:)=d_t_ajs(:,:)+d_t_ajsb(:,:)
    32923291          d_q_ajs(:,:)=d_q_ajs(:,:)+d_q_ajsb(:,:)
    … …  
    33393338    CALL add_phys_tend(du0,dv0,d_t_lsc,d_q_lsc,d_ql_lsc,d_qi_lsc,paprs, &
    33403339         'lsc',abortphy,flag_inhib_tend,itap,0)
    3341     call prt_enerbil('lsc',itap)
     3340    CALL prt_enerbil('lsc',itap)
    33423341    rain_num(:)=0.
    33433342    DO k = 1, klev
    … …  
    41554154
    41564155    CALL add_phys_tend(du0,dv0,d_t_swr,dq0,dql0,dqi0,paprs,'SW',abortphy,flag_inhib_tend,itap,0)
    4157     call prt_enerbil('SW',itap)
     4156    CALL prt_enerbil('SW',itap)
    41584157    CALL add_phys_tend(du0,dv0,d_t_lwr,dq0,dql0,dqi0,paprs,'LW',abortphy,flag_inhib_tend,itap,0)
    4159     call prt_enerbil('LW',itap)
     4158    CALL prt_enerbil('LW',itap)
    41604159
    41614160    !
    … …  
    42284227       CALL add_phys_tend(d_u_oro,d_v_oro,d_t_oro,dq0,dql0,dqi0,paprs,'oro', &
    42294228            abortphy,flag_inhib_tend,itap,0)
    4230        call prt_enerbil('oro',itap)
     4229       CALL prt_enerbil('oro',itap)
    42314230       !----------------------------------------------------------------------
    42324231       !
    … …  
    42754274       CALL add_phys_tend(d_u_lif, d_v_lif, d_t_lif, dq0, dql0, dqi0, paprs, &
    42764275            'lif', abortphy,flag_inhib_tend,itap,0)
    4277        call prt_enerbil('lif',itap)
     4276       CALL prt_enerbil('lif',itap)
    42784277    ENDIF ! fin de test sur ok_orolf
    42794278
    … …  
    42994298       CALL add_phys_tend(du_gwd_hines, dv_gwd_hines, d_t_hin, dq0, dql0, &
    43004299            dqi0, paprs, 'hin', abortphy,flag_inhib_tend,itap,0)
    4301        call prt_enerbil('hin',itap)
     4300       CALL prt_enerbil('hin',itap)
    43024301    ENDIF
    43034302
    … …  
    43184317       CALL add_phys_tend(du_gwd_front, dv_gwd_front, dt0, dq0, dql0, dqi0, &
    43194318            paprs, 'front_gwd_rando', abortphy,flag_inhib_tend,itap,0)
    4320        call prt_enerbil('front_gwd_rando',itap)
     4319       CALL prt_enerbil('front_gwd_rando',itap)
    43214320    ENDIF
    43224321
    … …  
    43274326       CALL add_phys_tend(du_gwd_rando, dv_gwd_rando, dt0, dq0, dql0, dqi0, &
    43284327            paprs, 'flott_gwd_rando', abortphy,flag_inhib_tend,itap,0)
    4329        call prt_enerbil('flott_gwd_rando',itap)
     4328       CALL prt_enerbil('flott_gwd_rando',itap)
    43304329       zustr_gwd_rando=0.
    43314330       zvstr_gwd_rando=0.
    … …  
    43734372    !IM cf. FLott END
    43744373    !DC Calcul de la tendance due au methane
    4375     IF(ok_qch4) THEN
     4374    IF (ok_qch4) THEN
    43764375       CALL METHOX(1,klon,klon,klev,q_seri,d_q_ch4,pplay)
    43774376       ! ajout de la tendance d'humidite due au methane
    … …  
    43904389! Inititialization
    43914390!------------------
    4392 
    4393    
    43944391
    43954392       addtkeoro=0   
    … …  
    44034400       alphatkeoro=min(max(0.,alphatkeoro),1.)
    44044401
    4405        smallscales_tkeoro=.false.   
     4402       smallscales_tkeoro=.FALSE.   
    44064403       CALL getin_p('smallscales_tkeoro',smallscales_tkeoro)
    44074404
    44084405
    4409         dtadd(:,:)=0.
    4410         duadd(:,:)=0.
    4411         dvadd(:,:)=0.
    4412 
    4413 
     4406       dtadd(:,:)=0.
     4407       duadd(:,:)=0.
     4408       dvadd(:,:)=0.
    44144409
    44154410! Choices for addtkeoro:
    … …  
    44264421
    44274422
    4428 
    44294423  IF (addtkeoro .EQ. 1 ) THEN
    44304424
    … …  
    44344428  ELSE IF (addtkeoro .EQ. 2) THEN
    44354429
    4436 
    4437 
    4438        IF (smallscales_tkeoro) THEN
     4430     IF (smallscales_tkeoro) THEN
    44394431       igwd=0
    44404432       DO i=1,klon
    … …  
    44594451             igwd=igwd+1
    44604452             idx(igwd)=i
    4461           ENDIF
    4462        ENDDO
    4463 
    4464        END IF
    4465 
    4466 
    4467 
    4468 
    4469        CALL drag_noro_strato(addtkeoro,klon,klev,phys_tstep,paprs,pplay, &
     4453        ENDIF
     4454       ENDDO
     4455
     4456     ENDIF
     4457
     4458     CALL drag_noro_strato(addtkeoro,klon,klev,phys_tstep,paprs,pplay, &
    44704459               zmea,zstd, zsig, zgam, zthe,zpic,zval, &
    44714460               igwd,idx,itest, &
    … …  
    44744463               d_t_oro_gw, d_u_oro_gw, d_v_oro_gw)
    44754464
    4476             zustrdr(:)=0.
    4477             zvstrdr(:)=0.
    4478             zulow(:)=0.
    4479             zvlow(:)=0.
    4480 
    4481             duadd(:,:)=alphatkeoro*d_u_oro_gw(:,:)
    4482             dvadd(:,:)=alphatkeoro*d_v_oro_gw(:,:)
    4483  END IF
    4484    
     4465     zustrdr(:)=0.
     4466     zvstrdr(:)=0.
     4467     zulow(:)=0.
     4468     zvlow(:)=0.
     4469
     4470     duadd(:,:)=alphatkeoro*d_u_oro_gw(:,:)
     4471     dvadd(:,:)=alphatkeoro*d_v_oro_gw(:,:)
     4472  ENDIF
    44854473
    44864474
    … …  
    44934481
    44944482
    4495 
    44964483       ENDIF
    44974484!      -----
    44984485!===============================================================
    4499 
    45004486
    45014487
    … …  
    49674953
    49684954! On remet des variables a .false. apres un premier appel
    4969     if (debut) then
     4955    IF (debut) THEN
    49704956#ifdef CPP_XIOS
    49714957      swaero_diag=.FALSE.
    … …  
    49754961!      write (lunout,*)'ok_4xCO2atm= ',swaero_diag, swaerofree_diag, dryaod_diag, ok_4xCO2atm
    49764962
    4977       IF (is_master) then
     4963      IF (is_master) THEN
    49784964        !--setting up swaero_diag to TRUE in XIOS case
    49794965        IF (xios_field_is_active("topswad").OR.xios_field_is_active("topswad0").OR. &
    … …  
    50064992           xios_field_is_active("rld4co2").OR.xios_field_is_active("rldcs4co2")) &
    50074993           ok_4xCO2atm=.TRUE.
    5008       endif
     4994      ENDIF
    50094995      !$OMP BARRIER
    5010       call bcast(swaero_diag)
    5011       call bcast(swaerofree_diag)
    5012       call bcast(dryaod_diag)
    5013       call bcast(ok_4xCO2atm)
     4996      CALL bcast(swaero_diag)
     4997      CALL bcast(swaerofree_diag)
     4998      CALL bcast(dryaod_diag)
     4999      CALL bcast(ok_4xCO2atm)
    50145000!      write (lunout,*)'ok_4xCO2atm= ',swaero_diag, swaerofree_diag, dryaod_diag, ok_4xCO2atm
    50155001#endif
    5016     endif
     5002    ENDIF
    50175003
    50185004    !====================================================================
    … …  
    50525038       IF (is_omp_master) CALL xios_context_finalize
    50535039#endif
    5054        print *,' physiq fin, nombre de steps ou cvpas = 1 : ', Ncvpaseq1
     5040       WRITE(lunout,*) ' physiq fin, nombre de steps ou cvpas = 1 : ', Ncvpaseq1
    50555041    ENDIF
    50565042
    50575043    !      first=.false.
    50585044
    5059 
    50605045  END SUBROUTINE physiq
    50615046
Note: See TracChangeset for help on using the changeset viewer.