Ignore:
Timestamp:
Dec 6, 2022, 12:01:16 AM (22 months ago)
Author:
lguez
Message:

Sync latest trunk changes to Ocean_skin

Location:
LMDZ6/branches/Ocean_skin
Files:
2 edited

Legend:

Unmodified
Added
Removed

  • LMDZ6/branches/Ocean_skin

  • LMDZ6/branches/Ocean_skin/libf/phylmdiso/isotopes_routines_mod.F90

    r3927 r4368  
    33
    44MODULE isotopes_routines_mod
     5  USE infotrac_phy, ONLY: niso, ntraciso=>ntiso, index_trac=>itZonIso, ntraceurs_zone=>nzone
    56IMPLICIT NONE
    67
    … …  
    1314&            zqs,zq_ancien,zqev_diag,zq)
    1415
    15 USE infotrac_phy, ONLY: ntraciso,niso, &
    16         ntraceurs_zone,index_trac
    1716USE isotopes_mod, ONLY: ridicule, ridicule_rain, thumxt1, no_pce,  &
    1817&       bidouille_anti_divergence, &
    … …  
    846845&    L, xtnu,Pveg)
    847846
    848 USE infotrac_phy, ONLY: niso
    849847USE isotopes_mod, ONLY: ridicule_qsol, ridicule, &
    850848&       ridicule_evap,P_veg,iso_HDO,iso_eau,iso_O17,iso_O18
    … …  
    13011299
    13021300subroutine calcul_kcin(Vsurf,KCIN)
    1303 USE infotrac_phy, ONLY: niso
    13041301USE isotopes_mod, ONLY: tv0cin,tkcin0,tkcin1,tkcin2
    13051302implicit none
    … …  
    13281325
    13291326     subroutine fractcalk(kt, ptin, pxtfra, pfraice)
    1330 !USE infotrac_phy, ONLY: use_iso
    13311327USE isotopes_mod, ONLY: talph1,talph2,talph3,pxtmin,iso_O17, &
    13321328&       fac_coeff_eq17_liq, pxtmelt, &
    … …  
    14571453      subroutine fractcalk_liq(kt, ptin, pxtfra)
    14581454
    1459 !      USE infotrac_phy, ONLY: use_iso
    14601455      USE isotopes_mod, ONLY: pxtmin,talph1,talph2,talph3, &
    14611456&       fac_coeff_eq17_liq, pxtice, &
    … …  
    15221517      subroutine fractcalk_glace(kt, ptin, pfraice)
    15231518
    1524 !      use infotrac_phy, ONLY: use_iso
    15251519      use isotopes_mod, ONLY: talps1,talps2, iso_O17,fac_coeff_eq17_ice, &
    15261520        & pxtmelt,musi, lambda_sursat, tdifrel, &
    … …  
    16311625      subroutine fractcalk_vectall(ptin, pxtfra, pfraice,n)
    16321626
    1633         USE infotrac_phy, ONLY: niso
    16341627        USE isotopes_mod, ONLY: talph1,talph2,talph3,tdifrel,pxtmin, &
    16351628&      iso_O17, iso_HTO, iso_eau, iso_O18, iso_HDO, musi, lambda_sursat, &
    … …  
    18031796      subroutine fractcalk_vectall_liq(ptin, pxtfra, n)
    18041797
    1805       USE infotrac_phy, ONLY: niso
    18061798      USE isotopes_mod, ONLY: pxtmin,talph1,talph2,talph3, &
    18071799&       iso_eau,iso_HDO, iso_O18, iso_O17,iso_HTO,fac_coeff_eq17_liq, &
    … …  
    18821874      subroutine fractcalk_vectall_ice(ptin, pfraice,n)
    18831875
    1884       use infotrac_phy, ONLY: niso
    18851876      use isotopes_mod, ONLY: talps1,talps2, fac_coeff_eq17_ice, &
    18861877        & pxtmelt,musi, lambda_sursat, tdifrel, &
    … …  
    20232014&            i,Rsol,klon)
    20242015
    2025   USE infotrac_phy, ONLY: niso,ntraciso
    20262016  USE isotopes_mod, ONLY: iso_eau, iso_HDO,ridicule, &
    20272017&        ridicule_qsol,iso_O17,iso_O18
    … …  
    22332223&          i,xtevap,klon) 
    22342224
    2235   USE infotrac_phy, ONLY: ntraciso,niso
    22362225  USE isotopes_mod, ONLY: iso_eau, iso_HDO,ridicule,ridicule_rain, &
    22372226        iso_O18,iso_O17
    … …  
    24442433&   )
    24452434
    2446   USE infotrac_phy, ONLY: ntraciso,niso
    24472435  USE isotopes_mod, ONLY: iso_eau, iso_HDO,expb_max,tdifrel,tdifexp, &
    24482436&       ridicule,thumxt1,ridicule_rain,bidouille_anti_divergence, &
    … …  
    45004488&           Tevap)
    45014489
    4502   USE infotrac_phy, ONLY: niso,ntraciso
    45034490  USE isotopes_mod, ONLY: iso_eau, iso_HDO,bidouille_anti_divergence, &
    45044491&       ridicule,ridicule_rain
    … …  
    46584645&           ,fac_ftmr)
    46594646
    4660   USE infotrac_phy, ONLY: niso,ntraciso
    46614647  USE isotopes_mod, ONLY: iso_eau, iso_HDO,bidouille_anti_divergence, &
    46624648&       Rdefault,ridicule,ridicule_rain
    … …  
    49044890     &           Pqiinf_cas,Pqiinf)
    49054891
    4906   USE infotrac_phy, ONLY: niso,ntraciso
    49074892  USE isotopes_mod, ONLY: iso_eau, iso_HDO
    49084893
    … …  
    50665051     &           xtnew_cas,xtnew,Pxtiinf_cas,Pxtiinf)
    50675052
    5068   USE infotrac_phy, ONLY: niso
    50695053  USE isotopes_mod, ONLY: iso_eau, iso_HDO
    50705054#ifdef ISOVERIF
    … …  
    51115095     &           ncum)
    51125096
    5113   USE infotrac_phy, ONLY: niso,ntraciso
    51145097  USE isotopes_mod, ONLY: iso_eau, iso_HDO
    51155098
    … …  
    51765159     &    nloc,ncum,nd,i)
    51775160
    5178   USE infotrac_phy, ONLY: niso, ntraciso
    51795161  USE isotopes_mod, ONLY: iso_eau, iso_HDO
    51805162
    … …  
    52525234     &    nloc,ncum,nd,i)
    52535235
    5254   USE infotrac_phy, ONLY: niso,ntraciso
    52555236  USE isotopes_mod, ONLY: iso_eau, iso_HDO
    52565237
    … …  
    53265307     &    nloc,ncum,nd,i)
    53275308
    5328   USE infotrac_phy, ONLY: niso,ntraciso
    53295309  USE isotopes_mod, ONLY: iso_eau, iso_HDO
    53305310
    … …  
    53965376     &    nloc,ncum,nd,i)
    53975377
    5398   USE infotrac_phy, ONLY: niso,ntraciso
    53995378  USE isotopes_mod, ONLY: iso_eau, iso_HDO,ridicule
    54005379
    … …  
    55665545     &    nloc,ncum,nd,i,frac_sublim)
    55675546
    5568   USE infotrac_phy, ONLY: niso,ntraciso
    55695547  USE isotopes_mod, ONLY: iso_eau, iso_HDO,ridicule
    55705548
    … …  
    57035681     &       zxtrfln_cas,zxt_cas,zxtrfl,zxtrfln,zxt,klon)
    57045682
    5705   USE infotrac_phy, ONLY: niso,ntraciso
    57065683  USE isotopes_mod, ONLY: iso_eau, iso_HDO
    57075684
    … …  
    57395716     &       delP,paprs,k,klon,klev)
    57405717
    5741   USE infotrac_phy, ONLY: niso
    57425718  USE isotopes_mod, ONLY: iso_eau, iso_HDO
    57435719         implicit none
    … …  
    57775753     &       delP,paprs,k,klon,klev)
    57785754
    5779   USE infotrac_phy, ONLY: niso,ntraciso
    57805755  USE isotopes_mod, ONLY: iso_eau, iso_HDO
    57815756         implicit none
    … …  
    58285803     &       delP,paprs,k,klon,klev,frac_sublim)
    58295804
    5830   USE infotrac_phy, ONLY: niso,ntraciso
    58315805  USE isotopes_mod, ONLY: iso_eau, iso_HDO
    58325806#ifdef ISOVERIF
    … …  
    59055879     &          qp0,A,m0,beta,gama,g0) 
    59065880
    5907   USE infotrac_phy, ONLY: niso
    59085881  USE isotopes_mod, ONLY: iso_eau, iso_HDO,ntot
    59095882#ifdef ISOVERIF
    … …  
    61006073 
    61016074
    6102   USE infotrac_phy, ONLY: ntraciso,niso,ntraceurs_zone, &
    6103 &       index_trac
    61046075  USE isotopes_mod, ONLY: iso_eau, iso_HDO,thumxt1, &
    61056076&       bidouille_anti_divergence,ridicule
    … …  
    63666337              else !if (qp(il,i).gt.0) then
    63676338                  ! si qp est négatif, on met les isos dedans à 0
    6368                 write(*,*) 'appel_stewart temporaire 230: qp=', &
    6369      &        qp(il,i)
    63706339                do ixt=1,ntraciso
    63716340                 xtp_avantevap(ixt,il)=0.0
    … …  
    76827651     &          )   
    76837652
    7684   USE infotrac_phy, ONLY: niso,ntraciso
    76857653  USE isotopes_mod, ONLY: iso_eau, iso_HDO, Rdefault, ridicule
    76867654#ifdef ISOVERIF
    … …  
    80508018     &  )
    80518019
    8052   USE infotrac_phy, ONLY: niso,ntraciso
    80538020  USE isotopes_mod, ONLY: iso_eau, iso_HDO, Rdefault,ridicule
    80548021#ifdef ISOVERIF
    … …  
    82558222     &          )
    82568223
    8257   USE infotrac_phy, ONLY: niso,ntraciso
    82588224  USE isotopes_mod, ONLY: iso_eau, iso_HDO, Rdefault,ridicule
    82598225#ifdef ISOVERIF
    … …  
    83928358     &          ,xtp_cas,xtwater_cas,xtevap_cas)
    83938359
    8394   USE infotrac_phy, ONLY: niso,ntraciso
    83958360  USE isotopes_mod, ONLY: iso_eau, iso_HDO,no_pce, Rdefault,ridicule       
    83968361#ifdef ISOVERIF
    … …  
    89298894     &          ,xtp_cas,xtwater_cas,xtevap_cas)
    89308895
    8931   USE infotrac_phy, ONLY: niso,ntraciso
    89328896  USE isotopes_mod, ONLY: iso_eau, iso_HDO, Rdefault, ridicule
    89338897#ifdef ISOVERIF
    … …  
    93199283 
    93209284
    9321   USE infotrac_phy, ONLY: niso,ntraciso, &
    9322 &       ntraceurs_zone,index_trac
    93239285  USE isotopes_mod, ONLY: iso_eau, iso_HDO,bidouille_anti_divergence, &
    93249286&       thumxt1, ridicule
    … …  
    95989560              else !if (qp(il,i).gt.0) then
    95999561                  ! si qp est négatif, on met les isos dedans à 0
    9600                 write(*,*) 'appel_stewart_np temporaire 230: qp=', &
    9601      &                         qp(il,i)
    96029562                do ixt=1,ntraciso
    96039563                 xtp_avantevap(ixt,il)=0.0
    … …  
    99819941#ifdef ISOVERIF
    99829942      ! vérif de la compression
    9983       write(*,*) 'appel_stewart_np tmp 506: ', &
    9984      &          'après compress_condensation_facftmr'
     9943!      write(*,*) 'appel_stewart_np tmp 506: ', &
     9944!     &          'après compress_condensation_facftmr'
    99859945!      write(*,*) 'sigd_cas(1:3)=',sigd_cas(1:3)
    9986       if (ncas_condensation_facftmr.ge.4) then
    9987           write(*,*) 'cas_condensation_facftmr(4)=', &
    9988      &          cas_condensation_facftmr(4)
    9989       endif
     9946!      if (ncas_condensation_facftmr.ge.4) then
     9947!          write(*,*) 'cas_condensation_facftmr(4)=', &
     9948!     &          cas_condensation_facftmr(4)
     9949!      endif
    99909950      do il=1,ncas_condensation_facftmr
    99919951        call iso_verif_egalite_choix((Pqisup_cas(il)), &
    … …  
    1015610116#ifdef ISOVERIF
    1015710117      ! vérif de la compression
    10158       write(*,*) 'appel_stewart_np tmp 616: ', &
    10159      &          'apres compress condensation_nofacftmr'
    10160       write(*,*) 'iso_routines 10153: sigd_cas(1:3)=', sigd_cas(1:3)
     10118!      write(*,*) 'appel_stewart_np tmp 616: ', &
     10119!     &          'apres compress condensation_nofacftmr'
     10120!      write(*,*) 'iso_routines 10153: sigd_cas(1:3)=', sigd_cas(1:3)
    1016110121      do il=1,ncas_condensation_nofacftmr
    1016210122        call iso_verif_egalite_choix((Pqisup_cas(il)), &
    … …  
    1024410204       enddo !do izone=1,ntraceurs_zone
    1024510205#ifdef ISOVERIF
    10246        write(*,*) 'appel_stewart_np tmp 690: ', &
    10247      &          'fin du cas condensation_nofacftmr'
     10206!       write(*,*) 'appel_stewart_np tmp 690: ', &
     10207!     &          'fin du cas condensation_nofacftmr'
    1024810208            do il=1,ncas_condensation_nofacftmr
    1024910209               call iso_verif_traceur(xtp &
    … …  
    1072710687
    1072810688#ifdef ISOVERIF
    10729       write(*,*) 'appel_stewart_np tmp 898 après compress glace'
     10689!      write(*,*) 'appel_stewart_np tmp 898 apres compress glace'
    1073010690!      write(*,*) 'qp_avantevap_cas(1),qp_avantevap(cas(1))=',
    1073110691!     &   qp_avantevap_cas(1),qp_avantevap(cas_evap_glace(1))   
    10732       write(*,*) 'Pqiinf_stewart(1)=',Pqiinf_stewart(1)
     10692      !write(*,*) 'Pqiinf_stewart(1)=',Pqiinf_stewart(1)
    1073310693      ! vérif de la compression
    1073410694      do il=1,ncas_evap_glace
    … …  
    1102610986     &          )
    1102710987
    11028   USE infotrac_phy, ONLY: niso,ntraciso
    1102910988  USE isotopes_mod, ONLY: iso_eau, iso_HDO,Rdefault,ridicule
    1103010989#ifdef ISOVERIF
    … …  
    1117611135     &          ,xtp_cas,xtwater_cas,xtevap_cas)
    1117711136
    11178   USE infotrac_phy, ONLY: niso,ntraciso
    1117911137  USE isotopes_mod, ONLY: iso_eau, iso_HDO,Rdefault,no_pce,ridicule
    1118011138#ifdef ISOVERIF
    … …  
    1131111269     &       fac_ftmr_cas(1))
    1131211270        else !if (no_pce.eq.1) then
    11313 #ifdef ISOVERIF                   
    11314             write(*,*) 'appel_stewart_np 1957 tmp'
    11315 #endif           
     11271           
    1131611272      call stewart_explicite_vectall(ncas, &
    1131711273     &       qp_avantevap_cas(1),xtp_avantevap_cas(1,1), &
    … …  
    1177611732     &          ,xtp_cas,xtwater_cas,xtevap_cas)
    1177711733
    11778   USE infotrac_phy, ONLY: niso,ntraciso
    1177911734  USE isotopes_mod, ONLY: iso_eau, iso_HDO, Rdefault, ridicule
    1178011735#ifdef ISOVERIF
    … …  
    1185711812     &       T_cas(1))
    1185811813      else !if (frac_sublim.eq.1) then
    11859 #ifdef ISOVERIF
    11860             write(*,*) 'appel_stewart_explicite_np 2269'
    11861             write(*,*) 'Pqiinf_stewart(1)=',Pqiinf_stewart(1)
    11862             write(*,*) 'Pqisup_cas(1)=',Pqisup_cas(1)
    11863             write(*,*) 'Eqi_cas(1)=',Eqi_cas(1)
    11864             write(*,*) 'Eqi_prime_cas(1)=',Eqi_prime_cas(1)
    11865             write(*,*) 'Eqi_stewart(1)=',Eqi_stewart(1)
    11866 #endif         
     11814!#ifdef ISOVERIF
     11815!            write(*,*) 'appel_stewart_explicite_np 2269'
     11816!            write(*,*) 'Pqiinf_stewart(1)=',Pqiinf_stewart(1)
     11817!            write(*,*) 'Pqisup_cas(1)=',Pqisup_cas(1)
     11818!            write(*,*) 'Eqi_cas(1)=',Eqi_cas(1)
     11819!            write(*,*) 'Eqi_prime_cas(1)=',Eqi_prime_cas(1)
     11820!            write(*,*) 'Eqi_stewart(1)=',Eqi_stewart(1)
     11821!#endif         
    1186711822            call stewart_sublim_nofrac_vectall( &
    1186811823     &        ncas,qp_avantevap_cas(1), &
    … …  
    1220412159     &          tcond,zfice,zxtice,zxtliq)
    1220512160
    12206     USE infotrac_phy, ONLY: ntraciso,niso
    1220712161    USE isotopes_mod, ONLY: iso_eau,iso_HDO,essai_convergence, &
    1220812162&       bidouille_anti_divergence,ridicule
    … …  
    1243812392     &          tcond,zfice,zxtice,zxtliq,n)
    1243912393
    12440     USE infotrac_phy, ONLY: ntraciso,niso
    1244112394    USE isotopes_mod, ONLY: iso_eau,iso_HDO,iso_O18,bidouille_anti_divergence, &
    1244212395&       ridicule
    … …  
    1289612849     &          tcond,zfice,zxtice,zxtliq)
    1289712850
    12898     USE infotrac_phy, ONLY: ntraciso
    1289912851    USE isotopes_mod, ONLY: iso_eau,iso_HDO,bidouille_anti_divergence, &
    1290012852&       ridicule,iso_O18
    … …  
    1309413046     &           xt1lay,q1lay,tsurf,t_coup,nisurf,Rland_ice)
    1309513047
    13096     USE infotrac_phy, ONLY: ntraciso,niso
    1309713048    USE isotopes_mod, ONLY: Rdefault,iso_eau,iso_HDO, &
    1309813049&       bidouille_anti_divergence, ridicule,ridicule_snow, &
    … …  
    1366413615     &   )
    1366513616
    13666     USE infotrac_phy, ONLY: ntraciso,niso
    1366713617    USE isotopes_mod, ONLY: iso_eau,iso_HDO,cste_surf_cond, &
    1366813618&       rh_cste_surf_cond,Rdefault,T_cste_surf_cond,iso_O17,iso_O18, &
    … …  
    1398813938     &   )
    1398913939
    13990     USE infotrac_phy, ONLY: ntraciso,niso
    1399113940    USE isotopes_mod, ONLY: tcorr, toce, alpha_liq_sol,ridicule_evap, &
    1399213941        iso_eau,iso_HDO
    … …  
    1405013999
    1405114000#ifdef ISOVERIF
    14052         write(*,*) 'calcul_iso_surf_sic 175: entree'
    14053 #endif 
    14054 #ifdef ISOVERIF
    1405514001        do i=1,knon
    1405614002         do ixt=1,ntraciso
    … …  
    1424714193     &   )
    1424814194
    14249     USE infotrac_phy, ONLY: ntraciso,niso
    1425014195    USE isotopes_mod, ONLY: h_land_ice, ridicule,ridicule_snow,ridicule_evap, &
    1425114196        iso_eau,iso_HDO,iso_O18
    … …  
    1458214527     &   )
    1458314528
    14584 USE infotrac_phy, ONLY: niso,ntraciso
    1458514529USE isotopes_mod, ONLY: tdifrel,tdifexp_sol, iso_eau, iso_HDO, &
    1458614530&       bidouille_anti_divergence,ruissellement_pluie, Rdefault,Kd, &
    … …  
    1471814662
    1471914663#ifdef ISOVERIF   
    14720       write(*,*) 'calcul_iso_surf_ter 494'
     14664!      write(*,*) 'calcul_iso_surf_ter 494'
    1472114665      do i=1,knon         
    1472214666        if (iso_eau.gt.0) then 
    … …  
    1476614710#endif     
    1476714711#ifdef ISOVERIF
    14768         write(*,*) 'calcul_iso_surf_ter 910'
     14712!        write(*,*) 'calcul_iso_surf_ter 910'
    1476914713        do i=1,knon
    1477014714          if (iso_eau.gt.0) then
    … …  
    1601015954      !USE write_field_phy
    1601115955      USE indice_sol_mod, only: nbsrf 
    16012   USE infotrac_phy, ONLY: ntraciso,niso
    1601315956  USE isotopes_mod, ONLY: initialisation_iso, iso_eau,iso_HDO, &
    1601415957        ridicule_qsol,tnat, P_veg,iso_O18,ridicule, ridicule_snow,iso_O17, &
    … …  
    1602515968#include "dimsoil.h"
    1602615969#include "clesphys.h"
    16027 #include "thermcell.h"
    1602815970#include "compbl.h"     
    1602915971
    … …  
    1607216014!      write(*,*) 'xtsnow(4,8,1)=',xtsnow(4,8,1)
    1607316015#ifdef ISOVERIF
    16074       write(*,*) 'phyiso_etat0 15993 tmp: xtsol(iso_eau,1),qsol(1)=',xtsol(iso_eau,1),qsol(1)
    1607516016      do i=1,klon
    1607616017         do ixt=1,niso
    … …  
    1619816139      !USE write_field_phy
    1619916140      USE indice_sol_mod, only: nbsrf
    16200   USE infotrac_phy, ONLY: ntraciso,niso
    1620116141  USE isotopes_mod, ONLY: tnat,iso_HDO,iso_O18,iso_HTO, iso_eau,toce, &
    1620216142&       Rdefault,iso_O17,ridicule,ridicule_qsol
    … …  
    1621116151#include "dimsoil.h"
    1621216152#include "clesphys.h"
    16213 #include "thermcell.h"
     16153!  #include "thermcell.h"
    1621416154#include "compbl.h"
    1621516155
    … …  
    1637716317     &                  *(1-kcin(ixt))/(1.0-kcin(ixt)*h0)
    1637816318           xt_ancien(ixt,i,k)=q_ancien(i,k)*RMerlivat(ixt) &
    16379      &                  *(q_ancien(i,k)/q0)**(alpha(ixt)-1.0)
     16319     &                  *(min(q0,q_ancien(i,k))/q0)**(alpha(ixt)-1.0)
    1638016320           if (q_ancien(i,k).gt.ridicule) then
    1638116321           xtl_ancien(ixt,i,k)=ql_ancien(i,k)*alpha(ixt) &
    … …  
    1647716417
    1647816418#ifdef ISOVERIF
     16419      write(*,*) 'phyisoetat0 16468: verif init dure'
    1647916420      do i=1,klon
    1648016421         do ixt=1,niso
    … …  
    1650816449           if (iso_eau.gt.0) then
    1650916450             call iso_verif_egalite(xt_ancien(iso_eau,i,k), &
    16510      &           q_ancien(i,k),'phyiso_etat0_dur 775')       
     16451     &           q_ancien(i,k),'phyiso_etat0_dur 775a')       
    1651116452           endif !if (iso_eau.gt.0) then
     16453           if (iso_HDO.gt.0) then
     16454             if (q_ancien(i,k).gt.ridicule) then
     16455              call iso_verif_aberrant_encadre( &
     16456     &           xt_ancien(iso_hdo,i,k)/q_ancien(i,k), &
     16457     &          'phyiso_etat0_dur 775b')
     16458             endif !if (q_ancien(i,k).gt.ridicule) then
     16459           endif !if (iso_HDO.gt.0) then
     16460           if ((iso_HDO.gt.0).and.(iso_O18.gt.0)) then
     16461              if (q_ancien(i,k).gt.ridicule) then
     16462                call iso_verif_O18_aberrant( &
     16463     &              xt_ancien(iso_hdo,i,k)/q_ancien(i,k), &
     16464     &              xt_ancien(iso_O18,i,k)/q_ancien(i,k), &
     16465     &              'phyiso_etat0_dur 775c')
     16466              endif ! if (q_ancien(i,k).gt.ridicule) then
     16467           endif ! if ((iso_HDO.gt.0).and.(iso_O18.gt.0)) then
    1651216468         enddo !do k=1,klev
    1651316469         do nsrf=1,nbsrf
    … …  
    1655216508      end subroutine phyiso_etat0_dur
    1655316509
    16554       subroutine phyiso_etat0_fichier( &
    16555      &           snow,run_off_lic_0, &
    16556      &           xtsnow,xtrun_off_lic_0, &
    16557      &           Rland_ice)
    16558       USE dimphy, only: klon,klev
    16559       !USE mod_grid_phy_lmdz
    16560       !USE mod_phys_lmdz_para
    16561       USE iophy
    16562       USE phys_state_var_mod, ONLY: q_ancien,xt_ancien,wake_deltaq,wake_deltaxt, &
    16563 #ifdef ISOVERIF
    16564         rain_fall,snow_fall,fevap,qsol, &
    16565 #endif
    16566         xtrain_fall,xtsnow_fall,ql_ancien,xtl_ancien,qs_ancien,xts_ancien, &
    16567         fxtevap,xtsol
    16568       !USE iostart
    16569       !USE write_field_phy
    16570       USE indice_sol_mod, only: nbsrf 
    16571   USE infotrac_phy, ONLY: ntraciso,niso
    16572   USE isotopes_mod, ONLY: striso,iso_HDO,iso_eau
    16573 #ifdef ISOVERIF
    16574   USE isotopes_verif_mod
     16510SUBROUTINE phyiso_etat0_fichier(snow, run_off_lic_0, xtsnow, xtrun_off_lic_0, Rland_ice)
     16511   USE dimphy,             ONLY: klon,klev
     16512   USE iophy
     16513   USE phys_state_var_mod, ONLY: q_ancien, xt_ancien, wake_deltaq, wake_deltaxt, &
     16514#ifdef ISOVERIF
     16515     rain_fall, snow_fall, fevap,qsol, &
     16516#endif
     16517     xtrain_fall, xtsnow_fall, ql_ancien, xtl_ancien, qs_ancien, xts_ancien, fxtevap, xtsol
     16518   USE indice_sol_mod,    ONLY: nbsrf 
     16519   USE isotopes_mod,      ONLY: isoName,iso_HDO,iso_eau
     16520   USE phyetat0_mod,      ONLY: phyetat0_get, phyetat0_srf
     16521   USE readTracFiles_mod, ONLY: new2oldH2O
     16522   USE strings_mod,       ONLY: strIdx, strHead, strTail, maxlen, msg, int2str
     16523#ifdef ISOVERIF
     16524   USE isotopes_verif_mod
    1657516525#endif
    1657616526#ifdef ISOTRAC
    16577  USE isotrac_mod, ONLY: strtrac,initialisation_isotrac,index_iso, &
    16578 &       index_zone,izone_init
    16579 #endif
    16580         implicit none
     16527   USE isotrac_mod, ONLY: strtrac, initialisation_isotrac, index_iso, index_zone, izone_init
     16528#endif
     16529   IMPLICIT NONE
    1658116530
    1658216531#include "netcdf.inc"
    1658316532#include "dimsoil.h"
    1658416533#include "clesphys.h"
    16585 #include "thermcell.h"
    1658616534#include "compbl.h"   
    1658716535
    16588         ! inputs
    16589         !REAL qsol(klon)
    16590         REAL snow(klon,nbsrf)
    16591         !REAL evap(klon,nbsrf)
    16592         REAL run_off_lic_0(klon)
    16593         ! outputs   
    16594         !REAL xtsol(niso,klon)
    16595         REAL xtsnow(niso,klon,nbsrf)
    16596         !REAL xtevap(ntraciso,klon,nbsrf)     
    16597         REAL xtrun_off_lic_0(niso,klon)
    16598         REAL Rland_ice(niso,klon)
    16599 
    16600         ! locals
    16601         real iso_tmp(klon)
    16602         real iso_tmp_lonlev(klon,klev)
    16603         real iso_tmp_lonsrf(klon,nbsrf)
    16604         INTEGER ierr
    16605         integer i,ixt,k,nsrf
    16606         INTEGER nid, nvarid
    16607         CHARACTER*2 str2
    16608         CHARACTER*5 str5
    16609         real xmin,xmax   
    16610         CHARACTER*50 striso_sortie 
    16611         integer lnblnk
    16612         LOGICAL :: found,phyetat0_get,phyetat0_srf
    16613 
    16614 !#ifdef ISOVERIF
    16615 !      integer iso_verif_egalite_nostop
    16616 !#endif
    16617 !#ifdef ISOVERIF
    16618 !        real deltaD
    16619 !        integer iso_verif_noNaN_nostop
    16620 !#endif
     16536   REAL, INTENT(IN) ::             snow     (klon,nbsrf)
     16537   REAL, INTENT(IN) ::    run_off_lic_0     (klon)
     16538   REAL, INTENT(OUT) ::          xtsnow(niso,klon,nbsrf)
     16539   REAL, INTENT(OUT) :: xtrun_off_lic_0(niso,klon)
     16540   REAL, INTENT(OUT) ::       Rland_ice(niso,klon)
     16541
     16542   INTEGER :: ierr, i, ixt, k, nsrf, nid, nvarid, lnblnk
     16543   CHARACTER(LEN=2) :: str2
     16544   CHARACTER(LEN=5) :: str5
     16545   CHARACTER(LEN=maxlen) :: outiso, oldIso, modname, nam(2)
     16546   REAL :: xmin, xmax
     16547   LOGICAL :: found
    1662116548#ifdef ISOTRAC
    16622         integer iiso,izone
    16623 #endif
    16624 
    16625 
    16626    write(*,*) 'phyiso_etat0_fichier 3'
    16627    write(*,*) 'niso=',niso
    16628    write(*,*) 'striso(1)=',striso(1)
    16629 
    16630    do ixt=1,ntraciso
    16631 
    16632      if (ixt.le.niso) then
    16633         striso_sortie=striso(ixt)
    16634      else
     16549   INTEGER :: iiso, izone
     16550#endif
     16551
     16552   modname = 'phyiso_etat0_fichier'
     16553   CALL msg('3', modname)
     16554   CALL msg('niso = '//TRIM(int2str(niso)), modname)
     16555   CALL msg('isoName(1) = '//TRIM(isoName(1)), modname)
     16556
     16557   DO ixt = 1, ntraciso
     16558
     16559      outiso = isoName(ixt)
     16560      oldIso = strTail(new2oldH2O(outiso), '_', lFirst=.TRUE.)
     16561      ! on lit seulement si ixt<=niso ou si on initialise les traceurs d'après fichier:
    1663516562#ifdef ISOTRAC
    16636         iiso=index_iso(ixt)
    16637         izone=index_zone(ixt)       
    16638         striso_sortie=striso(iiso)//strtrac(izone)
    16639 #else
    16640         write(*,*) 'phyredem 546: ixt,ntraciso=', ixt,ntraciso
    16641         stop
    16642 #endif
    16643      endif !if (ixt.le.niso) then
    16644      write(*,*) 'phyiso_etat0_fichier 16621: ixt,striso_sortie=',ixt,striso_sortie(1:lnblnk(striso_sortie))
    16645 
    16646            
    16647       ! on lit seulement si ixt<=niso ou si on initialise les traceurs d'après
    16648       ! fichier:
     16563      IF(ixt <= niso .OR. initialisation_isotrac == 0) THEN
     16564#endif
     16565      found = phyetat0iso_srf3(xtsnow,      "XTSNOW", "Surface snow", 0.)
     16566      if (.NOT.found) CALL abort_physic('isotopes_routines_mod', 'phyiso_etat0_fichier 16581: unfound isotopic variable',1)
     16567      found = phyetat0iso_srf3(fxtevap,     "XTEVAP", "evaporation",  0.)
     16568      found = phyetat0iso_get2(xtrain_fall, "xtrain_f", "xrain fall", 0.)
     16569      found = phyetat0iso_get2(xtrain_fall, "xtsnow_f", "xsnow fall", 0.)
     16570      found = phyetat0iso_get3(xt_ancien,   "XTANCIEN",  "QANCIEN",   0.)
     16571      found = phyetat0iso_get3(xtl_ancien,  "XTLANCIEN", "QLANCIEN",  0.)
     16572      found = phyetat0iso_get3(xts_ancien,  "XTASNCIEN", "QSANCIEN",  0.)
     16573      found = phyetat0iso_get2(xtrun_off_lic_0, "XTRUNOFFLIC0", "RUNOFFLIC0", 0.)
     16574      found = phyetat0iso_get3(wake_deltaxt,  "WAKE_DELTAXT", "Delta hum. wake/env",  0.)
     16575#ifdef ISOVERIF
     16576      IF(ixt == iso_eau .AND. iso_eau > 0) THEN
     16577         DO i=1,klon
     16578            CALL iso_verif_egalite(xtrain_fall(iso_eau,i),rain_fall(i),TRIM(modname)//' 231a')
     16579            CALL iso_verif_egalite(xtsnow_fall(iso_eau,i),snow_fall(i),TRIM(modname)//' 231b')
     16580            DO nsrf = 1, nbsrf
     16581               CALL iso_verif_egalite(fxtevap(iso_eau,i,nsrf),fevap(i,nsrf),TRIM(modname)//' 231c')
     16582               CALL iso_verif_egalite( xtsnow(iso_eau,i,nsrf), snow(i,nsrf),TRIM(modname)//' 231d')
     16583            END DO
     16584         END DO
     16585      END IF
     16586      IF(ixt == iso_HDO .AND. iso_HDO > 0) THEN
     16587         DO k=1,klev
     16588            DO i=1,klon
     16589               IF(q_ancien(i,k) > 2e-3) &
     16590                  CALL iso_verif_aberrant(xt_ancien(iso_hdo,i,k)/q_ancien(i,k),TRIM(modname)//' 312')
     16591            END DO
     16592         END DO
     16593      END IF
     16594      IF(iso_eau > 0 .AND. ixt == iso_eau) THEN
     16595         DO i=1,klon
     16596            IF(iso_verif_egalite_nostop(run_off_lic_0(i),xtrun_off_lic_0(iso_eau,i),TRIM(modname)//' 326') == 1) THEN
     16597               WRITE(*,*) 'i=',i
     16598               STOP
     16599            END IF
     16600         END DO
     16601      END IF
     16602#endif
     16603      ! ces variables n'ont pas de traceurs:
     16604      IF(ixt <= niso) THEN
     16605         found = phyetat0iso_get2(xtsol, "XTSOL", "Surface humidity / bucket", 0.)
     16606         found = phyetat0iso_get2(Rland_ice, "Rland_ice", "SR land ice", 0.)
     16607#ifdef ISOVERIF
     16608
     16609         DO i=1,klon
     16610            IF(iso_verif_noNaN_nostop(xtsol(ixt,i),TRIM(modname)//' 95') == 1) THEN
     16611               WRITE(*,*) 'ixt,i=',ixt,i
     16612               STOP
     16613            END IF
     16614         END DO
     16615#endif
     16616      END IF
    1664916617#ifdef ISOTRAC
    16650       if ((ixt.le.niso).or.(initialisation_isotrac.eq.0)) then
    16651 #endif
    16652 
    16653       found=phyetat0_srf(1,iso_tmp_lonsrf,"XTSNOW"//striso_sortie(1:lnblnk(striso_sortie)), &
    16654      &     "Surface snow",0.)
    16655       if (.NOT.found) then
    16656         CALL abort_physic('isotopes_routines_mod', &
    16657                 'phyiso_etat0_fichier 16581: variable isotopique not found',1)
    16658       endif
    16659       xtsnow(ixt,:,:)=iso_tmp_lonsrf(:,:)
    16660      
    16661       found=phyetat0_srf(1,iso_tmp_lonsrf,"XTEVAP"//striso_sortie &
    16662      &   (1:lnblnk(striso_sortie)),"evaporation",0.)
    16663       fxtevap(ixt,:,:)=iso_tmp_lonsrf(:,:)
    16664 
    16665       found=phyetat0_get(1,iso_tmp,"xtrain_f"//striso_sortie &
    16666      &   (1:lnblnk(striso_sortie)),"xrain fall",0.)
    16667       xtrain_fall(ixt,:)=iso_tmp(:)
    16668 
    16669       found=phyetat0_get(1,iso_tmp,"xtsnow_f"//striso_sortie &
    16670      &   (1:lnblnk(striso_sortie)),"snow fall",0.)
    16671       xtsnow_fall(ixt,:)=iso_tmp(:)
    16672 
    16673       found=phyetat0_get(klev,iso_tmp_lonlev,"XTANCIEN"//striso_sortie &
    16674      &       (1:lnblnk(striso_sortie)),"QANCIEN",0.)
    16675       xt_ancien(ixt,:,:)=iso_tmp_lonlev(:,:)
    16676 
    16677       found=phyetat0_get(klev,iso_tmp_lonlev,"XTLANCIEN"//striso_sortie &
    16678      &       (1:lnblnk(striso_sortie)),"QLANCIEN",0.)
    16679       xtl_ancien(ixt,:,:)=iso_tmp_lonlev(:,:)
    16680 
    16681       found=phyetat0_get(klev,iso_tmp_lonlev,"XTSANCIEN"//striso_sortie &
    16682      &       (1:lnblnk(striso_sortie)),"QSANCIEN",0.)
    16683       xts_ancien(ixt,:,:)=iso_tmp_lonlev(:,:)
    16684 
    16685 
    16686       found=phyetat0_get(1,iso_tmp,"XTRUNOFFLIC0"//striso_sortie(1:lnblnk(striso_sortie)), &
    16687      &          "RUNOFFLIC0",0.) 
    16688       xtrun_off_lic_0(ixt,:)=iso_tmp(:)
    16689 
    16690 
    16691       found=phyetat0_get(klev,iso_tmp_lonlev,"WAKE_DELTAXT"//striso_sortie &
    16692      &   (1:lnblnk(striso_sortie)),"Delta hum. wake/env",0.) 
    16693       wake_deltaxt(ixt,:,:)=iso_tmp_lonlev(:,:)
    16694 
    16695 #ifdef ISOVERIF           
    16696       if ((ixt.eq.iso_eau).and.(iso_eau.gt.0)) then
    16697         do i=1,klon
    16698          call iso_verif_egalite(xtrain_fall(iso_eau,i),rain_fall(i), &
    16699      &           'phyisoetat0_fichier 231a')
    16700          call iso_verif_egalite(xtsnow_fall(iso_eau,i),snow_fall(i), &
    16701      &           'phyisoetat0_fichier 231b')
    16702          DO nsrf = 1, nbsrf
    16703          call iso_verif_egalite(fxtevap(iso_eau,i,nsrf),fevap(i,nsrf), &
    16704      &           'phyisoetat0_fichier 231c')
    16705          call iso_verif_egalite(xtsnow(iso_eau,i,nsrf),snow(i,nsrf), &
    16706      &           'phyisoetat0_fichier 231d')
    16707          enddo !DO nsrf = 1, nbsrf
    16708         enddo !do i=1,klon       
    16709       endif !if (iso_eau.gt.0) then 
    16710         if ((iso_HDO.gt.0).and.(ixt.eq.iso_HDO)) then
    16711               do k=1,klev
    16712                do i=1,klon
    16713                 if (q_ancien(i,k).gt.2e-3) then
    16714                 call iso_verif_aberrant(xt_ancien(iso_hdo,i,k) &
    16715      &           /q_ancien(i,k),'phyisoetat0_fichier 312')
    16716                 endif !if (q_ancien(i,k).gt.2e-3) then
    16717                enddo !do i=1,klon
    16718               enddo !do k=1,klev
    16719       endif !if ((iso_HDO.gt.0).and.(ixt.eq.iso_HDO)) then
    16720       if (iso_eau.gt.0) then   
    16721         do i=1,klon
    16722           if (iso_verif_egalite_nostop(run_off_lic_0(i), &
    16723      &           xtrun_off_lic_0(iso_eau,i), &
    16724      &          'phyiso_etat0_fichier 326').eq.1) then
    16725             write(*,*) 'i=',i
    16726             stop
    16727           endif !if (iso_verif_egalite_nostop(run_off_lic_0(i),
    16728         enddo !do i=1,klon
    16729       endif !if (iso_eau.gt.0) then
    16730 #endif
    16731 
    16732        ! ces variables n'ont pas de traceurs:
    16733        if (ixt.le.niso) then
    16734         found=phyetat0_get(1,iso_tmp,"XTSOL"//striso_sortie(1:lnblnk(striso_sortie)), &
    16735      &     "Surface hmidity / bucket",0.) 
    16736         xtsol(ixt,:)=iso_tmp(:)
    16737 
    16738         found=phyetat0_get(1,iso_tmp,"Rland_ice"//striso_sortie &
    16739      &     (1:lnblnk(striso_sortie)),"R land ice",0.)
    16740         Rland_ice(ixt,:)=iso_tmp(:)
    16741 
    16742 #ifdef ISOVERIF
    16743       do i=1,klon
    16744           if (iso_verif_noNaN_nostop(xtsol(ixt,i), &
    16745      &          'phyiso_etat0_fichier 95').eq.1) then
    16746             write(*,*) 'ixt,i=',ixt,i
    16747             stop
    16748           endif       
    16749       enddo !do i=1,klon
    16750 #endif
    16751 
    16752        endif
     16618      END IF ! IF(ixt > niso .OR. initialisation_isotrac == 0))
     16619#endif
     16620
     16621   END DO
    1675316622
    1675416623#ifdef ISOTRAC
    16755      endif !if ((ixt.le.niso).or.(initialisation_isotrac.eq.0)) then
    16756 #endif
    16757 
    16758   enddo !do ixt=1,ntraciso
    16759 
    16760 #ifdef ISOTRAC
    16761         if (initialisation_isotrac.ne.0) then
    16762         ! on n'initialise pas d'après le fichier
    16763         ! l'eau normale est mise dans la zone izone_init
    16764 
    16765         do ixt=niso+1,ntraciso
    16766 
    16767              iiso=index_iso(ixt)
    16768 
    16769              if (index_zone(ixt).eq.izone_init) then
    16770                 do i=1,klon
    16771                  do nsrf = 1, nbsrf
    16772                   fxtevap(ixt,i,nsrf)=fxtevap(iiso,i,nsrf)
    16773                  enddo !do nsrf = 1, nbsrf
    16774                  xtsnow_fall(ixt,i)=xtsnow_fall(iiso,i)
    16775                  xtrain_fall(ixt,i)=xtrain_fall(iiso,i)
    16776                  do k=1,klev
    16777                    xt_ancien(ixt,i,k)=xt_ancien(iiso,i,k)
    16778                    xtl_ancien(ixt,i,k)=xtl_ancien(iiso,i,k)
    16779                    xts_ancien(ixt,i,k)=xts_ancien(iiso,i,k)
    16780                    wake_deltaxt(ixt,i,k)= wake_deltaxt(iiso,i,k)   
    16781                  enddo
    16782                 enddo !do i=1,klon
    16783              else !if (index_zone(ixt).eq.izone_init) then
    16784                 do i=1,klon
    16785                  do nsrf = 1, nbsrf
    16786                   fxtevap(ixt,i,nsrf)=0.0
    16787                  enddo !do nsrf = 1, nbsrf
    16788                  xtsnow_fall(ixt,i)=0.0
    16789                  xtrain_fall(ixt,i)=0.0
    16790                  do k=1,klev
    16791                    xt_ancien(ixt,i,k)=0.0
    16792                    xtl_ancien(ixt,i,k)=0.0
    16793                    xts_ancien(ixt,i,k)=0.0
    16794                  enddo
    16795                 enddo !do i=1,klon
    16796              endif !if (index_zone(ixt).eq.izone_init) then
    16797 
    16798          enddo  !do ixt=1,niso
    16799       endif !if (initialisation_isotrac.eq.0) then
    16800 
    16801 
    16802 #ifdef ISOVERIF
    16803         DO nsrf = 1, nbsrf
    16804          do i=1,klon
    16805                call iso_verif_traceur(fxtevap(1,i,nsrf), &
    16806      &                   'phyiso_etat0_fichier 426')
    16807          enddo !do i=1,klon
    16808         enddo !DO nsrf = 1, nbsrf
    16809         do i=1,klon
    16810            call iso_verif_traceur(xtrain_fall(1,i), &
    16811      &                   'phyiso_etat0_fichier 466')
    16812            call iso_verif_traceur(xtsnow_fall(1,i), &
    16813      &                   'phyiso_etat0_fichier 468')
    16814         enddo !do i=1,klon
    16815         do k=1,klev
    16816           do i=1,klon
    16817                call iso_verif_traceur(xt_ancien(1,i,k), &
    16818      &                   'phyiso_etat0_fichier 591')
    16819           enddo !do i=1,klon
    16820         enddo !do k=1,klev             
     16624   IF(initialisation_isotrac /= 0) THEN
     16625      ! On n'initialise pas d'apres le fichier. L'eau normale est mise dans la zone izone_init
     16626      DO ixt=niso+1,ntraciso
     16627         iiso=index_iso(ixt)
     16628         IF(index_zone(ixt) == izone_init) THEN
     16629            DO i = 1, klon
     16630               fxtevap(ixt,i,1:nsrf) = fxtevap(iiso,i,1:nsrf)
     16631               xtsnow_fall(ixt,i) = xtsnow_fall(iiso,i)
     16632               xtrain_fall(ixt,i) = xtrain_fall(iiso,i)
     16633               DO k = 1, klev
     16634                  xt_ancien   (ixt,i,k) = xt_ancien   (iiso,i,k)
     16635                  xtl_ancien  (ixt,i,k) = xtl_ancien  (iiso,i,k)
     16636                  xts_ancien  (ixt,i,k) = xts_ancien  (iiso,i,k)
     16637                  wake_deltaxt(ixt,i,k) = wake_deltaxt(iiso,i,k)   
     16638               END DO
     16639            END DO
     16640         ELSE
     16641            DO i = 1, klon
     16642               fxtevap(ixt,i,1:nbsrf)=0.0
     16643               xtsnow_fall(ixt,i)=0.0
     16644               xtrain_fall(ixt,i)=0.0
     16645               xt_ancien (ixt,i,1:klev) = 0.0
     16646               xtl_ancien(ixt,i,1:klev) = 0.0
     16647               xts_ancien(ixt,i,1:klev) = 0.0
     16648            END DO
     16649         END IF
     16650      END DO
     16651   END IF
     16652
     16653#ifdef ISOVERIF
     16654   DO nsrf = 1, nbsrf
     16655      DO i = 1, klon
     16656         CALL iso_verif_traceur(fxtevap(1,i,nsrf), 'phyiso_etat0_fichier 426')
     16657      END DO
     16658   END DO
     16659   DO i=1,klon
     16660      CALL iso_verif_traceur(xtrain_fall(1,i), 'phyiso_etat0_fichier 466')
     16661      CALL iso_verif_traceur(xtsnow_fall(1,i), 'phyiso_etat0_fichier 468')
     16662   END DO
     16663   DO k = 1, klev
     16664      DO i = 1, klon
     16665         CALL iso_verif_traceur(xt_ancien(1,i,k), 'phyiso_etat0_fichier 591')
     16666      END DO
     16667   END DO
    1682116668#endif
    1682216669        ! endif ISOVERIF       
    … …  
    1682416671        ! endif ISOTRAC     
    1682516672
    16826 ! on ferme le fichier
    16827 !      CALL close_startphy
    16828 ! déjà fermé dans phyetat0
     16673CONTAINS
     16674
     16675LOGICAL FUNCTION phyetat0iso_get2(field, pref, descr, default) RESULT(lFound)
     16676  REAL,             INTENT(INOUT) :: field(:,:)
     16677  CHARACTER(LEN=*), INTENT(IN)    :: pref, descr
     16678  REAL,             INTENT(IN)    :: default
     16679  REAL :: iso_tmp(klon)
     16680  nam(1) = TRIM(pref)//TRIM(outiso)
     16681  nam(2) = TRIM(pref)//TRIM(oldIso)
     16682  lFound = phyetat0_get(iso_tmp, nam, descr, default)
     16683  field(ixt,:) = iso_tmp
     16684END FUNCTION phyetat0iso_get2
     16685
     16686
     16687LOGICAL FUNCTION phyetat0iso_get3(field, pref, descr, default) RESULT(lFound)
     16688  REAL,             INTENT(INOUT) :: field(:,:,:)
     16689  CHARACTER(LEN=*), INTENT(IN)    :: pref, descr
     16690  REAL,             INTENT(IN)    :: default
     16691  REAL :: iso_tmp_lonlev(klon,klev)
     16692  nam(1) = TRIM(pref)//TRIM(outiso)
     16693  nam(2) = TRIM(pref)//TRIM(oldIso)
     16694  lFound = phyetat0_get(iso_tmp_lonlev, nam, descr, default)
     16695  field(ixt,:,:) = iso_tmp_lonlev(:,:)
     16696END FUNCTION phyetat0iso_get3
     16697
     16698LOGICAL FUNCTION phyetat0iso_srf3(field, pref, descr, default) RESULT(lFound)
     16699  REAL,             INTENT(INOUT) :: field(:,:,:)
     16700  CHARACTER(LEN=*), INTENT(IN)    :: pref, descr
     16701  REAL,             INTENT(IN)    :: default
     16702  REAL :: iso_tmp_lonsrf(klon,nbsrf)
     16703  nam(1) = TRIM(pref)//TRIM(outiso)
     16704  nam(2) = TRIM(pref)//TRIM(oldIso)
     16705  lFound = phyetat0_srf(iso_tmp_lonsrf, nam, descr, default)
     16706  field(ixt,:,:) = iso_tmp_lonsrf
     16707END FUNCTION phyetat0iso_srf3
    1682916708
    1683016709        end subroutine phyiso_etat0_fichier
     16710
     16711
    1683116712
    1683216713
    … …  
    1684416725     &           d_xt_decroiss, &
    1684516726     &           xt_seri)
    16846         USE infotrac_phy, only: ntraciso
    1684716727        USE isotopes_mod, only: iso_HTO,ok_prod_nucl_tritium
    1684816728        USE dimphy, only: klon,klev
    … …  
    1836618246!     &                         prod_nucl_HTO)
    1836718247
    18368         USE infotrac_phy, only: ntraciso
    1836918248        use isotopes_mod, only: nessai, lat_nucl, lon_nucl, &
    1837018249&               zmin_nucl, zmax_nucl, HTO_nucl
    … …  
    1858818467     &                                paprs, &
    1858918468     &                                prod_nucl)
    18590         USE infotrac_phy, only: ntraciso
    1859118469        USE isotopes_mod, ONLY: iso_HTO
    1859218470        use geometry_mod, only: cell_area
    … …  
    1873418612     &           tcond,zfice,zxtice,zxtliq)
    1873518613
    18736     USE infotrac_phy, ONLY: ntraciso,niso,index_trac,ntraceurs_zone
    1873718614    USE isotopes_mod, ONLY: iso_eau,iso_HDO,essai_convergence, &
    1873818615&       bidouille_anti_divergence,ridicule
    … …  
    1886418741     &           tcond,zfice,zxtice,zxtliq,n)
    1886518742
    18866     USE infotrac_phy, ONLY: ntraciso,niso,index_trac,ntraceurs_zone
    1886718743    USE isotopes_mod, ONLY: iso_eau,iso_HDO,iso_O18,bidouille_anti_divergence, &
    1886818744&       ridicule
Note: See TracChangeset for help on using the changeset viewer.