Ignore:
Timestamp:
Oct 22, 2020, 10:37:13 AM (4 years ago)
Author:
lmdz-user
Message:

Moved the opening and reading of the physics restart file before the
initialisation of the physics output as the output procedure needs to
know the timestep at the beginning of the run to initialise the calendar
correctly. That timestep is read in the physics restart file.
This should close tickets #109 and #114

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ6/trunk/libf/phylmd/physiq_mod.F90

    • Property svn:keywords changed from Author Date Id Revision to Id
    r3776 r3777  
    15351535       ! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    15361536       CALL init_iophy_new(latitude_deg,longitude_deg)
     1537       CALL create_etat0_limit_unstruct
     1538       CALL phyetat0 ("startphy.nc",clesphy0,tabcntr0)
     1539
    15371540
    15381541          !===================================================================
     
    16461649               pmflxr(:,1:klev),pmflxs(:,1:klev), &
    16471650               mr_ozone,cldtau, cldemi)
    1648       ENDIF
    1649 #endif
    16501651
    16511652#ifdef CPP_COSP2
     
    17311732#endif
    17321733       IF(read_climoz>=1 .AND. create_etat0_limit) CALL regr_horiz_time_climoz(read_climoz,ok_daily_climoz)
    1733        CALL create_etat0_limit_unstruct
    1734        CALL phyetat0 ("startphy.nc",clesphy0,tabcntr0)
    17351734
    17361735!jyg<
     
    19131912
    19141913#ifdef CPP_XIOS
    1915 ! Need to put this initialisation after phyetat0 as in the coupled model the XIOS context is only
    1916 ! initialised at that moment
    1917        ! Get "missing_val" value from XML files (from temperature variable)
    1918        !$OMP MASTER
    1919        CALL xios_get_field_attr("temp",default_value=missing_val_omp)
    1920        !$OMP END MASTER
    1921        !$OMP BARRIER
    1922        missing_val=missing_val_omp
     1914       !--setting up swaero_diag to TRUE in XIOS case
     1915       IF (xios_field_is_active("topswad").OR.xios_field_is_active("topswad0").OR. &
     1916           xios_field_is_active("solswad").OR.xios_field_is_active("solswad0").OR. &
     1917           xios_field_is_active("topswai").OR.xios_field_is_active("solswai").OR.  &
     1918             (iflag_rrtm==1.AND.(xios_field_is_active("toplwad").OR.xios_field_is_active("toplwad0").OR. &
     1919                                 xios_field_is_active("sollwad").OR.xios_field_is_active("sollwad0"))))  &
     1920           !!!--for now these fields are not in the XML files so they are omitted
     1921           !!!  xios_field_is_active("toplwai").OR.xios_field_is_active("sollwai") !))) &
     1922           swaero_diag=.TRUE.
     1923
     1924       !--setting up dryaod_diag to TRUE in XIOS case
     1925       DO naero = 1, naero_tot-1
     1926         IF (xios_field_is_active("dryod550_"//name_aero_tau(naero))) dryaod_diag=.TRUE.
     1927       ENDDO
     1928       !
     1929       !--setting up ok_4xCO2atm to TRUE in XIOS case
     1930       IF (xios_field_is_active("rsut4co2").OR.xios_field_is_active("rlut4co2").OR. &
     1931           xios_field_is_active("rsutcs4co2").OR.xios_field_is_active("rlutcs4co2").OR. &
     1932           xios_field_is_active("rsu4co2").OR.xios_field_is_active("rsucs4co2").OR. &
     1933           xios_field_is_active("rsd4co2").OR.xios_field_is_active("rsdcs4co2").OR. &
     1934           xios_field_is_active("rlu4co2").OR.xios_field_is_active("rlucs4co2").OR. &
     1935           xios_field_is_active("rld4co2").OR.xios_field_is_active("rldcs4co2")) &
     1936           ok_4xCO2atm=.TRUE.
    19231937#endif
    19241938
Note: See TracChangeset for help on using the changeset viewer.