Changeset 4741 for LMDZ6/branches


Ignore:
Timestamp:
Oct 27, 2023, 6:49:58 PM (7 months ago)
Author:
idelkadi
Message:

Initialization of the climatic fields for the first call to the COSP simulator. This, to solve the crash of LMDZ in debug mode when we activate COSP. This first call is only used to define the vertical axes for the COSP output files.

File:
1 edited

Legend:

Unmodified
Added
Removed

  • LMDZ6/branches/LMDZ_ECRad/libf/phylmd/physiq_mod.F90

    r4727 r4741  
    17771777      IF (.NOT. create_etat0_limit) CALL init_readaerosolstrato(flag_aerosol_strat)  !! initialise aero strato from file for XIOS interpolation (unstructured_grid)
    17781778
    1779       if (ok_cosp) then
     1779      IF (ok_cosp) THEN
    17801780#ifdef CPP_COSP
    17811781        ! A.I : Initialisations pour le 1er passage a Cosp
    … …  
    17891789               ecrit_mth,ecrit_day,ecrit_hf, ok_all_xml, missing_val, &
    17901790               klon,klev,longitude_deg,latitude_deg,presnivs,overlap, &
    1791                JrNt,ref_liq,ref_ice, &
    1792                pctsrf(:,is_ter)+pctsrf(:,is_lic), &
    1793                zu10m,zv10m,pphis, &
    1794                zphi,paprs(:,1:klev),pplay,zxtsol,t_seri, &
    1795                qx(:,:,ivap),zx_rh,cldfra,rnebcon,flwc,fiwc, &
    1796                prfl(:,1:klev),psfl(:,1:klev), &
    1797                pmflxr(:,1:klev),pmflxs(:,1:klev), &
    1798                mr_ozone,cldtau, cldemi)
     1791               JrNt_cosp0,ref_liq_cosp0,ref_ice_cosp0, &
     1792               pctsrf_cosp0, &
     1793               zu10m_cosp0,zv10m_cosp0,pphis, &
     1794               pphi,paprs(:,1:klev),pplay,zxtsol_cosp0,t, &
     1795               qx(:,:,ivap),zx_rh_cosp0,cldfra_cosp0,rnebcon_cosp0,flwc_cosp0,fiwc_cosp0, &
     1796               prfl_cosp0(:,1:klev),psfl_cosp0(:,1:klev), &
     1797               pmflxr_cosp0(:,1:klev),pmflxs_cosp0(:,1:klev), &
     1798               mr_ozone_cosp0,cldtau_cosp0, cldemi_cosp0)
     1799
    17991800#endif
    18001801
    18011802#ifdef CPP_COSP2
    1802           CALL phys_cosp2(itap,phys_tstep,freq_cosp, &
     1803              ! A.I : Initialisations pour le 1er passage a Cosp
     1804        CALL ini_COSP(ref_liq_cosp0,ref_ice_cosp0,pctsrf_cosp0,zu10m_cosp0,zv10m_cosp0, &
     1805               zxtsol_cosp0,zx_rh_cosp0,cldfra_cosp0,rnebcon_cosp0,flwc_cosp0, &
     1806               fiwc_cosp0,prfl_cosp0,psfl_cosp0,pmflxr_cosp0,pmflxs_cosp0, &
     1807               mr_ozone_cosp0,cldtau_cosp0,cldemi_cosp0,JrNt_cosp0)
     1808
     1809        CALL phys_cosp2(itap,phys_tstep,freq_cosp, &
    18031810               ok_mensuelCOSP,ok_journeCOSP,ok_hfCOSP, &
    18041811               ecrit_mth,ecrit_day,ecrit_hf, ok_all_xml, missing_val, &
    18051812               klon,klev,longitude_deg,latitude_deg,presnivs,overlap, &
    1806                JrNt,ref_liq,ref_ice, &
    1807                pctsrf(:,is_ter)+pctsrf(:,is_lic), &
    1808                zu10m,zv10m,pphis, &
    1809                zphi,paprs(:,1:klev),pplay,zxtsol,t_seri, &
    1810                qx(:,:,ivap),zx_rh,cldfra,rnebcon,flwc,fiwc, &
    1811                prfl(:,1:klev),psfl(:,1:klev), &
    1812                pmflxr(:,1:klev),pmflxs(:,1:klev), &
    1813                mr_ozone,cldtau, cldemi)
     1813               JrNt_cosp0,ref_liq_cosp0,ref_ice_cosp0, &
     1814               pctsrf_cosp0, &
     1815               zu10m_cosp0,zv10m_cosp0,pphis, &
     1816               pphi,paprs(:,1:klev),pplay,zxtsol_cosp0,t, &
     1817               qx(:,:,ivap),zx_rh_cosp0,cldfra_cosp0,rnebcon_cosp0,flwc_cosp0,fiwc_cosp0, &
     1818               prfl_cosp0(:,1:klev),psfl_cosp0(:,1:klev), &
     1819               pmflxr_cosp0(:,1:klev),pmflxs_cosp0(:,1:klev), &
     1820               mr_ozone_cosp0,cldtau_cosp0, cldemi_cosp0)       
    18141821#endif
    18151822
    18161823#ifdef CPP_COSPV2
    1817           CALL lmdz_cosp_interface(itap,phys_tstep,freq_cosp, &
     1824                 ! A.I : Initialisations pour le 1er passage a Cosp
     1825        CALL ini_COSP(ref_liq_cosp0,ref_ice_cosp0,pctsrf_cosp0,zu10m_cosp0,zv10m_cosp0, &
     1826               zxtsol_cosp0,zx_rh_cosp0,cldfra_cosp0,rnebcon_cosp0,flwc_cosp0, &
     1827               fiwc_cosp0,prfl_cosp0,psfl_cosp0,pmflxr_cosp0,pmflxs_cosp0, &
     1828               mr_ozone_cosp0,cldtau_cosp0,cldemi_cosp0,JrNt_cosp0)
     1829
     1830        CALL lmdz_cosp_interface(itap,phys_tstep,freq_cosp, &
    18181831               ok_mensuelCOSP,ok_journeCOSP,ok_hfCOSP, &
    18191832               ecrit_mth,ecrit_day,ecrit_hf, ok_all_xml, missing_val, &
    18201833               klon,klev,longitude_deg,latitude_deg,presnivs,overlap, &
    1821                JrNt,ref_liq,ref_ice, &
    1822                pctsrf(:,is_ter)+pctsrf(:,is_lic), &
    1823                zu10m,zv10m,pphis, &
    1824                phicosp,paprs(:,1:klev),pplay,zxtsol,t_seri, &
    1825                qx(:,:,ivap),zx_rh,cldfra,rnebcon,flwc,fiwc, &
    1826                prfl(:,1:klev),psfl(:,1:klev), &
    1827                pmflxr(:,1:klev),pmflxs(:,1:klev), &
    1828                mr_ozone,cldtau, cldemi)
     1834               JrNt_cosp0,ref_liq_cosp0,ref_ice_cosp0, &
     1835               pctsrf_cosp0, &
     1836               zu10m_cosp0,zv10m_cosp0,pphis, &
     1837               pphi,paprs(:,1:klev),pplay,zxtsol_cosp0,t, &
     1838               qx(:,:,ivap),zx_rh_cosp0,cldfra_cosp0,rnebcon_cosp0,flwc_cosp0,fiwc_cosp0, &
     1839               prfl_cosp0(:,1:klev),psfl_cosp0(:,1:klev), &
     1840               pmflxr_cosp0(:,1:klev),pmflxs_cosp0(:,1:klev), &
     1841               mr_ozone_cosp0,cldtau_cosp0, cldemi_cosp0)               
     1842               
    18291843#endif
    18301844      ENDIF
Note: See TracChangeset for help on using the changeset viewer.