! $Id: phyetat0_mod.F90 4581 2023-06-20 18:38:17Z evignon $

MODULE phyetat0_mod

  PRIVATE

  PUBLIC :: phyetat0

CONTAINS

  USE dimphy, only: klon, zmasq, klev

  USE iophy, ONLY : init_iophy_new

  USE ocean_cpl_mod,    ONLY : ocean_cpl_init

  USE fonte_neige_mod,  ONLY : fonte_neige_init

  USE pbl_surface_mod,  ONLY : pbl_surface_init

  USE surface_data,     ONLY : type_ocean, version_ocean

  USE phyetat0_get_mod, ONLY : phyetat0_get, phyetat0_srf

  USE phys_state_var_mod, ONLY : ancien_ok, clwcon, detr_therm, phys_tstep, &

       qsol, fevap, z0m, z0h, agesno, &

       du_gwd_rando, du_gwd_front, entr_therm, f0, fm_therm, &

       falb_dir, falb_dif, prw_ancien, prlw_ancien, prsw_ancien, prbsw_ancien, &

       ftsol, pbl_tke, pctsrf, q_ancien, ql_ancien, qs_ancien, qbs_ancien, rneb_ancien, radpas, radsol, rain_fall, ratqs, &

       rnebcon, rugoro, sig1, snow_fall, bs_fall, solaire_etat0, sollw, sollwdown, &

       solsw, solswfdiff, t_ancien, u_ancien, v_ancien, w01, wake_cstar, wake_deltaq, &

       wake_deltat, wake_delta_pbl_TKE, delta_tsurf, beta_aridity, wake_fip, wake_pe, &

       wake_s, wake_dens, awake_dens, cv_gen, zgam, zmax0, zmea, zpic, zsig, &

       zstd, zthe, zval, ale_bl, ale_bl_trig, alp_bl, u10m, v10m, treedrg, &

       ale_wake, ale_bl_stat, ds_ns, dt_ns, delta_sst, delta_sal, dter, dser, &

       dt_ds, ratqs_inter

!FC

  USE geometry_mod,     ONLY: longitude_deg, latitude_deg

  USE iostart,          ONLY: close_startphy, get_field, get_var, open_startphy

  USE infotrac_phy,     ONLY: nqtot, nbtr, type_trac, tracers

  USE readTracFiles_mod,ONLY: maxlen, new2oldH2O

  USE traclmdz_mod,     ONLY: traclmdz_from_restart

  USE carbon_cycle_mod, ONLY: carbon_cycle_init, carbon_cycle_cpl, carbon_cycle_tr, carbon_cycle_rad, co2_send, RCO2_glo

  USE indice_sol_mod,   ONLY: nbsrf, is_ter, epsfra, is_lic, is_oce, is_sic

  USE ocean_slab_mod,   ONLY: nslay, tslab, seaice, tice, ocean_slab_init

  USE time_phylmdz_mod, ONLY: init_iteration, pdtphys, itau_phy

#ifdef CPP_XIOS

  USE wxios, ONLY: missing_val

#else

  use netcdf, only: missing_val => nf90_fill_real

#endif

  use config_ocean_skin_m, only: activate_ocean_skin

  IMPLICIT none

  !======================================================================

  ! Auteur(s) Z.X. Li (LMD/CNRS) date: 19930818

  ! Objet: Lecture de l'etat initial pour la physique

  !======================================================================

  include "dimsoil.h"

  include "clesphys.h"

  include "alpale.h"

  include "compbl.h"

  include "YOMCST.h"

  !======================================================================

  CHARACTER*(*) fichnom

  ! les variables globales lues dans le fichier restart

  ! pour drag arbres FC

  CHARACTER*6 ocean_in

  LOGICAL ok_veget_in

  INTEGER        longcles

  PARAMETER    ( longcles = 20 )

  REAL clesphy0( longcles )

  REAL xmin, xmax

  INTEGER nid, nvarid

  INTEGER ierr, i, nsrf, isoil , k

  INTEGER length

  PARAMETER (length=100)

  INTEGER it, iq, isw

  REAL tab_cntrl(length), tabcntr0(length)

  CHARACTER*7 str7

  CHARACTER*2 str2

  LOGICAL :: found

  CHARACTER(LEN=maxlen) :: tname, t(2)

  ! FH1D

  !     real iolat(jjm+1)

  !real iolat(jjm+1-1/(iim*jjm))

  ! Ouvrir le fichier contenant l'etat initial:

  ! Lecture des parametres de controle:

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

  ! FH 2008/05/09 On elimine toutes les clefs physiques dans la dynamique

  ! Les constantes de la physiques sont lues dans la physique seulement.

  ! Les egalites du type

  !             tab_cntrl( 5 )=clesphy0(1)

  ! sont remplacees par

  !             clesphy0(1)=tab_cntrl( 5 )

  ! On inverse aussi la logique.

  ! On remplit les tab_cntrl avec les parametres lus dans les .def

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

  ENDDO

  ! co2_ppm : value from the previous time step

  ! co2_ppm0 : initial value of atmospheric CO2 (from create_etat0_limit.e .def)

  ! when no initial value is available e.g., from a restart

  ! this variable must be set  in a .def file which will then be

  ! used by the conf_phys_m.F90 routine.

  ! co2_ppm0 = 284.32 (illustrative example on how to set the variable in .def

  ! file, for a pre-industrial CO2 concentration value)

       ELSE

     ENDIF

     ! ELSE : keep value from .def

  ENDIF

  ! set time iteration

  ! read latitudes and make a sanity check (because already known from dyn)

      ! This is presumably serious enough to abort run

    ENDIF

    ENDIF

  ENDDO

  ! read longitudes and make a sanity check (because already known from dyn)

        ! This is presumably serious enough to abort run

      ENDIF

    ENDIF

      ENDIF

    ENDIF

  ENDDO

  ! Lecture du masque terre mer

  ENDIF

  ! Lecture des fractions pour chaque sous-surface

  ! initialisation des sous-surfaces

  ! fraction de terre

  ! fraction de glace de terre

  ! fraction d'ocean

  ! fraction glace de mer

  !  Verification de l'adequation entre le masque et les sous-surfaces

  fractint( 1 : klon) = pctsrf(1 : klon, is_ter)  &

     ENDIF

  ENDDO

  fractint (1 : klon) =  pctsrf(1 : klon, is_oce)  &

     ENDIF

  ENDDO

!===================================================================

! Lecture des temperatures du sol:

!===================================================================

     ENDDO

  ELSE

  ENDIF

!===================================================================

  ! Lecture des albedo difus et direct

!===================================================================

        ENDIF

     ENDDO

  ENDDO

!===================================================================

! Lecture dans le cas iflag_pbl_surface =1

!===================================================================

!===================================================================

  ! Lecture des temperatures du sol profond:

!===================================================================

        IF (isoil.GT.99) THEN

           PRINT*, "Trop de couches "

           call abort_physic("phyetat0", "", 1)

        ENDIF

        ENDIF

   ENDDO

!=======================================================================

! Lecture precipitation/evaporation

!=======================================================================

  ELSE

  ENDIF

!=======================================================================

! Radiation

!=======================================================================

     ENDDO

  ENDIF

  ! Lecture de la longueur de rugosite

  ELSE

  ENDIF

!FC

!CALL get_field("FTER", pctsrf(:, is_ter), found)

!  found=phyetat0_srf(treedrg,"treedrg","drag from vegetation" , 0.)

    !lecture du profile de freinage des arbres

    ELSE

!     found=phyetat0_get(treedrg,"treedrg","freinage arbres",0.)

    ENDIF

  ELSE

    ! initialize treedrg(), because it will be written in restartphy.nc

  ENDIF

  endif ! iflag_physiq <= 1

  ! Lecture de l'age de la neige:

  ! cas specifique des variables de la neige soufflee

  ELSE

  ENDIF

  ! Ehouarn: addtional tests to check if t_ancien, q_ancien contain

  !          dummy values (as is the case when generated by ce0l,

  !          or by iniaqua)

  IF ( (maxval(q_ancien).EQ.minval(q_ancien))       .OR. &

       (maxval(ql_ancien).EQ.minval(ql_ancien))     .OR. &

       (maxval(qs_ancien).EQ.minval(qs_ancien))     .OR. &

       (maxval(rneb_ancien).EQ.minval(rneb_ancien)) .OR. &

       (maxval(prw_ancien).EQ.minval(prw_ancien))   .OR. &

       (maxval(prlw_ancien).EQ.minval(prlw_ancien)) .OR. &

       (maxval(t_ancien).EQ.minval(t_ancien)) ) THEN

  ENDIF

         (maxval(prbsw_ancien).EQ.minval(prbsw_ancien)) ) THEN

    ENDIF

  ENDIF

!==================================

!  TKE

!==================================

!

  ENDIF

!!    found=phyetat0_srf(delta_tsurf,"DELTA_TSURF","Delta Ts wk/env ",0.)

!!    found=phyetat0_srf(beta_aridity,"BETA_S","Aridity factor ",1.)

  ENDIF   !(iflag_pbl>1 .AND. iflag_wake>=1 .AND. iflag_pbl_split >=1 )

!==================================

!  thermiques, poches, convection

!==================================

! Emanuel

! Wake

!jyg<

!  Set wake_dens to -1000. when there is no restart so that the actual

!  initialization is made in calwake.

!!  found=phyetat0_get(1,wake_dens,"WAKE_DENS","Wake num. /unit area",0.)

!>jyg

! Thermiques

! ALE/ALP

! fisrtilp/Clouds 0.002 could be ratqsbas. But can stay like this as well

!===========================================

  ! Read and send field trs to traclmdz

!===========================================

!--OB now this is for co2i - ThL: and therefore also for inco

     ENDIF

     it = 0

     END DO

  ENDIF

!===========================================

!  ondes de gravite / relief

!===========================================

!  ondes de gravite non orographiques

!  prise en compte du relief sous-maille

!===========================================

! Initialize ocean

!===========================================

      ELSE

          ENDDO

      ENDIF

          ENDDO

      ENDIF

      ! Sea ice variables

          ENDIF

                  seaice=917.

              ENDWHERE

          ENDIF

      ENDIF !sea ice INT

  ENDIF ! Slab

        found = phyetat0_get(delta_sal, "delta_sal", &

        found = phyetat0_get(delta_sst, "delta_SST", &

        found = phyetat0_get(dter, "dter", &

        found = phyetat0_get(dser, "dser", &

           delta_sst = missing_val

           delta_sal = missing_val

           dter = missing_val

           dser = missing_val

           dt_ds = missing_val

        end where

     end if

     found = phyetat0_get(dt_ns, "dT_ns", "delta temperature near surface", &

        ds_ns = missing_val

        dt_ns = missing_val

        delta_sst = missing_val

        delta_sal = missing_val

     end where

  end if

  ! on ferme le fichier

  ! Initialize module pbl_surface_mod

  endif

  ! Initialize module ocean_cpl_mod for the case of coupled ocean

  ENDIF

!  CALL init_iophy_new(latitude_deg, longitude_deg)

  ! Initilialize module fonte_neige_mod

END MODULE phyetat0_mod