Ignore:
Timestamp:
Jun 11, 2008, 1:14:41 PM (16 years ago)
Author:
Laurent Fairhead
Message:

Pour des raisons pratiques (besoin de tester facilement les parametrisations
physiques en 1D et 3D), les arguments des routines de lecture et d'ecriture
de l'etat initial de la physique ont disparu des appels à ces routines et sont
maintenant passés par le module phys_state_var_mod
LF

File:
1 edited

Legend:

Unmodified
Added
Removed

  • LMDZ4/trunk/libf/phylmd/phyredem.F

    r937 r967  
    33!
    44c
    5       SUBROUTINE phyredem (fichnom,dtime,radpas,ocean,
    6      .           rlat_p,rlon_p, pctsrf_p,tsol_p,
    7      .           alb1_p, alb2_p,
    8      .           rain_fall_p, snow_fall_p,solsw_p, sollw_p,
    9      .           radsol_p,zmea_p,zstd_p,zsig_p,
    10      .           zgam_p,zthe_p,zpic_p,zval_p,rugsrel_p,
    11      .           t_ancien_p, q_ancien_p, rnebcon_p, ratqs_p, clwcon_p,
    12      .           pbl_tke_p, zmax0_p, f0_p, ema_work1_p, ema_work2_p)
     5      SUBROUTINE phyredem (fichnom)
    136
    147      USE dimphy
    … …  
    1811      USE fonte_neige_mod,  ONLY : fonte_neige_final
    1912      USE pbl_surface_mod,  ONLY : pbl_surface_final
     13      USE surface_data,     ONLY : ocean, ok_veget
     14      USE phys_state_var_mod
    2015
    2116      IMPLICIT none
    … …  
    3429c======================================================================
    3530      CHARACTER*(*) fichnom
    36       REAL dtime
    37       INTEGER radpas
    38       REAL rlat_p(klon), rlon_p(klon)
    39       REAL tsol_p(klon,nbsrf)
    40       REAL pbl_tke_p(klon,klev,nbsrf)
    41       REAL zmax0_p(klon), f0_p(klon)
    42       REAL ema_work1_p(klon, klev), ema_work2_p(klon, klev)
     31
     32c les variables globales ecrites dans le fichier restart
     33      REAL rlat_glo(klon_glo), rlon_glo(klon_glo)
     34      REAL pctsrf_glo(klon_glo, nbsrf)
     35      REAL tsol_glo(klon_glo,nbsrf)
     36      REAL alb1_glo(klon_glo,nbsrf)
     37      REAL alb2_glo(klon_glo,nbsrf)
     38      REAL rain_fall_glo(klon_glo)
     39      REAL snow_fall_glo(klon_glo)
     40      real solsw_glo(klon_glo)
     41      real sollw_glo(klon_glo)
     42      REAL radsol_glo(klon_glo)
     43      REAL zmea_glo(klon_glo)
     44      REAL zstd_glo(klon_glo)
     45      REAL zsig_glo(klon_glo)
     46      REAL zgam_glo(klon_glo)
     47      REAL zthe_glo(klon_glo)
     48      REAL zpic_glo(klon_glo)
     49      REAL zval_glo(klon_glo)
     50      REAL rugsrel_glo(klon_glo)
     51      REAL t_ancien_glo(klon_glo,klev), q_ancien_glo(klon_glo,klev)
     52      REAL clwcon_glo(klon_glo,klev)
     53      REAL rnebcon_glo(klon_glo,klev)
     54      REAL ratqs_glo(klon_glo,klev)
     55      REAL pbl_tke_glo(klon_glo,klev,nbsrf)
     56      REAL zmax0_glo(klon_glo), f0_glo(klon)
     57      REAL ema_work1_glo(klon_glo, klev), ema_work2_glo(klon_glo, klev)
     58
     59cIM "slab" ocean
    4360      REAL tsoil_p(klon,nsoilmx,nbsrf)
    44       CHARACTER*6 ocean
    45 cIM "slab" ocean
    4661      REAL tslab_p(klon), seaice_p(klon)
    4762      REAL qsurf_p(klon,nbsrf)
    4863      REAL qsol_p(klon)
    4964      REAL snow_p(klon,nbsrf)
    50       REAL alb1_p(klon,nbsrf)   ! albedo in visible SW interval 
    51       REAL alb2_p(klon,nbsrf)   ! albedo in near IR interval
    5265      REAL evap_p(klon,nbsrf)
    53       REAL rain_fall_p(klon)
    54       REAL snow_fall_p(klon)
    55       real solsw_p(klon)
    56       real sollw_p(klon)
    5766      real fder_p(klon)
    58       REAL radsol_p(klon)
    5967      REAL frugs_p(klon,nbsrf)
    6068      REAL agesno_p(klon,nbsrf)
    61       REAL zmea_p(klon)
    62       REAL zstd_p(klon)
    63       REAL zsig_p(klon)
    64       REAL zgam_p(klon)
    65       REAL zthe_p(klon)
    66       REAL zpic_p(klon)
    67       REAL zval_p(klon)
    68       REAL rugsrel_p(klon)
    69       REAL pctsrf_p(klon, nbsrf)
    70       REAL t_ancien_p(klon,klev), q_ancien_p(klon,klev)
    71       real clwcon_p(klon,klev),rnebcon_p(klon,klev),ratqs_p(klon,klev)
    7269      REAL run_off_lic_0_p(klon)
    7370     
    74       REAL rlat(klon_glo), rlon(klon_glo)
    75       REAL tsol(klon_glo,nbsrf)
    76       REAL pbl_tke(klon_glo,klev,nbsrf)
    77       REAL zmax0(klon_glo), f0(klon)
    78       REAL ema_work1(klon, klev), ema_work2(klon, klev)
    7971      REAL tsoil(klon_glo,nsoilmx,nbsrf)
    8072      REAL tslab(klon_glo), seaice(klon_glo)
    … …  
    8274      REAL qsol(klon_glo)
    8375      REAL snow(klon_glo,nbsrf)
    84       REAL alb1(klon_glo,nbsrf)
    85       REAL alb2(klon_glo,nbsrf)
    8676      REAL evap(klon_glo,nbsrf)
    87       REAL rain_fall(klon_glo)
    88       REAL snow_fall(klon_glo)
    89       real solsw(klon_glo)
    90       real sollw(klon_glo)
    9177      real fder(klon_glo)
    92       REAL radsol(klon_glo)
    9378      REAL frugs(klon_glo,nbsrf)
    9479      REAL agesno(klon_glo,nbsrf)
    95       REAL zmea(klon_glo)
    96       REAL zstd(klon_glo)
    97       REAL zsig(klon_glo)
    98       REAL zgam(klon_glo)
    99       REAL zthe(klon_glo)
    100       REAL zpic(klon_glo)
    101       REAL zval(klon_glo)
    102       REAL rugsrel(klon_glo)
    103       REAL pctsrf(klon_glo, nbsrf)
    104       REAL t_ancien(klon_glo,klev), q_ancien(klon_glo,klev)
    105       REAL clwcon(klon_glo,klev)
    106       REAL rnebcon(klon_glo,klev)
    107       REAL ratqs(klon_glo,klev)
    10880      REAL run_off_lic_0(klon_glo)
    10981      REAL masq(klon_glo)
    … …  
    11688c
    11789      INTEGER isoil, nsrf
    118       CHARACTER*7 str7
    119       CHARACTER*2 str2
     90      CHARACTER (len=7) :: str7
     91      CHARACTER (len=2) :: str2
    12092
    12193c======================================================================
    … …  
    139111c======================================================================
    140112
    141       call Gather( rlat_p,rlat)
    142       call Gather( rlon_p,rlon)
    143       call Gather( tsol_p,tsol)
    144       call Gather( pbl_tke_p,pbl_tke)
    145       call Gather( zmax0_p,zmax0)
    146       call Gather( f0_p,f0)
    147       call Gather( ema_work1_p, ema_work1)
    148       call Gather( ema_work2_p, ema_work2)
     113      call Gather( rlat,rlat_glo)
     114      call Gather( rlon,rlon_glo)
     115      call Gather( pctsrf,pctsrf_glo)
     116      call Gather( ftsol,tsol_glo)
     117      call Gather( falb1,alb1_glo)
     118      call Gather( falb2,alb2_glo)
     119      call Gather( rain_fall,rain_fall_glo)
     120      call Gather( snow_fall,snow_fall_glo)
     121      call Gather( sollw,sollw_glo)
     122      call Gather( solsw,solsw_glo)
     123      call Gather( radsol,radsol_glo)
     124      call Gather( zmea,zmea_glo)
     125      call Gather( zstd,zstd_glo)
     126      call Gather( zsig,zsig_glo)
     127      call Gather( zgam,zgam_glo)
     128      call Gather( zthe,zthe_glo)
     129      call Gather( zpic,zpic_glo)
     130      call Gather( zval,zval_glo)
     131      call Gather( rugoro,rugsrel_glo)
     132      call Gather( t_ancien,t_ancien_glo)
     133      call Gather( q_ancien,q_ancien_glo)
     134      call Gather( clwcon,clwcon_glo)
     135      call Gather( rnebcon,rnebcon_glo)
     136      call Gather( ratqs,ratqs_glo)
     137      call Gather( pbl_tke,pbl_tke_glo)
     138      call Gather( zmax0,zmax0_glo)
     139      call Gather( f0,f0_glo)
     140      call Gather( ema_work1, ema_work1_glo)
     141      call Gather( ema_work2, ema_work2_glo)
     142
    149143      call Gather( tsoil_p,tsoil)
    150144      call Gather( tslab_p,tslab)
    … …  
    153147      call Gather( qsol_p,qsol)
    154148      call Gather( snow_p,snow)
    155       call Gather( alb1_p,alb1)
    156       call Gather( alb2_p,alb2)
    157149      call Gather( evap_p,evap)
    158       call Gather( radsol_p,radsol)
    159       call Gather( rain_fall_p,rain_fall)
    160       call Gather( snow_fall_p,snow_fall)
    161       call Gather( sollw_p,sollw)
    162       call Gather( solsw_p,solsw)
    163150      call Gather( fder_p,fder)
    164151      call Gather( frugs_p,frugs)
    165152      call Gather( agesno_p,agesno)
    166       call Gather( zmea_p,zmea)
    167       call Gather( zstd_p,zstd)
    168       call Gather( zsig_p,zsig)
    169       call Gather( zgam_p,zgam)
    170       call Gather( zthe_p,zthe)
    171       call Gather( zpic_p,zpic)
    172       call Gather( zval_p,zval)
    173       call Gather( rugsrel_p,rugsrel)
    174       call Gather( pctsrf_p,pctsrf)
    175153      call Gather( run_off_lic_0_p,run_off_lic_0)
    176       call Gather( t_ancien_p,t_ancien)
    177       call Gather( q_ancien_p,q_ancien)
    178       call Gather( rnebcon_p,rnebcon)
    179       call Gather( clwcon_p,clwcon)
    180       call Gather( ratqs_p,ratqs)
    181154      call Gather( zmasq,masq)
    182155     
    … …  
    245218      ierr = NF_ENDDEF(nid)
    246219#ifdef NC_DOUBLE
    247       ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,rlon)
    248 #else
    249       ierr = NF_PUT_VAR_REAL (nid,nvarid,rlon)
     220      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,rlon_glo)
     221#else
     222      ierr = NF_PUT_VAR_REAL (nid,nvarid,rlon_glo)
    250223#endif
    251224c
    … …  
    260233      ierr = NF_ENDDEF(nid)
    261234#ifdef NC_DOUBLE
    262       ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,rlat)
    263 #else
    264       ierr = NF_PUT_VAR_REAL (nid,nvarid,rlat)
     235      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,rlat_glo)
     236#else
     237      ierr = NF_PUT_VAR_REAL (nid,nvarid,rlat_glo)
    265238#endif
    266239c
    … …  
    295268      ierr = NF_ENDDEF(nid)
    296269#ifdef NC_DOUBLE
    297       ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,pctsrf(1 : klon_glo, is_ter))
    298 #else
    299       ierr = NF_PUT_VAR_REAL (nid,nvarid,pctsrf(1 : klon_glo, is_ter))
     270      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,pctsrf_glo(1 : klon_glo,     &
     271     &   is_ter))
     272#else
     273      ierr = NF_PUT_VAR_REAL (nid,nvarid,pctsrf_glo(1 : klon_glo,       &
     274     &    is_ter))
    300275#endif
    301276C
    … …  
    312287      ierr = NF_ENDDEF(nid)
    313288#ifdef NC_DOUBLE
    314       ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,pctsrf(1 : klon_glo,is_lic))
    315 #else
    316       ierr = NF_PUT_VAR_REAL (nid,nvarid,pctsrf(1 : klon_glo, is_lic))
     289      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,pctsrf_glo(1 : klon_glo,
     290     &                          is_lic))
     291#else
     292      ierr = NF_PUT_VAR_REAL (nid,nvarid,pctsrf_glo(1 : klon_glo,
     293     &                        is_lic))
    317294#endif
    318295C
    … …  
    329306      ierr = NF_ENDDEF(nid)
    330307#ifdef NC_DOUBLE
    331       ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,pctsrf(1 : klon_glo, is_oce))
    332 #else
    333       ierr = NF_PUT_VAR_REAL (nid,nvarid,pctsrf(1 : klon_glo, is_oce))
     308      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,pctsrf_glo(1 : klon_glo,
     309     &                          is_oce))
     310#else
     311      ierr = NF_PUT_VAR_REAL (nid,nvarid,pctsrf_glo(1 : klon_glo,
     312     &                        is_oce))
    334313#endif
    335314C
    … …  
    346325      ierr = NF_ENDDEF(nid)
    347326#ifdef NC_DOUBLE
    348       ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,pctsrf(1 : klon_glo, is_sic))
    349 #else
    350       ierr = NF_PUT_VAR_REAL (nid,nvarid,pctsrf(1 : klon_glo, is_sic))
     327      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,pctsrf_glo(1 : klon_glo,
     328     &                          is_sic))
     329#else
     330      ierr = NF_PUT_VAR_REAL (nid,nvarid,pctsrf_glo(1 : klon_glo,
     331     &                         is_sic))
    351332#endif
    352333C
    … …  
    370351        ENDIF
    371352#ifdef NC_DOUBLE
    372         ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,tsol(1,nsrf))
    373 #else
    374         ierr = NF_PUT_VAR_REAL (nid,nvarid,tsol(1,nsrf))
     353        ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,tsol_glo(1,nsrf))
     354#else
     355        ierr = NF_PUT_VAR_REAL (nid,nvarid,tsol_glo(1,nsrf))
    375356#endif
    376357      ENDDO
    … …  
    487468        ENDIF
    488469#ifdef NC_DOUBLE
    489       ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,alb1(1,nsrf))
    490 #else
    491       ierr = NF_PUT_VAR_REAL (nid,nvarid,alb1(1,nsrf))
     470      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,alb1_glo(1,nsrf))
     471#else
     472      ierr = NF_PUT_VAR_REAL (nid,nvarid,alb1_glo(1,nsrf))
    492473#endif
    493474      ENDDO
    … …  
    510491        ENDIF
    511492#ifdef NC_DOUBLE
    512       ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,alb2(1,nsrf))
    513 #else
    514       ierr = NF_PUT_VAR_REAL (nid,nvarid,alb2(1,nsrf))
     493      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,alb2_glo(1,nsrf))
     494#else
     495      ierr = NF_PUT_VAR_REAL (nid,nvarid,alb2_glo(1,nsrf))
    515496#endif
    516497      ENDDO
    … …  
    575556      ierr = NF_ENDDEF(nid)
    576557#ifdef NC_DOUBLE
    577       ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,radsol)
    578 #else
    579       ierr = NF_PUT_VAR_REAL (nid,nvarid,radsol)
     558      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,radsol_glo)
     559#else
     560      ierr = NF_PUT_VAR_REAL (nid,nvarid,radsol_glo)
    580561#endif
    581562c
    … …  
    590571      ierr = NF_ENDDEF(nid)
    591572#ifdef NC_DOUBLE
    592       ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,solsw)
    593 #else
    594       ierr = NF_PUT_VAR_REAL (nid,nvarid,solsw)
     573      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,solsw_glo)
     574#else
     575      ierr = NF_PUT_VAR_REAL (nid,nvarid,solsw_glo)
    595576#endif
    596577c
    … …  
    605586      ierr = NF_ENDDEF(nid)
    606587#ifdef NC_DOUBLE
    607       ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,sollw)
    608 #else
    609       ierr = NF_PUT_VAR_REAL (nid,nvarid,sollw)
     588      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,sollw_glo)
     589#else
     590      ierr = NF_PUT_VAR_REAL (nid,nvarid,sollw_glo)
    610591#endif
    611592c
    … …  
    635616      ierr = NF_ENDDEF(nid)
    636617#ifdef NC_DOUBLE
    637       ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,rain_fall)
    638 #else
    639       ierr = NF_PUT_VAR_REAL (nid,nvarid,rain_fall)
     618      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,rain_fall_glo)
     619#else
     620      ierr = NF_PUT_VAR_REAL (nid,nvarid,rain_fall_glo)
    640621#endif
    641622c
    … …  
    650631      ierr = NF_ENDDEF(nid)
    651632#ifdef NC_DOUBLE
    652       ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,snow_fall)
    653 #else
    654       ierr = NF_PUT_VAR_REAL (nid,nvarid,snow_fall)
     633      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,snow_fall_glo)
     634#else
     635      ierr = NF_PUT_VAR_REAL (nid,nvarid,snow_fall_glo)
    655636#endif
    656637c
    … …  
    716697      ierr = NF_ENDDEF(nid)
    717698#ifdef NC_DOUBLE
    718       ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,zmea)
    719 #else
    720       ierr = NF_PUT_VAR_REAL (nid,nvarid,zmea)
     699      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,zmea_glo)
     700#else
     701      ierr = NF_PUT_VAR_REAL (nid,nvarid,zmea_glo)
    721702#endif
    722703c
    … …  
    729710      ierr = NF_ENDDEF(nid)
    730711#ifdef NC_DOUBLE
    731       ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,zstd)
    732 #else
    733       ierr = NF_PUT_VAR_REAL (nid,nvarid,zstd)
     712      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,zstd_glo)
     713#else
     714      ierr = NF_PUT_VAR_REAL (nid,nvarid,zstd_glo)
    734715#endif
    735716      ierr = NF_REDEF (nid)
    … …  
    741722      ierr = NF_ENDDEF(nid)
    742723#ifdef NC_DOUBLE
    743       ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,zsig)
    744 #else
    745       ierr = NF_PUT_VAR_REAL (nid,nvarid,zsig)
     724      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,zsig_glo)
     725#else
     726      ierr = NF_PUT_VAR_REAL (nid,nvarid,zsig_glo)
    746727#endif
    747728      ierr = NF_REDEF (nid)
    … …  
    753734      ierr = NF_ENDDEF(nid)
    754735#ifdef NC_DOUBLE
    755       ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,zgam)
    756 #else
    757       ierr = NF_PUT_VAR_REAL (nid,nvarid,zgam)
     736      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,zgam_glo)
     737#else
     738      ierr = NF_PUT_VAR_REAL (nid,nvarid,zgam_glo)
    758739#endif
    759740      ierr = NF_REDEF (nid)
    … …  
    765746      ierr = NF_ENDDEF(nid)
    766747#ifdef NC_DOUBLE
    767       ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,zthe)
    768 #else
    769       ierr = NF_PUT_VAR_REAL (nid,nvarid,zthe)
     748      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,zthe_glo)
     749#else
     750      ierr = NF_PUT_VAR_REAL (nid,nvarid,zthe_glo)
    770751#endif
    771752      ierr = NF_REDEF (nid)
    … …  
    777758      ierr = NF_ENDDEF(nid)
    778759#ifdef NC_DOUBLE
    779       ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,zpic)
    780 #else
    781       ierr = NF_PUT_VAR_REAL (nid,nvarid,zpic)
     760      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,zpic_glo)
     761#else
     762      ierr = NF_PUT_VAR_REAL (nid,nvarid,zpic_glo)
    782763#endif
    783764      ierr = NF_REDEF (nid)
    … …  
    789770      ierr = NF_ENDDEF(nid)
    790771#ifdef NC_DOUBLE
    791       ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,zval)
    792 #else
    793       ierr = NF_PUT_VAR_REAL (nid,nvarid,zval)
     772      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,zval_glo)
     773#else
     774      ierr = NF_PUT_VAR_REAL (nid,nvarid,zval_glo)
    794775#endif
    795776      ierr = NF_REDEF (nid)
    … …  
    801782      ierr = NF_ENDDEF(nid)
    802783#ifdef NC_DOUBLE
    803       ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,rugsrel)
    804 #else
    805       ierr = NF_PUT_VAR_REAL (nid,nvarid,rugsrel)
     784      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,rugsrel_glo)
     785#else
     786      ierr = NF_PUT_VAR_REAL (nid,nvarid,rugsrel_glo)
    806787#endif
    807788c
    … …  
    814795      ierr = NF_ENDDEF(nid)
    815796#ifdef NC_DOUBLE
    816       ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,t_ancien)
    817 #else
    818       ierr = NF_PUT_VAR_REAL (nid,nvarid,t_ancien)
     797      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,t_ancien_glo)
     798#else
     799      ierr = NF_PUT_VAR_REAL (nid,nvarid,t_ancien_glo)
    819800#endif
    820801c
    … …  
    827808      ierr = NF_ENDDEF(nid)
    828809#ifdef NC_DOUBLE
    829       ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,q_ancien)
    830 #else
    831       ierr = NF_PUT_VAR_REAL (nid,nvarid,q_ancien)
     810      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,q_ancien_glo)
     811#else
     812      ierr = NF_PUT_VAR_REAL (nid,nvarid,q_ancien_glo)
    832813#endif
    833814c
    … …  
    857838      ierr = NF_ENDDEF(nid)
    858839#ifdef NC_DOUBLE
    859       ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,clwcon)
    860 #else
    861       ierr = NF_PUT_VAR_REAL (nid,nvarid,clwcon)
     840      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,clwcon_glo)
     841#else
     842      ierr = NF_PUT_VAR_REAL (nid,nvarid,clwcon_glo)
    862843#endif
    863844c
    … …  
    872853      ierr = NF_ENDDEF(nid)
    873854#ifdef NC_DOUBLE
    874       ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,rnebcon)
    875 #else
    876       ierr = NF_PUT_VAR_REAL (nid,nvarid,rnebcon)
     855      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,rnebcon_glo)
     856#else
     857      ierr = NF_PUT_VAR_REAL (nid,nvarid,rnebcon_glo)
    877858#endif
    878859c
    … …  
    887868      ierr = NF_ENDDEF(nid)
    888869#ifdef NC_DOUBLE
    889       ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,ratqs)
    890 #else
    891       ierr = NF_PUT_VAR_REAL (nid,nvarid,ratqs)
     870      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,ratqs_glo)
     871#else
     872      ierr = NF_PUT_VAR_REAL (nid,nvarid,ratqs_glo)
    892873#endif
    893874c
    … …  
    932913        ENDIF
    933914#ifdef NC_DOUBLE
    934         ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,pbl_tke(:,:,nsrf))
    935 #else
    936       ierr = NF_PUT_VAR_REAL (nid,nvarid,pbl_tke(:,:,nsrf))
     915        ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,pbl_tke_glo(:,:,nsrf))
     916#else
     917      ierr = NF_PUT_VAR_REAL (nid,nvarid,pbl_tke_glo(:,:,nsrf))
    937918#endif
    938919      ENDDO
Note: See TracChangeset for help on using the changeset viewer.