source: LMDZ5/branches/testing/libf/phylmd/phyetat0.F90 @ 2641

Last change on this file since 2641 was 2641, checked in by Laurent Fairhead, 8 years ago

Merged trunk changes r2593:2640 into testing branch

  • Property copyright set to
    Name of program: LMDZ
    Creation date: 1984
    Version: LMDZ5
    License: CeCILL version 2
    Holder: Laboratoire de m\'et\'eorologie dynamique, CNRS, UMR 8539
    See the license file in the root directory
  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 19.9 KB
RevLine 
[1403]1! $Id: phyetat0.F90 2641 2016-09-29 21:26:46Z fairhead $
[782]2
[1827]3SUBROUTINE phyetat0 (fichnom, clesphy0, tabcntr0)
[1279]4
[2073]5  USE dimphy, only: klon, zmasq, klev, nslay
[1999]6  USE iophy, ONLY : init_iophy_new
[1827]7  USE ocean_cpl_mod,    ONLY : ocean_cpl_init
8  USE fonte_neige_mod,  ONLY : fonte_neige_init
9  USE pbl_surface_mod,  ONLY : pbl_surface_init
[2220]10  USE surface_data,     ONLY : type_ocean, version_ocean
[1999]11  USE phys_state_var_mod, ONLY : ancien_ok, clwcon, detr_therm, dtime, &
[2298]12       qsol, fevap, z0m, z0h, agesno, &
[2408]13       du_gwd_rando, du_gwd_front, entr_therm, f0, fm_therm, &
[2542]14       falb_dir, falb_dif, prw_ancien, prlw_ancien, prsw_ancien, &
15       ftsol, pbl_tke, pctsrf, q_ancien, ql_ancien, qs_ancien, radpas, radsol, rain_fall, ratqs, &
[2435]16       rnebcon, rugoro, sig1, snow_fall, solaire_etat0, sollw, sollwdown, &
[1999]17       solsw, t_ancien, u_ancien, v_ancien, w01, wake_cstar, wake_deltaq, &
[2187]18       wake_deltat, wake_delta_pbl_TKE, delta_tsurf, wake_fip, wake_pe, &
[2641]19       wake_s, wake_dens, zgam, zmax0, zmea, zpic, zsig, &
[2594]20       zstd, zthe, zval, ale_bl, ale_bl_trig, alp_bl, u10m, v10m
[2435]21  USE geometry_mod, ONLY : longitude_deg, latitude_deg
[1999]22  USE iostart, ONLY : close_startphy, get_field, get_var, open_startphy
[2408]23  USE infotrac_phy, only: nbtr, nqo, type_trac, tname, niadv
[1827]24  USE traclmdz_mod,    ONLY : traclmdz_from_restart
25  USE carbon_cycle_mod, ONLY : carbon_cycle_tr, carbon_cycle_cpl, co2_send
[1999]26  USE indice_sol_mod, only: nbsrf, is_ter, epsfra, is_lic, is_oce, is_sic
[2220]27  USE ocean_slab_mod, ONLY: tslab, seaice, tice, ocean_slab_init
[2408]28  USE time_phylmdz_mod, ONLY: init_iteration, pdtphys, itau_phy
[967]29
[1827]30  IMPLICIT none
31  !======================================================================
32  ! Auteur(s) Z.X. Li (LMD/CNRS) date: 19930818
33  ! Objet: Lecture de l'etat initial pour la physique
34  !======================================================================
35  include "netcdf.inc"
36  include "dimsoil.h"
37  include "clesphys.h"
38  include "thermcell.h"
39  include "compbl.h"
[2220]40  include "YOMCST.h"
[1827]41  !======================================================================
42  CHARACTER*(*) fichnom
[524]43
[1827]44  ! les variables globales lues dans le fichier restart
[1001]45
[1827]46  REAL tsoil(klon, nsoilmx, nbsrf)
47  REAL qsurf(klon, nbsrf)
48  REAL snow(klon, nbsrf)
49  real fder(klon)
50  REAL run_off_lic_0(klon)
51  REAL fractint(klon)
52  REAL trs(klon, nbtr)
[2220]53  REAL zts(klon)
[651]54
[1827]55  CHARACTER*6 ocean_in
56  LOGICAL ok_veget_in
[879]57
[1827]58  INTEGER        longcles
59  PARAMETER    ( longcles = 20 )
60  REAL clesphy0( longcles )
[766]61
[1827]62  REAL xmin, xmax
[766]63
[1827]64  INTEGER nid, nvarid
65  INTEGER ierr, i, nsrf, isoil , k
66  INTEGER length
67  PARAMETER (length=100)
[2258]68  INTEGER it, iiq, isw
[1827]69  REAL tab_cntrl(length), tabcntr0(length)
70  CHARACTER*7 str7
71  CHARACTER*2 str2
[2298]72  LOGICAL :: found,phyetat0_get,phyetat0_srf
[2435]73  REAL :: lon_startphy(klon), lat_startphy(klon)
[1827]74
75  ! FH1D
76  !     real iolat(jjm+1)
[2408]77  !real iolat(jjm+1-1/(iim*jjm))
[1827]78
79  ! Ouvrir le fichier contenant l'etat initial:
80
81  CALL open_startphy(fichnom)
82
83  ! Lecture des parametres de controle:
84
85  CALL get_var("controle", tab_cntrl)
86
[956]87!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
[1827]88  ! FH 2008/05/09 On elimine toutes les clefs physiques dans la dynamique
89  ! Les constantes de la physiques sont lues dans la physique seulement.
90  ! Les egalites du type
91  !             tab_cntrl( 5 )=clesphy0(1)
92  ! sont remplacees par
93  !             clesphy0(1)=tab_cntrl( 5 )
94  ! On inverse aussi la logique.
95  ! On remplit les tab_cntrl avec les parametres lus dans les .def
[956]96!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
97
[1827]98  DO i = 1, length
99     tabcntr0( i ) = tab_cntrl( i )
100  ENDDO
[1279]101
[2408]102  tab_cntrl(1)=pdtphys
[1827]103  tab_cntrl(2)=radpas
[1279]104
[1827]105  ! co2_ppm : value from the previous time step
106  IF (carbon_cycle_tr .OR. carbon_cycle_cpl) THEN
107     co2_ppm = tab_cntrl(3)
108     RCO2    = co2_ppm * 1.0e-06  * 44.011/28.97
109     ! ELSE : keep value from .def
110  END IF
[1279]111
[1827]112  ! co2_ppm0 : initial value of atmospheric CO2 (from create_etat0_limit.e .def)
113  co2_ppm0   = tab_cntrl(16)
[524]114
[1827]115  solaire_etat0      = tab_cntrl(4)
116  tab_cntrl(5)=iflag_con
117  tab_cntrl(6)=nbapp_rad
[524]118
[1827]119  if (cycle_diurne) tab_cntrl( 7) =1.
120  if (soil_model) tab_cntrl( 8) =1.
121  if (new_oliq) tab_cntrl( 9) =1.
122  if (ok_orodr) tab_cntrl(10) =1.
123  if (ok_orolf) tab_cntrl(11) =1.
124  if (ok_limitvrai) tab_cntrl(12) =1.
[956]125
[1827]126  itau_phy = tab_cntrl(15)
[956]127
[1827]128  clesphy0(1)=tab_cntrl( 5 )
129  clesphy0(2)=tab_cntrl( 6 )
130  clesphy0(3)=tab_cntrl( 7 )
131  clesphy0(4)=tab_cntrl( 8 )
132  clesphy0(5)=tab_cntrl( 9 )
133  clesphy0(6)=tab_cntrl( 10 )
134  clesphy0(7)=tab_cntrl( 11 )
135  clesphy0(8)=tab_cntrl( 12 )
[956]136
[2408]137  ! set time iteration
138   CALL init_iteration(itau_phy)
139
[2435]140  ! read latitudes and make a sanity check (because already known from dyn)
141  CALL get_field("latitude",lat_startphy)
142  DO i=1,klon
143    IF (ABS(lat_startphy(i)-latitude_deg(i))>=1) THEN
144      WRITE(*,*) "phyetat0: Error! Latitude discrepancy wrt startphy file:",&
145                 " i=",i," lat_startphy(i)=",lat_startphy(i),&
146                 " latitude_deg(i)=",latitude_deg(i)
147      ! This is presumably serious enough to abort run
148      CALL abort_physic("phyetat0","discrepancy in latitudes!",1)
149    ENDIF
150    IF (ABS(lat_startphy(i)-latitude_deg(i))>=0.0001) THEN
151      WRITE(*,*) "phyetat0: Warning! Latitude discrepancy wrt startphy file:",&
152                 " i=",i," lat_startphy(i)=",lat_startphy(i),&
153                 " latitude_deg(i)=",latitude_deg(i)
154    ENDIF
155  ENDDO
[766]156
[2435]157  ! read longitudes and make a sanity check (because already known from dyn)
158  CALL get_field("longitude",lon_startphy)
159  DO i=1,klon
160    IF (ABS(lon_startphy(i)-longitude_deg(i))>=1) THEN
161      WRITE(*,*) "phyetat0: Error! Longitude discrepancy wrt startphy file:",&
162                 " i=",i," lon_startphy(i)=",lon_startphy(i),&
163                 " longitude_deg(i)=",longitude_deg(i)
164      ! This is presumably serious enough to abort run
165      CALL abort_physic("phyetat0","discrepancy in longitudes!",1)
166    ENDIF
167    IF (ABS(lon_startphy(i)-longitude_deg(i))>=0.0001) THEN
168      WRITE(*,*) "phyetat0: Warning! Longitude discrepancy wrt startphy file:",&
169                 " i=",i," lon_startphy(i)=",lon_startphy(i),&
170                 " longitude_deg(i)=",longitude_deg(i)
171    ENDIF
172  ENDDO
[1001]173
[1827]174  ! Lecture du masque terre mer
[766]175
[1827]176  CALL get_field("masque", zmasq, found)
177  IF (.NOT. found) THEN
178     PRINT*, 'phyetat0: Le champ <masque> est absent'
179     PRINT *, 'fichier startphy non compatible avec phyetat0'
180  ENDIF
[1001]181
[1827]182  ! Lecture des fractions pour chaque sous-surface
[766]183
[1827]184  ! initialisation des sous-surfaces
[766]185
[1827]186  pctsrf = 0.
[766]187
[1827]188  ! fraction de terre
[766]189
[1827]190  CALL get_field("FTER", pctsrf(:, is_ter), found)
191  IF (.NOT. found) PRINT*, 'phyetat0: Le champ <FTER> est absent'
[766]192
[1827]193  ! fraction de glace de terre
[766]194
[1827]195  CALL get_field("FLIC", pctsrf(:, is_lic), found)
196  IF (.NOT. found) PRINT*, 'phyetat0: Le champ <FLIC> est absent'
[1001]197
[1827]198  ! fraction d'ocean
[1001]199
[1827]200  CALL get_field("FOCE", pctsrf(:, is_oce), found)
201  IF (.NOT. found) PRINT*, 'phyetat0: Le champ <FOCE> est absent'
[1001]202
[1827]203  ! fraction glace de mer
[1001]204
[1827]205  CALL get_field("FSIC", pctsrf(:, is_sic), found)
206  IF (.NOT. found) PRINT*, 'phyetat0: Le champ <FSIC> est absent'
[1001]207
[1827]208  !  Verification de l'adequation entre le masque et les sous-surfaces
209
210  fractint( 1 : klon) = pctsrf(1 : klon, is_ter)  &
211       + pctsrf(1 : klon, is_lic)
212  DO i = 1 , klon
213     IF ( abs(fractint(i) - zmasq(i) ) .GT. EPSFRA ) THEN
214        WRITE(*, *) 'phyetat0: attention fraction terre pas ',  &
215             'coherente ', i, zmasq(i), pctsrf(i, is_ter) &
216             , pctsrf(i, is_lic)
217        WRITE(*, *) 'Je force la coherence zmasq=fractint'
218        zmasq(i) = fractint(i)
219     ENDIF
220  END DO
221  fractint (1 : klon) =  pctsrf(1 : klon, is_oce)  &
222       + pctsrf(1 : klon, is_sic)
223  DO i = 1 , klon
224     IF ( abs( fractint(i) - (1. - zmasq(i))) .GT. EPSFRA ) THEN
225        WRITE(*, *) 'phyetat0 attention fraction ocean pas ',  &
226             'coherente ', i, zmasq(i) , pctsrf(i, is_oce) &
227             , pctsrf(i, is_sic)
[2056]228        WRITE(*, *) 'Je force la coherence zmasq=1.-fractint'
229        zmasq(i) = 1. - fractint(i)
[1827]230     ENDIF
231  END DO
232
[2298]233!===================================================================
234! Lecture des temperatures du sol:
235!===================================================================
[1827]236
[2298]237  found=phyetat0_get(1,ftsol(:,1),"TS","Surface temperature",283.)
238  IF (found) THEN
239     DO nsrf=2,nbsrf
240        ftsol(:,nsrf)=ftsol(:,1)
[1827]241     ENDDO
242  ELSE
[2298]243     found=phyetat0_srf(1,ftsol,"TS","Surface temperature",283.)
[1827]244  ENDIF
[524]245
[2258]246!===================================================================
247  ! Lecture des albedo difus et direct
[2298]248!===================================================================
[2258]249
250  DO nsrf = 1, nbsrf
251     DO isw=1, nsw
[2298]252        IF (isw.GT.99) THEN
253           PRINT*, "Trop de bandes SW"
[2408]254           call abort_physic("phyetat0", "", 1)
[2258]255        ENDIF
[2298]256        WRITE(str2, '(i2.2)') isw
257        found=phyetat0_srf(1,falb_dir(:, isw,:),"A_dir_SW"//str2//"srf","Direct Albedo",0.2)
258        found=phyetat0_srf(1,falb_dif(:, isw,:),"A_dif_SW"//str2//"srf","Direct Albedo",0.2)
[2258]259     ENDDO
260  ENDDO
261
[2594]262  found=phyetat0_srf(1,u10m,"U10M","u a 10m",0.)
263  found=phyetat0_srf(1,v10m,"V10M","v a 10m",0.)
[2568]264
[2258]265!===================================================================
[1827]266  ! Lecture des temperatures du sol profond:
[2298]267!===================================================================
[524]268
[2298]269   DO isoil=1, nsoilmx
270        IF (isoil.GT.99) THEN
271           PRINT*, "Trop de couches "
[2408]272           call abort_physic("phyetat0", "", 1)
[1827]273        ENDIF
[2298]274        WRITE(str2,'(i2.2)') isoil
275        found=phyetat0_srf(1,tsoil(:, isoil,:),"Tsoil"//str2//"srf","Temp soil",0.)
[1827]276        IF (.NOT. found) THEN
277           PRINT*, "phyetat0: Le champ <Tsoil"//str7//"> est absent"
278           PRINT*, "          Il prend donc la valeur de surface"
[2298]279           tsoil(:, isoil, :)=ftsol(:, :)
[1827]280        ENDIF
[2298]281   ENDDO
[524]282
[2298]283!=======================================================================
284! Lecture precipitation/evaporation
285!=======================================================================
[1001]286
[2298]287  found=phyetat0_srf(1,qsurf,"QS","Near surface hmidity",0.)
288  found=phyetat0_get(1,qsol,"QSOL","Surface hmidity / bucket",0.)
289  found=phyetat0_srf(1,snow,"SNOW","Surface snow",0.)
290  found=phyetat0_srf(1,fevap,"EVAP","evaporation",0.)
291  found=phyetat0_get(1,snow_fall,"snow_f","snow fall",0.)
292  found=phyetat0_get(1,rain_fall,"rain_f","rain fall",0.)
[1001]293
[2298]294!=======================================================================
295! Radiation
296!=======================================================================
[1001]297
[2298]298  found=phyetat0_get(1,solsw,"solsw","net SW radiation surf",0.)
299  found=phyetat0_get(1,sollw,"sollw","net LW radiation surf",0.)
300  found=phyetat0_get(1,sollwdown,"sollwdown","down LW radiation surf",0.)
[1827]301  IF (.NOT. found) THEN
[2298]302     sollwdown = 0. ;  zts=0.
[2220]303     do nsrf=1,nbsrf
304        zts(:)=zts(:)+ftsol(:,nsrf)*pctsrf(:,nsrf)
305     enddo
306     sollwdown(:)=sollw(:)+RSIGMA*zts(:)**4
307  ENDIF
308
[2298]309  found=phyetat0_get(1,radsol,"RADS","Solar radiation",0.)
310  found=phyetat0_get(1,fder,"fder","Flux derivative",0.)
[2220]311
[1827]312
313  ! Lecture de la longueur de rugosite
[2298]314  found=phyetat0_srf(1,z0m,"RUG","Z0m ancien",0.001)
315  IF (found) THEN
316     z0h(:,1:nbsrf)=z0m(:,1:nbsrf)
[1827]317  ELSE
[2298]318     found=phyetat0_srf(1,z0m,"Z0m","Roughness length, momentum ",0.001)
319     found=phyetat0_srf(1,z0h,"Z0h","Roughness length, enthalpy ",0.001)
[1827]320  ENDIF
321
322  ! Lecture de l'age de la neige:
[2298]323  found=phyetat0_srf(1,agesno,"AGESNO","SNOW AGE",0.001)
[1827]324
[2298]325  ancien_ok=.true.
326  ancien_ok=ancien_ok.AND.phyetat0_get(klev,t_ancien,"TANCIEN","TANCIEN",0.)
327  ancien_ok=ancien_ok.AND.phyetat0_get(klev,q_ancien,"QANCIEN","QANCIEN",0.)
[2542]328  ancien_ok=ancien_ok.AND.phyetat0_get(klev,ql_ancien,"QLANCIEN","QLANCIEN",0.)
329  ancien_ok=ancien_ok.AND.phyetat0_get(klev,qs_ancien,"QSANCIEN","QSANCIEN",0.)
[2298]330  ancien_ok=ancien_ok.AND.phyetat0_get(klev,u_ancien,"UANCIEN","UANCIEN",0.)
331  ancien_ok=ancien_ok.AND.phyetat0_get(klev,v_ancien,"VANCIEN","VANCIEN",0.)
[2542]332  ancien_ok=ancien_ok.AND.phyetat0_get(1,prw_ancien,"PRWANCIEN","PRWANCIEN",0.)
333  ancien_ok=ancien_ok.AND.phyetat0_get(1,prlw_ancien,"PRLWANCIEN","PRLWANCIEN",0.)
334  ancien_ok=ancien_ok.AND.phyetat0_get(1,prsw_ancien,"PRSWANCIEN","PRSWANCIEN",0.)
[1827]335
[2542]336  ! Ehouarn: addtional tests to check if t_ancien, q_ancien contain
337  !          dummy values (as is the case when generated by ce0l,
338  !          or by iniaqua)
339  if ( (maxval(q_ancien).eq.minval(q_ancien))       .or. &
340       (maxval(ql_ancien).eq.minval(ql_ancien))     .or. &
341       (maxval(qs_ancien).eq.minval(qs_ancien))     .or. &
342       (maxval(prw_ancien).eq.minval(prw_ancien))   .or. &
343       (maxval(prlw_ancien).eq.minval(prlw_ancien)) .or. &
344       (maxval(prsw_ancien).eq.minval(prsw_ancien)) .or. &
345       (maxval(t_ancien).eq.minval(t_ancien)) ) then
346    ancien_ok=.false.
347  endif
348
[2298]349  found=phyetat0_get(klev,clwcon,"CLWCON","CLWCON",0.)
350  found=phyetat0_get(klev,rnebcon,"RNEBCON","RNEBCON",0.)
351  found=phyetat0_get(klev,ratqs,"RATQS","RATQS",0.)
[1827]352
[2568]353
[2298]354  found=phyetat0_get(1,run_off_lic_0,"RUNOFFLIC0","RUNOFFLIC0",0.)
[1827]355
[2298]356!==================================
357!  TKE
358!==================================
359!
[1827]360  IF (iflag_pbl>1) then
[2298]361     found=phyetat0_srf(klev+1,pbl_tke,"TKE","Turb. Kinetic. Energ. ",1.e-8)
[1827]362  ENDIF
[1403]363
[2298]364  IF (iflag_pbl>1 .AND. iflag_wake>=1  .AND. iflag_pbl_split >=1 ) then
365    found=phyetat0_srf(klev+1,wake_delta_pbl_tke,"DELTATKE","Del TKE wk/env",0.)
366    found=phyetat0_srf(1,delta_tsurf,"DELTA_TSURF","Delta Ts wk/env ",0.)
[2187]367  ENDIF   !(iflag_pbl>1 .AND. iflag_wake>=1 .AND. iflag_pbl_split >=1 )
368
[2298]369!==================================
370!  thermiques, poches, convection
371!==================================
[1403]372
[2298]373! Emanuel
374  found=phyetat0_get(klev,sig1,"sig1","sig1",0.)
375  found=phyetat0_get(klev,w01,"w01","w01",0.)
[1403]376
[2298]377! Wake
378  found=phyetat0_get(klev,wake_deltat,"WAKE_DELTAT","Delta T wake/env",0.)
379  found=phyetat0_get(klev,wake_deltaq,"WAKE_DELTAQ","Delta hum. wake/env",0.)
[2641]380  found=phyetat0_get(1,wake_s,"WAKE_S","Wake frac. area",0.)
381  found=phyetat0_get(1,wake_dens,"WAKE_DENS","Wake num. /unit area",0.)
[2298]382  found=phyetat0_get(1,wake_cstar,"WAKE_CSTAR","WAKE_CSTAR",0.)
383  found=phyetat0_get(1,wake_pe,"WAKE_PE","WAKE_PE",0.)
384  found=phyetat0_get(1,wake_fip,"WAKE_FIP","WAKE_FIP",0.)
[1403]385
[2298]386! Thermiques
387  found=phyetat0_get(1,zmax0,"ZMAX0","ZMAX0",40.)
388  found=phyetat0_get(1,f0,"F0","F0",1.e-5)
389  found=phyetat0_get(klev+1,fm_therm,"FM_THERM","Thermals mass flux",0.)
390  found=phyetat0_get(klev,entr_therm,"ENTR_THERM","Thermals Entrain.",0.)
391  found=phyetat0_get(klev,detr_therm,"DETR_THERM","Thermals Detrain.",0.)
[1403]392
[2298]393! ALE/ALP
394  found=phyetat0_get(1,ale_bl,"ALE_BL","ALE BL",0.)
395  found=phyetat0_get(1,ale_bl_trig,"ALE_BL_TRIG","ALE BL_TRIG",0.)
396  found=phyetat0_get(1,alp_bl,"ALP_BL","ALP BL",0.)
[1403]397
[2298]398!===========================================
[1827]399  ! Read and send field trs to traclmdz
[2298]400!===========================================
[1827]401
402  IF (type_trac == 'lmdz') THEN
[2298]403     DO it=1, nbtr                                                                 
404!!        iiq=niadv(it+2)                                                           ! jyg
405        iiq=niadv(it+nqo)                                                           ! jyg
406        found=phyetat0_get(1,trs(:,it),"trs_"//tname(iiq), &
407              "Surf trac"//tname(iiq),0.)
[1827]408     END DO
409     CALL traclmdz_from_restart(trs)
410
411     IF (carbon_cycle_cpl) THEN
412        ALLOCATE(co2_send(klon), stat=ierr)
[2408]413        IF (ierr /= 0) CALL abort_physic('phyetat0', 'pb allocation co2_send', 1)
[2298]414        found=phyetat0_get(1,co2_send,"co2_send","co2 send",0.)
[1827]415     END IF
416  END IF
417
[2298]418!===========================================
419!  ondes de gravite / relief
420!===========================================
421
422!  ondes de gravite non orographiques
[2408]423  if (ok_gwd_rando) found = &
424       phyetat0_get(klev,du_gwd_rando,"du_gwd_rando","du_gwd_rando",0.)
425  IF (.not. ok_hines .and. ok_gwd_rando) found &
426       = phyetat0_get(klev,du_gwd_front,"du_gwd_front","du_gwd_front",0.)
[1999]427
[2298]428!  prise en compte du relief sous-maille
429  found=phyetat0_get(1,zmea,"ZMEA","sub grid orography",0.)
430  found=phyetat0_get(1,zstd,"ZSTD","sub grid orography",0.)
431  found=phyetat0_get(1,zsig,"ZSIG","sub grid orography",0.)
432  found=phyetat0_get(1,zgam,"ZGAM","sub grid orography",0.)
433  found=phyetat0_get(1,zthe,"ZTHE","sub grid orography",0.)
434  found=phyetat0_get(1,zpic,"ZPIC","sub grid orography",0.)
435  found=phyetat0_get(1,zval,"ZVAL","sub grid orography",0.)
436  found=phyetat0_get(1,zmea,"ZMEA","sub grid orography",0.)
437  found=phyetat0_get(1,rugoro,"RUGSREL","sub grid orography",0.)
438
439!===========================================
440! Initialize ocean
441!===========================================
442
[2073]443  IF ( type_ocean == 'slab' ) THEN
[2220]444      CALL ocean_slab_init(dtime, pctsrf)
[2298]445      found=phyetat0_get(nslay,tslab,"tslab","tslab",0.)
[2073]446      IF (.NOT. found) THEN
447          PRINT*, "phyetat0: Le champ <tslab> est absent"
448          PRINT*, "Initialisation a tsol_oce"
449          DO i=1,nslay
[2220]450              tslab(:,i)=MAX(ftsol(:,is_oce),271.35)
[2073]451          END DO
452      END IF
[2298]453
[2220]454      ! Sea ice variables
[2298]455      found=phyetat0_get(1,tice,"slab_tice","slab_tice",0.)
[2220]456      IF (version_ocean == 'sicINT') THEN
457          IF (.NOT. found) THEN
458              PRINT*, "phyetat0: Le champ <tice> est absent"
459              PRINT*, "Initialisation a tsol_sic"
460                  tice(:)=ftsol(:,is_sic)
461          END IF
462          IF (.NOT. found) THEN
463              PRINT*, "phyetat0: Le champ <seaice> est absent"
464              PRINT*, "Initialisation a 0/1m suivant fraction glace"
465              seaice(:)=0.
466              WHERE (pctsrf(:,is_sic).GT.EPSFRA)
467                  seaice=917.
468              END WHERE
469          END IF
470      END IF !sea ice INT
[2073]471  END IF ! Slab       
472
[1827]473  ! on ferme le fichier
474  CALL close_startphy
475
476  ! Initialize module pbl_surface_mod
477
[2298]478  CALL pbl_surface_init(fder, snow, qsurf, tsoil)
[1827]479
480  ! Initialize module ocean_cpl_mod for the case of coupled ocean
481  IF ( type_ocean == 'couple' ) THEN
[2435]482     CALL ocean_cpl_init(dtime, longitude_deg, latitude_deg)
[1827]483  ENDIF
484
[2435]485  CALL init_iophy_new(latitude_deg, longitude_deg)
[2056]486
[1827]487  ! Initilialize module fonte_neige_mod     
488  CALL fonte_neige_init(run_off_lic_0)
489
490END SUBROUTINE phyetat0
[2298]491
492!===================================================================
493FUNCTION phyetat0_get(nlev,field,name,descr,default)
494!===================================================================
495! Lecture d'un champ avec contrôle
496! Function logique dont le resultat indique si la lecture
497! s'est bien passée
498! On donne une valeur par defaut dans le cas contraire
499!===================================================================
500
501USE iostart, ONLY : get_field
502USE dimphy, only: klon
[2408]503USE print_control_mod, ONLY: lunout
[2298]504
505IMPLICIT NONE
506
507LOGICAL phyetat0_get
508
509! arguments
510INTEGER,INTENT(IN) :: nlev
511CHARACTER*(*),INTENT(IN) :: name,descr
512REAL,INTENT(IN) :: default
513REAL,DIMENSION(klon,nlev),INTENT(INOUT) :: field
514
515! Local variables
516LOGICAL found
517
518   CALL get_field(name, field, found)
519   IF (.NOT. found) THEN
520     WRITE(lunout,*) "phyetat0: Le champ <",name,"> est absent"
521     WRITE(lunout,*) "Depart legerement fausse. Mais je continue"
522     field(:,:)=default
523   ENDIF
524   WRITE(lunout,*) name, descr, MINval(field),MAXval(field)
525   phyetat0_get=found
526
527RETURN
528END FUNCTION phyetat0_get
529
530!================================================================
531FUNCTION phyetat0_srf(nlev,field,name,descr,default)
532!===================================================================
533! Lecture d'un champ par sous-surface avec contrôle
534! Function logique dont le resultat indique si la lecture
535! s'est bien passée
536! On donne une valeur par defaut dans le cas contraire
537!===================================================================
538
539USE iostart, ONLY : get_field
540USE dimphy, only: klon
541USE indice_sol_mod, only: nbsrf
[2408]542USE print_control_mod, ONLY: lunout
[2298]543
544IMPLICIT NONE
545
546LOGICAL phyetat0_srf
547! arguments
548INTEGER,INTENT(IN) :: nlev
549CHARACTER*(*),INTENT(IN) :: name,descr
550REAL,INTENT(IN) :: default
551REAL,DIMENSION(klon,nlev,nbsrf),INTENT(INOUT) :: field
552
553! Local variables
554LOGICAL found,phyetat0_get
555INTEGER nsrf
556CHARACTER*2 str2
557 
558     IF (nbsrf.GT.99) THEN
559        WRITE(lunout,*) "Trop de sous-mailles"
[2408]560        call abort_physic("phyetat0", "", 1)
[2298]561     ENDIF
562
563     DO nsrf = 1, nbsrf
564        WRITE(str2, '(i2.2)') nsrf
565        found= phyetat0_get(nlev,field(:,:, nsrf), &
566        name//str2,descr//" srf:"//str2,default)
567     ENDDO
568
569     phyetat0_srf=found
570
571RETURN
572END FUNCTION phyetat0_srf
573
Note: See TracBrowser for help on using the repository browser.