Changeset 6445 for CONFIG


Ignore:
Timestamp:
05/11/23 14:00:20 (13 months ago)
Author:
aclsce
Message:

Modified to be able to work with

  • LMDZ testing rev 4515
  • ORCHIDEE rev 7991
  • ICOSA_LMDZ rev 4501
  • DYNAMICO rev 3558d704

Removed SOURCES not needed anymore due to the use above revisions.

Location:
CONFIG/UNIFORM/v7/IPSLCM7
Files:
3 deleted
15 edited

Legend:

Unmodified
Added
Removed

  • CONFIG/UNIFORM/v7/IPSLCM7/ARCH/arch-X64_IRENE-AMD.path

    r5479 r6445  
    22NETCDF_LIBDIR="-L $NETCDFC_LIBDIR -L $NETCDFFORTRAN_LIBDIR" 
    33NETCDF_LIB="-lnetcdf -lnetcdff" 
     4 
     5NETCDF95_INCDIR=-I/ccc/work/cont003/gencmip6/guezl/build/Libraries_ifort_17_release/NetCDF95/modules 
     6NETCDF95_LIBDIR=-L/ccc/work/cont003/gencmip6/guezl/build/Libraries_ifort_17_release/NetCDF95 
     7NETCDF95_LIB=-lnetcdf95  
    48 
    59HDF5_INCDIR="" 

  • CONFIG/UNIFORM/v7/IPSLCM7/ARCH/arch-X64_IRENE.path

    r5479 r6445  
    22NETCDF_LIBDIR="-L $NETCDFC_LIBDIR -L $NETCDFFORTRAN_LIBDIR" 
    33NETCDF_LIB="-lnetcdf -lnetcdff" 
     4 
     5NETCDF95_INCDIR=-I/ccc/work/cont003/gencmip6/guezl/build/Libraries_ifort_17_release/NetCDF95/modules 
     6NETCDF95_LIBDIR=-L/ccc/work/cont003/gencmip6/guezl/build/Libraries_ifort_17_release/NetCDF95 
     7NETCDF95_LIB=-lnetcdf95  
    48 
    59HDF5_INCDIR="" 

  • CONFIG/UNIFORM/v7/IPSLCM7/EXPERIMENTS/ICOLMDZ/CREATE_amip/COMP/dynamico.card

    r6362 r6445  
    2929        (${SUBMIT_DIR}/PARAM/run_dynamico.def_nbp${RESOL_NBP}, run_dynamico.def),\ 
    3030        (${SUBMIT_DIR}/PARAM/iodef.xml, .),\ 
    31         (${SUBMIT_DIR}/PARAM/sponge_dynamico.xml, .),\ 
     31        (${MODIPSL}/modeles/DYNAMICO/xml/sponge_dynamico.xml, .),\ 
    3232        (${MODIPSL}/modeles/DYNAMICO/xml/context_dynamico.xml, .),\ 
    3333        (${MODIPSL}/modeles/DYNAMICO/xml/context_input_dynamico.xml, .),\ 

  • CONFIG/UNIFORM/v7/IPSLCM7/EXPERIMENTS/ICOLMDZ/CREATE_hadisst/COMP/dynamico.card

    r6362 r6445  
    2929        (${SUBMIT_DIR}/PARAM/run_dynamico.def_nbp${RESOL_NBP}, run_dynamico.def),\ 
    3030        (${SUBMIT_DIR}/PARAM/iodef.xml, .),\ 
    31         (${SUBMIT_DIR}/PARAM/sponge_dynamico.xml, .),\ 
     31        (${MODIPSL}/modeles/DYNAMICO/xml/sponge_dynamico.xml, .),\ 
    3232        (${MODIPSL}/modeles/DYNAMICO/xml/context_dynamico.xml, .),\ 
    3333        (${MODIPSL}/modeles/DYNAMICO/xml/context_input_dynamico.xml, .),\ 

  • CONFIG/UNIFORM/v7/IPSLCM7/EXPERIMENTS/ICOLMDZ/aquaplanet/CREATE_aquaplanet/COMP/dynamico.card

    r6362 r6445  
    2929        (${SUBMIT_DIR}/PARAM/run_dynamico.def_nbp${RESOL_NBP}, run_dynamico.def),\ 
    3030        (${SUBMIT_DIR}/PARAM/iodef.xml, .),\ 
    31         (${SUBMIT_DIR}/PARAM/sponge_dynamico.xml, .),\ 
     31        (${MODIPSL}/modeles/DYNAMICO/xml/sponge_dynamico.xml, .),\ 
    3232        (${MODIPSL}/modeles/DYNAMICO/xml/context_dynamico.xml, .),\ 
    3333        (${MODIPSL}/modeles/DYNAMICO/xml/context_input_dynamico.xml, .),\ 

  • CONFIG/UNIFORM/v7/IPSLCM7/EXPERIMENTS/ICOLMDZ/aquaplanet/climAQUA_360d/COMP/dynamico.card

    r6362 r6445  
    2929        (${SUBMIT_DIR}/PARAM/run_dynamico.def_nbp${RESOL_NBP}, run_dynamico.def),\ 
    3030        (${SUBMIT_DIR}/PARAM/iodef.xml, .),\ 
    31         (${SUBMIT_DIR}/PARAM/sponge_dynamico.xml, .),\ 
     31        (${MODIPSL}/modeles/DYNAMICO/xml/sponge_dynamico.xml, .),\ 
    3232        (${MODIPSL}/modeles/DYNAMICO/xml/context_dynamico.xml, .),\ 
    3333        (${MODIPSL}/modeles/DYNAMICO/xml/context_input_dynamico.xml, .),\ 

  • CONFIG/UNIFORM/v7/IPSLCM7/EXPERIMENTS/ICOLMDZ/clim_noleap/COMP/dynamico.card

    r6362 r6445  
    3838        (${SUBMIT_DIR}/PARAM/run_dynamico.def_nbp${RESOL_NBP}, run_dynamico.def),\ 
    3939        (${SUBMIT_DIR}/PARAM/iodef.xml, .),\ 
    40         (${SUBMIT_DIR}/PARAM/sponge_dynamico.xml, .),\ 
     40        (${MODIPSL}/modeles/DYNAMICO/xml/sponge_dynamico.xml, .),\ 
    4141        (${MODIPSL}/modeles/DYNAMICO/xml/context_dynamico.xml, .),\ 
    4242        (${MODIPSL}/modeles/DYNAMICO/xml/context_input_dynamico.xml, .),\ 

  • CONFIG/UNIFORM/v7/IPSLCM7/EXPERIMENTS/ICOLMDZOR/clim_noleap/COMP/dynamico.card

    r6417 r6445  
    3838        (${SUBMIT_DIR}/PARAM/run_dynamico.def_nbp${RESOL_NBP}, run_dynamico.def),\ 
    3939        (${SUBMIT_DIR}/PARAM/iodef.xml, .),\ 
    40         (${SUBMIT_DIR}/PARAM/sponge_dynamico.xml, .),\ 
     40        (${MODIPSL}/modeles/DYNAMICO/xml/sponge_dynamico.xml, .),\ 
    4141        (${MODIPSL}/modeles/DYNAMICO/xml/context_dynamico.xml, .),\ 
    4242        (${MODIPSL}/modeles/DYNAMICO/xml/context_input_dynamico.xml, .),\ 

  • CONFIG/UNIFORM/v7/IPSLCM7/EXPERIMENTS/ICOLMDZOR/clim_noleap_ORCHv4/COMP/dynamico.card

    r6407 r6445  
    3838        (${SUBMIT_DIR}/PARAM/run_dynamico.def_nbp${RESOL_NBP}, run_dynamico.def),\ 
    3939        (${SUBMIT_DIR}/PARAM/iodef.xml, .),\ 
    40         (${SUBMIT_DIR}/PARAM/sponge_dynamico.xml, .),\ 
     40        (${MODIPSL}/modeles/DYNAMICO/xml/sponge_dynamico.xml, .),\ 
    4141        (${MODIPSL}/modeles/DYNAMICO/xml/context_dynamico.xml, .),\ 
    4242        (${MODIPSL}/modeles/DYNAMICO/xml/context_input_dynamico.xml, .),\ 

  • CONFIG/UNIFORM/v7/IPSLCM7/EXPERIMENTS/ICOLMDZOR/highresmip/COMP/dynamico.card

    r6407 r6445  
    3838        (${SUBMIT_DIR}/PARAM/run_dynamico.def_nbp${RESOL_NBP}, run_dynamico.def),\ 
    3939        (${SUBMIT_DIR}/PARAM/iodef.xml, .),\ 
    40         (${SUBMIT_DIR}/PARAM/sponge_dynamico.xml, .),\ 
     40        (${MODIPSL}/modeles/DYNAMICO/xml/sponge_dynamico.xml, .),\ 
    4141        (${MODIPSL}/modeles/DYNAMICO/xml/context_dynamico.xml, .),\ 
    4242        (${MODIPSL}/modeles/DYNAMICO/xml/context_input_dynamico.xml, .),\ 

  • CONFIG/UNIFORM/v7/IPSLCM7/EXPERIMENTS/ICOLMDZOR/highresmip_CMIP6/COMP/dynamico.card

    r6407 r6445  
    3838        (${SUBMIT_DIR}/PARAM/run_dynamico.def_nbp${RESOL_NBP}, run_dynamico.def),\ 
    3939        (${SUBMIT_DIR}/PARAM/iodef.xml, .),\ 
    40         (${SUBMIT_DIR}/PARAM/sponge_dynamico.xml, .),\ 
     40        (${MODIPSL}/modeles/DYNAMICO/xml/sponge_dynamico.xml, .),\ 
    4141        (${MODIPSL}/modeles/DYNAMICO/xml/context_dynamico.xml, .),\ 
    4242        (${MODIPSL}/modeles/DYNAMICO/xml/context_input_dynamico.xml, .),\ 

  • CONFIG/UNIFORM/v7/IPSLCM7/EXPERIMENTS/IPSLCM/piControl_TEST/COMP/dynamico.card

    r6417 r6445  
    3434        (${SUBMIT_DIR}/PARAM/run_dynamico.def_nbp${RESOL_NBP}, run_dynamico.def),\ 
    3535        (${SUBMIT_DIR}/PARAM/iodef.xml, .),\ 
    36         (${SUBMIT_DIR}/PARAM/sponge_dynamico.xml, .),\ 
     36        (${MODIPSL}/modeles/DYNAMICO/xml/sponge_dynamico.xml, .),\ 
    3737        (${MODIPSL}/modeles/DYNAMICO/xml/context_dynamico.xml, .),\ 
    3838        (${MODIPSL}/modeles/DYNAMICO/xml/context_input_dynamico.xml, .),\ 

  • CONFIG/UNIFORM/v7/IPSLCM7/SOURCES/LMDZ/oasis.F90

    r6296 r6445  
    6262  INTEGER, PARAMETER :: ids_qraiic = 28 
    6363  INTEGER, PARAMETER :: ids_qsnoic = 29 
    64   INTEGER, PARAMETER :: ids_delta_sst = 30, ids_delta_sal = 31 
    65    
    66   INTEGER, PARAMETER :: maxsend    = 31  ! Maximum number of fields to send 
     64  INTEGER, PARAMETER :: ids_delta_sst = 30, ids_delta_sal = 31, ids_dter = 32, & 
     65       ids_dser = 33, ids_dt_ds = 34 
     66   
     67  INTEGER, PARAMETER :: maxsend    = 34  ! Maximum number of fields to send 
    6768   
    6869  ! Id for fields received from ocean 
    … …  
    207208       infosend(ids_delta_sal)%action = .TRUE. 
    208209       infosend(ids_delta_sal)%name = 'CODELSSS' 
     210       infosend(ids_dter)%action = .TRUE. 
     211       infosend(ids_dter)%name = 'CODELTER' 
     212       infosend(ids_dser)%action = .TRUE. 
     213       infosend(ids_dser)%name = 'CODELSER' 
     214       infosend(ids_dt_ds)%action = .TRUE. 
     215       infosend(ids_dt_ds)%name = 'CODTDS' 
    209216    end if 
    210217            

  • CONFIG/UNIFORM/v7/IPSLCM7/SOURCES/LMDZ/physiq_mod.F90

    r6296 r6445  
    3434    USE FLOTT_GWD_rando_m, only: FLOTT_GWD_rando 
    3535    USE fonte_neige_mod, ONLY  : fonte_neige_get_vars 
    36     USE geometry_mod, ONLY: cell_area, latitude_deg, longitude_deg, longitude,latitude, & 
    37          boundslon,boundslat, dx, dy, ind_cell_glo 
     36    USE geometry_mod, ONLY: cell_area, latitude_deg, longitude_deg 
    3837    USE ioipsl, only: histbeg, histvert, histdef, histend, histsync, & 
    3938         histwrite, ju2ymds, ymds2ju, getin 
    4039    USE ioipsl_getin_p_mod, ONLY : getin_p 
    4140    USE indice_sol_mod 
    42     USE infotrac_phy, ONLY: nqtot, nbtr, nqo, tracers, types_trac, nqCO2 
     41    USE infotrac_phy, ONLY: nqtot, nbtr, nqo, tracers, type_trac 
    4342    USE readTracFiles_mod, ONLY: addPhase 
    4443    USE strings_mod,  ONLY: strIdx 
    … …  
    5453    USE pbl_surface_mod, ONLY : pbl_surface 
    5554    USE phyaqua_mod, only: zenang_an 
     55    USE phyetat0_mod, only: phyetat0 
    5656    USE phystokenc_mod, ONLY: offline, phystokenc 
    5757    USE phys_cal_mod, only: year_len, mth_len, days_elapsed, jh_1jan, & 
    … …  
    6969    USE regr_horiz_time_climoz_m, ONLY: regr_horiz_time_climoz 
    7070    USE regr_pr_time_av_m, only: regr_pr_time_av 
    71     USE surface_data,     ONLY : type_ocean, ok_veget, landice_opt 
    72     USE time_phylmdz_mod, only: annee_ref, current_time, day_ini, day_ref, & 
    73           day_step_phy, itau_phy, pdtphys, raz_date, start_time, update_time, ndays 
     71    USE surface_data,     ONLY : type_ocean, ok_veget 
     72    USE time_phylmdz_mod, only: current_time, itau_phy, pdtphys, raz_date, update_time 
    7473    USE tracinca_mod, ONLY: config_inca 
    7574    USE tropopause_m,     ONLY: dyn_tropopause 
    7675    USE ice_sursat_mod,  ONLY: flight_init, airplane 
    7776    USE vampir 
    78     USE VERTICAL_LAYERS_MOD, ONLY: aps,bps, ap, bp 
    7977    USE write_field_phy 
    8078#ifdef CPP_XIOS 
    … …  
    8381    USE lscp_mod, ONLY : lscp 
    8482    USE wake_ini_mod, ONLY : wake_ini 
     83    USE yamada_ini_mod, ONLY : yamada_ini 
     84    USE atke_turbulence_ini_mod, ONLY : atke_ini 
    8585    USE thermcell_ini_mod, ONLY : thermcell_ini 
     86    USE lscp_ini_mod, ONLY : lscp_ini 
    8687 
    8788    !USE cmp_seri_mod 
    … …  
    101102 
    102103 
     104#ifdef INCA 
     105    USE geometry_mod,      ONLY: longitude, latitude, boundslon, boundslat, ind_cell_glo 
     106    USE time_phylmdz_mod,  ONLY: ndays 
     107    USE infotrac_phy,      ONLY: nqCO2 
     108#endif 
    103109#ifdef REPROBUS 
    104     USE CHEM_REP, ONLY : Init_chem_rep_xjour, & 
    105          d_q_rep,d_ql_rep,d_qi_rep,ptrop,ttrop, & 
    106          ztrop, gravit,itroprep, Z1,Z2,fac,B 
     110    USE chem_rep, ONLY: Init_chem_rep_xjour, d_q_rep, d_ql_rep, d_qi_rep, & 
     111                        ptrop, ttrop, ztrop, gravit, itroprep, Z1, Z2, fac, B 
     112#endif 
     113#if defined INCA || defined REPROBUS 
     114    USE time_phylmdz_mod,    ONLY: annee_ref, day_ini, day_ref, start_time 
     115    USE vertical_layers_mod, ONLY: aps, bps, ap, bp 
    107116#endif 
    108117 
    … …  
    110119#ifdef CPP_RRTM 
    111120    USE YOERAD, ONLY : NRADLP 
    112     USE YOESW, ONLY : RSUN 
     121!    USE YOESW, ONLY : RSUN 
    113122#endif 
    114123 
    … …  
    136145       ! Variables locales pour effectuer les appels en serie 
    137146       t_seri,q_seri,ql_seri,qs_seri,u_seri,v_seri,tr_seri,rneb_seri, & 
     147       rhcl, &         
    138148       ! Dynamic tendencies (diagnostics) 
    139149       d_t_dyn,d_q_dyn,d_ql_dyn,d_qs_dyn,d_u_dyn,d_v_dyn,d_tr_dyn,d_rneb_dyn, & 
    … …  
    146156       d_t_ajsb,d_q_ajsb, & 
    147157       d_t_ajs,d_q_ajs,d_u_ajs,d_v_ajs, & 
    148        d_t_ajs_w,d_q_ajs_w, & 
    149        d_t_ajs_x,d_q_ajs_x, & 
     158!       d_t_ajs_w,d_q_ajs_w, & 
     159!       d_t_ajs_x,d_q_ajs_x, & 
    150160       ! 
    151161       d_t_eva,d_q_eva,d_ql_eva,d_qi_eva, & 
    … …  
    160170       d_ts, & 
    161171       ! 
    162        d_t_oli,d_u_oli,d_v_oli, & 
     172!       d_t_oli,d_u_oli,d_v_oli, & 
    163173       d_t_oro,d_u_oro,d_v_oro, & 
    164174       d_t_oro_gw,d_u_oro_gw,d_v_oro_gw, & 
    … …  
    279289       sij, & 
    280290       ! 
     291       rneblsvol, & 
     292       zqsatl, zqsats, & 
     293       qclr, qcld, qss, qvc, rnebclr, rnebss, gamma_ss, & 
     294       Tcontr, qcontr, qcontr2, fcontrN, fcontrP, & 
    281295       cldemi,  & 
    282296       cldfra, cldtau, fiwc,  & 
    … …  
    491505    REAL dtadd(klon,klev) 
    492506 
    493 !#ifdef CPP_XIOS 
    494 !    TYPE(xios_context), SAVE :: g_ctx 
    495 !#endif 
    496  
    497507#ifndef CPP_XIOS 
    498508    REAL, SAVE :: missing_val=nf90_fill_real 
    … …  
    521531    ! 
    522532    ! 
    523     INTEGER debug 
    524533    INTEGER n 
    525534    !ym      INTEGER npoints 
    … …  
    578587    ! Upmost level reached by deep convection and related variable (jyg) 
    579588    ! 
    580     INTEGER izero 
     589!    INTEGER izero 
    581590    INTEGER k_upper_cv 
    582591    !------------------------------------------------------------------ 
    … …  
    748757    REAL beta_prec_fisrt(klon,klev) ! taux de conv de l'eau cond (fisrt) 
    749758    ! RomP <<< 
    750     REAL          :: calday 
    751759 
    752760    !IM cf FH pour Tiedtke 080604 
    … …  
    802810    !C      EXTERNAL o3cm      ! initialiser l'ozone 
    803811    EXTERNAL orbite    ! calculer l'orbite terrestre 
    804     EXTERNAL phyetat0  ! lire l'etat initial de la physique 
    805812    EXTERNAL phyredem  ! ecrire l'etat de redemarrage de la physique 
    806813    EXTERNAL suphel    ! initialiser certaines constantes 
    … …  
    820827    ! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 
    821828    ! 
    822     REAL rhcl(klon,klev)    ! humiditi relative ciel clair 
     829!    REAL rhcl(klon,klev)    ! humiditi relative ciel clair 
    823830    REAL dialiq(klon,klev)  ! eau liquide nuageuse 
    824831    REAL diafra(klon,klev)  ! fraction nuageuse 
    825     REAL cldliq(klon,klev)  ! eau liquide nuageuse 
     832    REAL radocond(klon,klev)  ! eau condensee nuageuse 
    826833    ! 
    827834    !XXX PB  
    … …  
    846853    REAL conv_t(klon,klev) ! convergence de la temperature(K/s) 
    847854    ! 
    848 #ifdef INCA 
    849     REAL zxsnow_dummy(klon) 
    850 #endif 
    851855    REAL zsav_tsol(klon) 
    852856    ! 
    … …  
    863867    real zqsat(klon,klev) 
    864868    ! 
    865     INTEGER i, k, iq, j, nsrf, ll, l, itr 
     869    INTEGER i, k, iq, nsrf, l, itr 
    866870    ! 
    867871    REAL t_coup 
    … …  
    972976 
    973977    REAL picefra(klon,klev) 
     978    REAL zrel_oro(klon) 
    974979    !IM cf. AM 081204 END 
    975980    ! 
    … …  
    11161121    !     climatology and the daylight climatology 
    11171122    INTEGER,SAVE :: ncid_climoz                ! NetCDF file identifier 
    1118     REAL, POINTER, SAVE :: press_cen_climoz(:) ! Pressure levels 
    1119     REAL, POINTER, SAVE :: press_edg_climoz(:) ! Edges of pressure intervals 
    1120     REAL, POINTER, SAVE :: time_climoz(:)      ! Time vector 
     1123    REAL, allocatable, SAVE :: press_cen_climoz(:) ! Pressure levels 
     1124    REAL, allocatable, SAVE :: press_edg_climoz(:) ! Edges of pressure intervals 
     1125    REAL, allocatable, SAVE :: time_climoz(:)      ! Time vector 
    11211126    CHARACTER(LEN=13), PARAMETER :: vars_climoz(2) & 
    11221127                                  = ["tro3         ","tro3_daylight"] 
    … …  
    11781183 
    11791184#ifdef INCA 
     1185    REAL :: calday, zxsnow_dummy(klon) 
    11801186    ! set de variables utilisees pour l'initialisation des valeurs provenant de INCA 
    11811187    REAL, DIMENSION(klon,klev,naero_grp,nbands) :: init_tauinca 
    … …  
    12241230#ifdef CPP_XIOS 
    12251231! switch to XIOS LMDZ physics context  
    1226 !!!!$OMP MASTER 
    1227 !!!!    WRITE(*,*)'PHYSICS XIOS Context :', g_ctx 
    1228 !!!!    CALL wxios_set_context() 
    1229 !!!!$OMP END MASTER 
    12301232    IF (.NOT. debut .AND. is_omp_master) THEN 
    12311233       CALL wxios_set_context() 
    … …  
    13031305       CALL strataer_init 
    13041306#endif 
    1305  
    1306        !!CALL flight_init 
    13071307 
    13081308       print*, '=================================================' 
    … …  
    14711471       tau_overturning_th(:)=0. 
    14721472 
    1473        IF (ANY(types_trac == 'inca') .OR. ANY(types_trac == 'inco')) THEN 
     1473       IF (ANY(type_trac == ['inca','inco'])) THEN 
    14741474          ! jg : initialisation jusqu'au ces variables sont dans restart 
    14751475          ccm(:,:,:) = 0. 
    … …  
    17501750!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 
    17511751       CALL wake_ini(rg,rd,rv,prt_level) 
     1752       CALL yamada_ini(klon,lunout,prt_level) 
     1753       CALL atke_ini(prt_level, lunout, RG, RD, RPI, RCPD) 
    17521754       CALL thermcell_ini(iflag_thermals,prt_level,tau_thermals,lunout, & 
    17531755   &    RG,RD,RCPD,RKAPPA,RLVTT,RETV) 
    1754  
     1756       IF (ok_new_lscp) then 
     1757           CALL lscp_ini(pdtphys,ok_ice_sursat) 
     1758       endif 
    17551759!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 
    17561760 
    … …  
    20382042       !c         ENDDO 
    20392043       ! 
    2040        IF (ANY(types_trac == 'inca') .OR. ANY(types_trac == 'inco')) THEN ! ModThL 
     2044       IF (ANY(type_trac == ['inca','inco'])) THEN ! ModThL 
    20412045#ifdef INCA 
    20422046          CALL VTe(VTphysiq) 
    … …  
    21202124       ENDIF 
    21212125       ! 
    2122        IF (ANY(types_trac == 'repr')) THEN 
     2126       IF (type_trac == 'repr') THEN 
    21232127#ifdef REPROBUS 
    21242128          CALL chemini_rep(  & 
    … …  
    22222226      ENDIF 
    22232227    ENDIF 
    2224 ! switch to XIOS LMDZ physics context just in case 
    2225 !$OMP MASTER 
    2226 !!!!#ifdef CPP_XIOS 
    2227 !!!!    WRITE(*,*)'PHYSICS XIOS Context :', g_ctx 
    2228 !!!!    CALL xios_set_current_context(g_ctx) 
    2229 !!!!#endif 
    2230 !$OMP END MASTER 
    2231  
    22322228    ! 
    22332229    ! 
    … …  
    22412237 
    22422238    ! Update time and other variables in Reprobus 
    2243     IF (ANY(types_trac == 'repr')) THEN 
     2239    IF (type_trac == 'repr') THEN 
    22442240#ifdef REPROBUS 
    22452241       CALL Init_chem_rep_xjour(jD_cur-jD_ref+day_ref) 
    … …  
    30323028          ! 
    30333029          !>jyg 
    3034           IF (ANY(types_trac == 'repr')) THEN 
     3030          IF (type_trac == 'repr') THEN 
    30353031             nbtr_tmp=ntra 
    30363032          ELSE 
    … …  
    36453641    ENDIF 
    36463642 
    3647     CALL lscp(phys_tstep,missing_val,paprs,pplay, & 
     3643    CALL lscp(klon,klev,phys_tstep,missing_val,paprs,pplay, & 
    36483644         t_seri, q_seri,ptconv,ratqs, & 
    3649          d_t_lsc, d_q_lsc, d_ql_lsc, d_qi_lsc, rneb, rneb_seri, &  
    3650          cldliq, picefra, rain_lsc, snow_lsc, & 
    3651          pfrac_impa, pfrac_nucl, pfrac_1nucl, & 
     3645         d_t_lsc, d_q_lsc, d_ql_lsc, d_qi_lsc, rneb, rneblsvol, rneb_seri, &  
     3646         radocond, picefra, rain_lsc, snow_lsc, & 
    36523647         frac_impa, frac_nucl, beta_prec_fisrt, & 
    36533648         prfl, psfl, rhcl,  & 
    36543649         zqasc, fraca,ztv,zpspsk,ztla,zthl,iflag_cld_th, & 
    3655          iflag_ice_thermo, ok_ice_sursat) 
     3650         iflag_ice_thermo, ok_ice_sursat, zqsatl, zqsats, & 
     3651         qclr, qcld, qss, qvc, rnebclr, rnebss, gamma_ss, & 
     3652         Tcontr, qcontr, qcontr2, fcontrN, fcontrP ) 
    36563653 
    36573654    ELSE 
    … …  
    36593656    CALL fisrtilp(phys_tstep,paprs,pplay, & 
    36603657         t_seri, q_seri,ptconv,ratqs, & 
    3661          d_t_lsc, d_q_lsc, d_ql_lsc, d_qi_lsc, rneb, cldliq, & 
     3658         d_t_lsc, d_q_lsc, d_ql_lsc, d_qi_lsc, rneb, radocond, & 
    36623659         rain_lsc, snow_lsc, & 
    36633660         pfrac_impa, pfrac_nucl, pfrac_1nucl, & 
    … …  
    37063703          cldfra(i,k) = rneb(i,k) 
    37073704          !CR: a quoi ca sert? Faut-il ajouter qs_seri? 
    3708           IF (.NOT.new_oliq) cldliq(i,k) = ql_seri(i,k) 
     3705          IF (.NOT.new_oliq) radocond(i,k) = ql_seri(i,k) 
    37093706       ENDDO 
    37103707    ENDDO 
    … …  
    37703767          DO i = 1, klon 
    37713768             IF (diafra(i,k).GT.cldfra(i,k)) THEN 
    3772                 cldliq(i,k) = dialiq(i,k) 
     3769                radocond(i,k) = dialiq(i,k) 
    37733770                cldfra(i,k) = diafra(i,k) 
    37743771             ENDIF 
    … …  
    38073804                DO i=1,klon 
    38083805                   IF (ptconv(i,k).AND.ptconvth(i,k)) THEN 
    3809                       cldliq(i,k)=cldliq(i,k)+rnebcon(i,k)*clwcon(i,k) 
     3806                      radocond(i,k)=radocond(i,k)+rnebcon(i,k)*clwcon(i,k) 
    38103807                      cldfra(i,k)=min(cldfra(i,k)+rnebcon(i,k),1.) 
    38113808                   ELSE IF (ptconv(i,k)) THEN 
    38123809                      cldfra(i,k)=rnebcon(i,k) 
    3813                       cldliq(i,k)=rnebcon(i,k)*clwcon(i,k) 
     3810                      radocond(i,k)=rnebcon(i,k)*clwcon(i,k) 
    38143811                   ENDIF 
    38153812                ENDDO 
    … …  
    38203817                DO i=1,klon 
    38213818                   cldfra(i,k)=min(cldfra(i,k)+rnebcon(i,k),1.) 
    3822                    cldliq(i,k)=cldliq(i,k)+rnebcon(i,k)*clwcon(i,k) 
     3819                   radocond(i,k)=radocond(i,k)+rnebcon(i,k)*clwcon(i,k) 
    38233820                ENDDO 
    38243821             ENDDO 
    … …  
    38383835                   IF (ptconv(i,k).AND. .NOT.ptconvth(i,k)) THEN 
    38393836                      cldfra(i,k)=rnebcon(i,k) 
    3840                       cldliq(i,k)=rnebcon(i,k)*clwcon(i,k) 
     3837                      radocond(i,k)=rnebcon(i,k)*clwcon(i,k) 
    38413838                   ENDIF 
    38423839                ENDDO 
    … …  
    38493846          ! Ancienne version 
    38503847          cldfra(:,:)=min(max(cldfra(:,:),rnebcon(:,:)),1.) 
    3851           cldliq(:,:)=cldliq(:,:)+rnebcon(:,:)*clwcon(:,:) 
     3848          radocond(:,:)=radocond(:,:)+rnebcon(:,:)*clwcon(:,:) 
    38523849       ENDIF 
    38533850 
    … …  
    38693866          DO i = 1, klon 
    38703867             IF (diafra(i,k).GT.cldfra(i,k)) THEN 
    3871                 cldliq(i,k) = dialiq(i,k) 
     3868                radocond(i,k) = dialiq(i,k) 
    38723869                cldfra(i,k) = diafra(i,k) 
    38733870             ENDIF 
    … …  
    39373934    ENDDO 
    39383935 
    3939     IF (ANY(types_trac == 'inca') .OR. ANY(types_trac == 'inco')) THEN ! ModThL 
     3936    IF (ANY(type_trac == ['inca','inco'])) THEN ! ModThL 
    39403937#ifdef INCA 
    39413938       CALL VTe(VTphysiq) 
    … …  
    39933990#endif 
    39943991    ENDIF !type_trac = inca or inco 
    3995     IF (ANY(types_trac == 'repr')) THEN 
     3992    IF (type_trac == 'repr') THEN 
    39963993#ifdef REPROBUS 
    39973994    !CALL chemtime_rep(itap+itau_phy-1, date0, dtime, itap) 
    … …  
    40194016                     tausum_aero, tau3d_aero) 
    40204017             ENDIF 
    4021           ELSE                       ! RRTM radiation 
     4018          ELSE IF (iflag_rrtm .EQ.1) THEN  ! RRTM radiation 
    40224019             IF (aerosol_couple .AND. config_inca == 'aero' ) THEN 
    40234020                abort_message='config_inca=aero et rrtm=1 impossible' 
    … …  
    40854082                ! 
    40864083             ENDIF 
     4084          ELSE IF (iflag_rrtm .EQ.2) THEN    ! ecrad RADIATION 
     4085#ifdef CPP_ECRAD 
     4086             !--climatologies or INCA aerosols 
     4087             CALL readaerosol_optic_ecrad( debut, aerosol_couple, ok_alw, ok_volcan, & 
     4088                  flag_aerosol, flag_bc_internal_mixture, itap, jD_cur-jD_ref, & 
     4089                  pdtphys, pplay, paprs, t_seri, rhcl, presnivs,  & 
     4090                  tr_seri, mass_solu_aero, mass_solu_aero_pi,  & 
     4091                  tau_aero_sw_rrtm, piz_aero_sw_rrtm, cg_aero_sw_rrtm,  & 
     4092                  tausum_aero, drytausum_aero, tau3d_aero) 
     4093#else 
     4094                abort_message='You should compile with -rad ecrad if running with iflag_rrtm=2' 
     4095                CALL abort_physic(modname,abort_message,1) 
     4096#endif 
    40874097          ENDIF 
     4098 
    40884099       ELSE   !--flag_aerosol = 0  
    40894100          tausum_aero(:,:,:) = 0. 
    … …  
    41784189          ENDIF 
    41794190          CALL newmicro (flag_aerosol, ok_cdnc, bl95_b0, bl95_b1, & 
    4180                paprs, pplay, t_seri, cldliq, picefra, cldfra, & 
     4191               paprs, pplay, t_seri, radocond, picefra, cldfra, & 
    41814192               cldtau, cldemi, cldh, cldl, cldm, cldt, cldq, & 
    41824193               flwp, fiwp, flwc, fiwc, & 
    … …  
    41864197       ELSE 
    41874198          CALL nuage (paprs, pplay, & 
    4188                t_seri, cldliq, picefra, cldfra, cldtau, cldemi, & 
     4199               t_seri, radocond, picefra, cldfra, cldtau, cldemi, & 
    41894200               cldh, cldl, cldm, cldt, cldq, & 
    41904201               ok_aie, & 
    … …  
    43784389  
    43794390#ifndef CPP_XIOS 
    4380           !--OB 30/05/2016 modified 21/10/2016 
    4381           !--here we return swaero_diag and dryaod_diag to FALSE 
    4382           !--and histdef will switch it back to TRUE if necessary 
    4383           !--this is necessary to get the right swaero at first step 
    4384           !--but only in the case of no XIOS as XIOS is covered elsewhere  
    4385           IF (debut) swaerofree_diag = .FALSE. 
    4386           IF (debut) swaero_diag = .FALSE. 
    4387           IF (debut) dryaod_diag = .FALSE. 
    4388           !--IM 15/09/2017 here we return ok_4xCO2atm to FALSE 
    4389           !--as for swaero_diag, see above 
    4390           IF (debut) ok_4xCO2atm = .FALSE. 
    4391  
    43924391          ! 
    43934392          !IM 2eme calcul radiatif pour le cas perturbe ou au moins un 
    … …  
    45494548       DO i=1,klon 
    45504549          itest(i)=0 
    4551           !        IF ((zstd(i).gt.10.0)) THEN 
    4552           IF (((zpic(i)-zmea(i)).GT.100.).AND.(zstd(i).GT.10.0)) THEN 
     4550          zrel_oro(i)=zstd(i)/(max(zsig(i),1.E-8)*sqrt(cell_area(i))) 
     4551          !zrel_oro: relative mountain height wrt relief explained by mean slope 
     4552          ! -> condition on zrel_oro can deactivate the drag on tilted planar terrains 
     4553          !    such as ice sheets (work by V. Wiener) 
     4554          ! zpmm_orodr_t and zstd_orodr_t are activation thresholds set by F. Lott to 
     4555          ! earn computation time but they are not physical. 
     4556          IF (((zpic(i)-zmea(i)).GT.zpmm_orodr_t).AND.(zstd(i).GT.zstd_orodr_t).AND.(zrel_oro(i).LE.zrel_oro_t)) THEN 
    45534557             itest(i)=1 
    45544558             igwd=igwd+1 
    … …  
    45994603       DO i=1,klon 
    46004604          itest(i)=0 
    4601           IF ((zpic(i)-zmea(i)).GT.100.) THEN 
     4605          !zrel_oro: relative mountain height wrt relief explained by mean slope 
     4606          ! -> condition on zrel_oro can deactivate the lifting on tilted planar terrains 
     4607          !    such as ice sheets (work by V. Wiener) 
     4608          zrel_oro(i)=zstd(i)/(max(zsig(i),1.E-8)*sqrt(cell_area(i))) 
     4609          IF (((zpic(i)-zmea(i)).GT.zpmm_orolf_t).AND.(zrel_oro(i).LE.zrel_oro_t)) THEN 
    46024610             itest(i)=1 
    46034611             igwd=igwd+1 
    … …  
    47894797! car on peut s'attendre a ce que les petites echelles produisent aussi de la TKE 
    47904798! Mais attention, cela ne va pas dans le sens de la conservation de l'energie!  
    4791           IF (zstd(i).GT.1.0) THEN 
     4799          IF ((zstd(i).GT.1.0) .AND.(zrel_oro(i).LE.zrel_oro_t)) THEN 
    47924800             itest(i)=1 
    47934801             igwd=igwd+1 
    … …  
    48014809       DO i=1,klon 
    48024810          itest(i)=0 
    4803         IF (((zpic(i)-zmea(i)).GT.100.).AND.(zstd(i).GT.10.0)) THEN 
     4811        IF (((zpic(i)-zmea(i)).GT.zpmm_orodr_t).AND.(zstd(i).GT.zstd_orodr_t).AND.(zrel_oro(i).LE.zrel_oro_t)) THEN 
    48044812             itest(i)=1 
    48054813             igwd=igwd+1 
    … …  
    49684976    ! 
    49694977 
    4970     IF (ANY(types_trac=='repr')) THEN 
     4978    IF (type_trac == 'repr') THEN 
    49714979!MM pas d'impact, car on recupere q_seri,tr_seri,t_seri via phys_local_var_mod 
    49724980!MM                               dans Reprobus 
    … …  
    50185026         presnivs, pphis,     pphi,     albsol1, & 
    50195027         sh_in,   ch_in,    rhcl,      cldfra,   rneb, & 
    5020          diafra,   cldliq,    itop_con, ibas_con, & 
     5028         diafra,   radocond,    itop_con, ibas_con, & 
    50215029         pmflxr,   pmflxs,    prfl,     psfl, & 
    50225030         da,       phi,       mp,       upwd, & 
    … …  
    51215129    ENDDO 
    51225130    ! 
    5123     IF (ANY(types_trac == 'inca') .OR. ANY(types_trac == 'inco')) THEN 
     5131    IF (ANY(type_trac == ['inca','inco'])) THEN 
    51245132#ifdef INCA 
    51255133       CALL VTe(VTphysiq) 
    … …  
    51455153    ENDIF 
    51465154 
    5147     IF (ANY(types_trac == 'repr')) THEN 
     5155    IF (type_trac == 'repr') THEN 
    51485156#ifdef REPROBUS 
    51495157        CALL coord_hyb_rep(paprs, pplay, aps, bps, ap, bp, cell_area) 
    … …  
    51835191    ENDDO 
    51845192    ! 
    5185     IF (nqtot > nqo) THEN 
    5186        itr = 0 
    5187        DO iq = 1, nqtot 
    5188           IF(.NOT.tracers(iq)%isInPhysics) CYCLE 
    5189           itr = itr+1 
    5190           DO  k = 1, klev 
    5191              DO  i = 1, klon 
    5192                 d_qx(i,k,iq) = ( tr_seri(i,k,itr) - qx(i,k,iq) ) / phys_tstep 
    5193              ENDDO 
     5193    ! DC: All iterations are cycled if nqtot==nqo, so no nqtot>nqo condition required 
     5194    itr = 0 
     5195    DO iq = 1, nqtot 
     5196       IF(.NOT.tracers(iq)%isInPhysics) CYCLE 
     5197       itr = itr+1 
     5198       DO  k = 1, klev 
     5199          DO  i = 1, klon 
     5200             d_qx(i,k,iq) = ( tr_seri(i,k,itr) - qx(i,k,iq) ) / phys_tstep 
    51945201          ENDDO 
    51955202       ENDDO 
    5196     ENDIF 
     5203    ENDDO 
    51975204    ! 
    51985205    !IM rajout diagnostiques bilan KP pour analyse MJO par Jun-Ichi Yano 
    … …  
    53985405         IF (read_climoz >= 1) THEN 
    53995406           IF (is_mpi_root) CALL nf95_close(ncid_climoz) 
    5400             DEALLOCATE(press_edg_climoz) ! pointer 
    5401             DEALLOCATE(press_cen_climoz) ! pointer 
     5407            DEALLOCATE(press_edg_climoz) 
     5408            DEALLOCATE(press_cen_climoz) 
    54025409         ENDIF 
    54035410        
    … …  
    54075414 
    54085415#ifdef INCA 
    5409        if (ANY(types_trac == 'inca' )) then  
     5416       if (type_trac == 'inca') then  
    54105417          IF (is_omp_master .and. grid_type==unstructured) THEN  
    54115418             CALL finalize_inca 

  • CONFIG/UNIFORM/v7/IPSLCM7/compile_ipslcm7.sh

    r6421 r6445  
    219219echo; echo "NOW COMPILE ORCHIDEE" 
    220220echo >> $outfile ; echo " NOW COMPILE ORCHIDEE"   >> $outfile  
    221 if [ -d $modipsl/config/IPSLCM7/SOURCES/ORCHIDEE ] ; then 
    222     cd $modipsl/config/IPSLCM7/SOURCES/ORCHIDEE 
    223     for dir in $( find . -mindepth 1 -maxdepth 1 -type d ); do 
    224         for file in $( find $dir -type f ! -path '*.svn*' ); do 
    225             echo cp $file $modipsl/modeles/ORCHIDEE/$file >> $outfile 
    226             cp $file $modipsl/modeles/ORCHIDEE/$file  >> $outfile 2>&1 
    227         done 
    228     done 
    229     cd $modipsl/modeles/ORCHIDEE 
    230 fi 
    231221 
    232222# Check if the compilation of ORCHIDEE was previsouly interupted prematured.  
    … …  
    361351echo; echo "NOW COMPILE DYNAMICO " 
    362352echo >> $outfile ; echo " NOW COMPILE DYNAMICO"   >> $outfile  
    363  
    364         if [ -d $modipsl/config/IPSLCM7/SOURCES/DYNAMICO ] ; then 
    365          cd $modipsl/config/IPSLCM7/SOURCES/DYNAMICO 
    366          for dir in $( find . -mindepth 1 -maxdepth 1 -type d ); do 
    367              for file in $( find $dir -type f ! -path '*.svn*' ); do 
    368                  echo cp $file $modipsl/modeles/DYNAMICO/$file >> $outfile 
    369                  cp $file $modipsl/modeles/DYNAMICO/$file  >> $outfile 2>&1 
    370              done 
    371          done 
    372          cd $modipsl/modeles/DYNAMICO 
    373         fi 
    374  
    375353 
    376354echo ./make_icosa -$optmode -with_fcm1 -parallel $parallel -external_ioipsl -with_oasis -with_xios -arch $fcm_arch -arch_path $arch_path -job 8 $full_dyna    >> $outfile  
Note: See TracChangeset for help on using the changeset viewer.