Ignore:
Timestamp:
May 24, 2002, 4:10:05 PM (22 years ago)
Author:
lmdzadmin
Message:

Les temperatures de sol en sortie selon les surfaces n'étaient pas bonnes
dans le fichier histday
Pb de date initiale pour les variables en "once"
LF

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ.3.3/branches/rel-LF/libf/phylmd/physiq.F

    r353 r361  
    16571657
    16581658c$$$PB Positionner date0 pour initialisation de ORCHIDEE
    1659 c$$$      date0 = zjulian
    1660       date0 = day_ini
     1659      date0 = zjulian
     1660C      date0 = day_ini
    16611661      WRITE(*,*) 'physiq date0 : ',date0
    16621662c
     
    24032403         i = NINT(zout/zsto)
    24042404         CALL gr_fi_ecrit(1,klon,iim,jjmp1,pphis,zx_tmp_2d)
    2405          CALL histwrite(nid_day,"phis",i,zx_tmp_2d,iim*jjmp1,ndex2d)
     2405       CALL histwrite(nid_day,"phis",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    24062406         varname = 'phis'
    24072407         vartitle= 'Surface geop. height'
     
    24122412         i = NINT(zout/zsto)
    24132413         CALL gr_fi_ecrit(1,klon,iim,jjmp1,paire,zx_tmp_2d)
    2414          CALL histwrite(nid_day,"aire",i,zx_tmp_2d,iim*jjmp1,ndex2d)
     2414       CALL histwrite(nid_day,"aire",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    24152415         varname = 'aire'
    24162416         vartitle= 'Grid area'
     
    24322432C
    24332433      zx_tmp_fi2d(1 : klon) = ftsol(1 : klon, is_lic)
    2434       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxtsol,zx_tmp_2d)
     2434      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    24352435      CALL histwrite(nid_day,"tlic",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    24362436c     call writephy(fid_day,prof2d_av,'tlic',ftsol(1 : klon, is_lic),
     
    24382438C
    24392439      zx_tmp_fi2d(1 : klon) = ftsol(1 : klon, is_oce)
    2440       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxtsol,zx_tmp_2d)
     2440      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    24412441      CALL histwrite(nid_day,"toce",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    24422442c     call writephy(fid_day,prof2d_av,'toce',ftsol(1 : klon, is_oce),
     
    24442444C
    24452445      zx_tmp_fi2d(1 : klon) = ftsol(1 : klon, is_sic)
    2446       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxtsol,zx_tmp_2d)
     2446      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    24472447      CALL histwrite(nid_day,"tsic",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    24482448c     call writephy(fid_day,prof2d_av,'tsic',ftsol(1 : klon, is_sic),
     
    26962696         itau_w = itau_phy + itap
    26972697
    2698          i = NINT(zout/zsto)
    2699          CALL gr_fi_ecrit(1,klon,iim,jjmp1,pphis,zx_tmp_2d)
    2700          CALL histwrite(nid_mth,"phis",i,zx_tmp_2d,iim*jjmp1,ndex2d)
    2701 C
    2702          i = NINT(zout/zsto)
    2703          CALL gr_fi_ecrit(1,klon,iim,jjmp1,paire,zx_tmp_2d)
    2704          CALL histwrite(nid_mth,"aire",i,zx_tmp_2d,iim*jjmp1,ndex2d)
     2698      i = NINT(zout/zsto)
     2699      CALL gr_fi_ecrit(1,klon,iim,jjmp1,pphis,zx_tmp_2d)
     2700      CALL histwrite(nid_mth,"phis",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     2701C
     2702      i = NINT(zout/zsto)
     2703      CALL gr_fi_ecrit(1,klon,iim,jjmp1,paire,zx_tmp_2d)
     2704      CALL histwrite(nid_mth,"aire",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    27052705
    27062706      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxtsol,zx_tmp_2d)
     
    31043104
    31053105         i = NINT(zout/zsto)
    3106          CALL gr_fi_ecrit(1,klon,iim,jjmp1,pphis,zx_tmp_2d)
    3107          CALL histwrite(nid_ins,"phis",i,zx_tmp_2d,iim*jjmp1,ndex2d)
     3106      CALL gr_fi_ecrit(1,klon,iim,jjmp1,pphis,zx_tmp_2d)
     3107      CALL histwrite(nid_ins,"phis",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    31083108c
    31093109         i = NINT(zout/zsto)
    3110          CALL gr_fi_ecrit(1,klon,iim,jjmp1,paire,zx_tmp_2d)
    3111          CALL histwrite(nid_ins,"aire",i,zx_tmp_2d,iim*jjmp1,ndex2d)
     3110      CALL gr_fi_ecrit(1,klon,iim,jjmp1,paire,zx_tmp_2d)
     3111      CALL histwrite(nid_ins,"aire",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    31123112
    31133113      DO i = 1, klon
Note: See TracChangeset for help on using the changeset viewer.