Changeset 825 for LMDZ4/trunk


Ignore:
Timestamp:
Aug 3, 2007, 3:08:17 PM (17 years ago)
Author:
Laurent Fairhead
Message:

Probleme de calcul des pas de temps d'ecriture dans le fichier paramLMDZ_dyn.nc
LF

Location:
LMDZ4/trunk/libf/dyn3d
Files:
2 edited

Legend:

Unmodified
Added
Removed

  • LMDZ4/trunk/libf/dyn3d/leapfrog.F

    r773 r825  
    185185      PARAMETER (testita = 9)
    186186
     187      integer itau_w   ! pas de temps ecriture = itap + itau_phy
     188
     189
    187190      itaufin   = nday*day_step
    188191      itaufinp1 = itaufin +1

  • LMDZ4/trunk/libf/dyn3d/write_paramLMDZ_dyn.h

    r792 r825  
    11c
    22      ndex2d = 0
    3       itau_dyn=itau_dyn+itau
     3      itau_w=itau_dyn+itau
    44c
    55      zx_tmp_2d(1:iip1,1:jjp1)=FLOAT(prt_level)
    6       CALL histwrite(nid_ctesGCM, "prt_level", itau_dyn,
     6      CALL histwrite(nid_ctesGCM, "prt_level", itau_w,
    77     .               zx_tmp_2d,iip1*jjp1,ndex2d)
    88c
    99      zx_tmp_2d(1:iip1,1:jjp1)=FLOAT(dayref)
    10       CALL histwrite(nid_ctesGCM, "dayref", itau_dyn,
     10      CALL histwrite(nid_ctesGCM, "dayref", itau_w,
    1111     .               zx_tmp_2d,iip1*jjp1,ndex2d)
    1212c
    1313      zx_tmp_2d(1:iip1,1:jjp1)=FLOAT(anneeref)
    14       CALL histwrite(nid_ctesGCM, "anneeref", itau_dyn,
     14      CALL histwrite(nid_ctesGCM, "anneeref", itau_w,
    1515     .               zx_tmp_2d,iip1*jjp1,ndex2d)
    1616c
    1717      zx_tmp_2d(1:iip1,1:jjp1)=FLOAT(raz_date)
    18       CALL histwrite(nid_ctesGCM, "raz_date", itau_dyn,
     18      CALL histwrite(nid_ctesGCM, "raz_date", itau_w,
    1919     .               zx_tmp_2d,iip1*jjp1,ndex2d)
    2020c
    2121      zx_tmp_2d(1:iip1,1:jjp1)=FLOAT(nday)
    22       CALL histwrite(nid_ctesGCM, "nday", itau_dyn,
     22      CALL histwrite(nid_ctesGCM, "nday", itau_w,
    2323     .               zx_tmp_2d,iip1*jjp1,ndex2d)
    2424c
    2525      zx_tmp_2d(1:iip1,1:jjp1)=FLOAT(day_step)
    26       CALL histwrite(nid_ctesGCM, "day_step", itau_dyn,
     26      CALL histwrite(nid_ctesGCM, "day_step", itau_w,
    2727     .               zx_tmp_2d,iip1*jjp1,ndex2d)
    2828c
    2929      zx_tmp_2d(1:iip1,1:jjp1)=FLOAT(iperiod)
    30       CALL histwrite(nid_ctesGCM, "iperiod", itau_dyn,
     30      CALL histwrite(nid_ctesGCM, "iperiod", itau_w,
    3131     .               zx_tmp_2d,iip1*jjp1,ndex2d)
    3232c
    3333      zx_tmp_2d(1:iip1,1:jjp1)=FLOAT(iapp_tracvl)
    34       CALL histwrite(nid_ctesGCM, "iapp_tracvl", itau_dyn,
     34      CALL histwrite(nid_ctesGCM, "iapp_tracvl", itau_w,
    3535     .               zx_tmp_2d,iip1*jjp1,ndex2d)
    3636c
    3737      zx_tmp_2d(1:iip1,1:jjp1)=FLOAT(iconser)
    38       CALL histwrite(nid_ctesGCM, "iconser", itau_dyn,
     38      CALL histwrite(nid_ctesGCM, "iconser", itau_w,
    3939     .               zx_tmp_2d,iip1*jjp1,ndex2d)
    4040c
    4141      zx_tmp_2d(1:iip1,1:jjp1)=FLOAT(iecri)
    42       CALL histwrite(nid_ctesGCM, "iecri", itau_dyn,
     42      CALL histwrite(nid_ctesGCM, "iecri", itau_w,
    4343     .               zx_tmp_2d,iip1*jjp1,ndex2d)
    4444c
    4545      zx_tmp_2d(1:iip1,1:jjp1)=periodav
    46       CALL histwrite(nid_ctesGCM, "periodav", itau_dyn,
     46      CALL histwrite(nid_ctesGCM, "periodav", itau_w,
    4747     .               zx_tmp_2d,iip1*jjp1,ndex2d)
    4848c
    4949      zx_tmp_2d(1:iip1,1:jjp1)=FLOAT(idissip)
    50       CALL histwrite(nid_ctesGCM, "idissip", itau_dyn,
     50      CALL histwrite(nid_ctesGCM, "idissip", itau_w,
    5151     .               zx_tmp_2d,iip1*jjp1,ndex2d)
    5252c
    … …  
    5656       zx_tmp_2d(1:iip1,1:jjp1)=0.
    5757      ENDIF
    58       CALL histwrite(nid_ctesGCM, "lstardis", itau_dyn,
     58      CALL histwrite(nid_ctesGCM, "lstardis", itau_w,
    5959     .               zx_tmp_2d,iip1*jjp1,ndex2d)
    6060c
    6161      zx_tmp_2d(1:iip1,1:jjp1)=FLOAT(nitergdiv)
    62       CALL histwrite(nid_ctesGCM, "nitergdiv", itau_dyn,
     62      CALL histwrite(nid_ctesGCM, "nitergdiv", itau_w,
    6363     .               zx_tmp_2d,iip1*jjp1,ndex2d)
    6464c
    6565      zx_tmp_2d(1:iip1,1:jjp1)=FLOAT(nitergrot)
    66       CALL histwrite(nid_ctesGCM, "nitergrot", itau_dyn,
     66      CALL histwrite(nid_ctesGCM, "nitergrot", itau_w,
    6767     .               zx_tmp_2d,iip1*jjp1,ndex2d)
    6868c
    6969      zx_tmp_2d(1:iip1,1:jjp1)=FLOAT(niterh)
    70       CALL histwrite(nid_ctesGCM, "niterh", itau_dyn,
     70      CALL histwrite(nid_ctesGCM, "niterh", itau_w,
    7171     .               zx_tmp_2d,iip1*jjp1,ndex2d)
    7272c
    7373      zx_tmp_2d(1:iip1,1:jjp1)=tetagdiv
    74       CALL histwrite(nid_ctesGCM, "tetagdiv", itau_dyn,
     74      CALL histwrite(nid_ctesGCM, "tetagdiv", itau_w,
    7575     .               zx_tmp_2d,iip1*jjp1,ndex2d)
    7676c
    7777      zx_tmp_2d(1:iip1,1:jjp1)=tetagrot
    78       CALL histwrite(nid_ctesGCM, "tetagrot", itau_dyn,
     78      CALL histwrite(nid_ctesGCM, "tetagrot", itau_w,
    7979     .               zx_tmp_2d,iip1*jjp1,ndex2d)
    8080c
    8181      zx_tmp_2d(1:iip1,1:jjp1)=tetatemp
    82       CALL histwrite(nid_ctesGCM, "tetatemp", itau_dyn,
     82      CALL histwrite(nid_ctesGCM, "tetatemp", itau_w,
    8383     .               zx_tmp_2d,iip1*jjp1,ndex2d)
    8484c
    8585      zx_tmp_2d(1:iip1,1:jjp1)=coefdis
    86       CALL histwrite(nid_ctesGCM, "coefdis", itau_dyn,
     86      CALL histwrite(nid_ctesGCM, "coefdis", itau_w,
    8787     .               zx_tmp_2d,iip1*jjp1,ndex2d)
    8888c
    … …  
    9292       zx_tmp_2d(1:iip1,1:jjp1)=0.
    9393      ENDIF
    94       CALL histwrite(nid_ctesGCM, "purmats", itau_dyn,
     94      CALL histwrite(nid_ctesGCM, "purmats", itau_w,
    9595     .               zx_tmp_2d,iip1*jjp1,ndex2d)
    9696c
    … …  
    100100       zx_tmp_2d(1:iip1,1:jjp1)=0.
    101101      ENDIF
    102       CALL histwrite(nid_ctesGCM, "ok_guide", itau_dyn,
     102      CALL histwrite(nid_ctesGCM, "ok_guide", itau_w,
    103103     .               zx_tmp_2d,iip1*jjp1,ndex2d)
    104104c
    … …  
    108108       zx_tmp_2d(1:iip1,1:jjp1)=0.
    109109      ENDIF
    110       CALL histwrite(nid_ctesGCM, "true_calendar", itau_dyn,
     110      CALL histwrite(nid_ctesGCM, "true_calendar", itau_w,
    111111     .               zx_tmp_2d,iip1*jjp1,ndex2d)
    112112c
    113113      zx_tmp_2d(1:iip1,1:jjp1)=FLOAT(iflag_phys)
    114       CALL histwrite(nid_ctesGCM, "iflag_phys", itau_dyn,
     114      CALL histwrite(nid_ctesGCM, "iflag_phys", itau_w,
    115115     .               zx_tmp_2d,iip1*jjp1,ndex2d)
    116116c
    117117      zx_tmp_2d(1:iip1,1:jjp1)=FLOAT(iphysiq)
    118       CALL histwrite(nid_ctesGCM, "iphysiq", itau_dyn,
     118      CALL histwrite(nid_ctesGCM, "iphysiq", itau_w,
    119119     .               zx_tmp_2d,iip1*jjp1,ndex2d)
    120120c
    121121      zx_tmp_2d(1:iip1,1:jjp1)=ecritphy
    122       CALL histwrite(nid_ctesGCM, "ecritphy", itau_dyn,
     122      CALL histwrite(nid_ctesGCM, "ecritphy", itau_w,
    123123     .               zx_tmp_2d,iip1*jjp1,ndex2d)
    124124c
    … …  
    128128       zx_tmp_2d(1:iip1,1:jjp1)=0.
    129129      ENDIF
    130       CALL histwrite(nid_ctesGCM, "cycle_diurne", itau_dyn,
     130      CALL histwrite(nid_ctesGCM, "cycle_diurne", itau_w,
    131131     .               zx_tmp_2d,iip1*jjp1,ndex2d)
    132132c
    … …  
    136136       zx_tmp_2d(1:iip1,1:jjp1)=0.
    137137      ENDIF
    138       CALL histwrite(nid_ctesGCM, "soil_model", itau_dyn,
     138      CALL histwrite(nid_ctesGCM, "soil_model", itau_w,
    139139     .               zx_tmp_2d,iip1*jjp1,ndex2d)
    140140c
    … …  
    144144       zx_tmp_2d(1:iip1,1:jjp1)=0.
    145145      ENDIF
    146       CALL histwrite(nid_ctesGCM, "new_oliq", itau_dyn,
     146      CALL histwrite(nid_ctesGCM, "new_oliq", itau_w,
    147147     .               zx_tmp_2d,iip1*jjp1,ndex2d)
    148148c
    … …  
    152152       zx_tmp_2d(1:iip1,1:jjp1)=0.
    153153      ENDIF
    154       CALL histwrite(nid_ctesGCM, "ok_orodr", itau_dyn,
     154      CALL histwrite(nid_ctesGCM, "ok_orodr", itau_w,
    155155     .               zx_tmp_2d,iip1*jjp1,ndex2d)
    156156c
    … …  
    160160       zx_tmp_2d(1:iip1,1:jjp1)=0.
    161161      ENDIF
    162       CALL histwrite(nid_ctesGCM, "ok_orolf", itau_dyn,
     162      CALL histwrite(nid_ctesGCM, "ok_orolf", itau_w,
    163163     .               zx_tmp_2d,iip1*jjp1,ndex2d)
    164164c
    … …  
    168168       zx_tmp_2d(1:iip1,1:jjp1)=0.
    169169      ENDIF
    170       CALL histwrite(nid_ctesGCM, "ok_limitvrai", itau_dyn,
     170      CALL histwrite(nid_ctesGCM, "ok_limitvrai", itau_w,
    171171     .               zx_tmp_2d,iip1*jjp1,ndex2d)
    172172c
    173173      zx_tmp_2d(1:iip1,1:jjp1)=nbapp_rad
    174       CALL histwrite(nid_ctesGCM, "nbapp_rad", itau_dyn,
     174      CALL histwrite(nid_ctesGCM, "nbapp_rad", itau_w,
    175175     .               zx_tmp_2d,iip1*jjp1,ndex2d)
    176176c
    177177      zx_tmp_2d(1:iip1,1:jjp1)=iflag_con
    178       CALL histwrite(nid_ctesGCM, "iflag_con", itau_dyn,
     178      CALL histwrite(nid_ctesGCM, "iflag_con", itau_w,
    179179     .               zx_tmp_2d,iip1*jjp1,ndex2d)
    180180c
    181181      zx_tmp_2d(1:iip1,1:jjp1)=clon
    182       CALL histwrite(nid_ctesGCM, "clon", itau_dyn,
     182      CALL histwrite(nid_ctesGCM, "clon", itau_w,
    183183     .               zx_tmp_2d,iip1*jjp1,ndex2d)
    184184c
    185185      zx_tmp_2d(1:iip1,1:jjp1)=clat
    186       CALL histwrite(nid_ctesGCM, "clat", itau_dyn,
     186      CALL histwrite(nid_ctesGCM, "clat", itau_w,
    187187     .               zx_tmp_2d,iip1*jjp1,ndex2d)
    188188c
    189189      zx_tmp_2d(1:iip1,1:jjp1)=grossismx
    190       CALL histwrite(nid_ctesGCM, "grossismx", itau_dyn,
     190      CALL histwrite(nid_ctesGCM, "grossismx", itau_w,
    191191     .               zx_tmp_2d,iip1*jjp1,ndex2d)
    192192c
    193193      zx_tmp_2d(1:iip1,1:jjp1)=grossismy
    194       CALL histwrite(nid_ctesGCM, "grossismy", itau_dyn,
     194      CALL histwrite(nid_ctesGCM, "grossismy", itau_w,
    195195     .               zx_tmp_2d,iip1*jjp1,ndex2d)
    196196c
    … …  
    200200       zx_tmp_2d(1:iip1,1:jjp1)=0.
    201201      ENDIF
    202       CALL histwrite(nid_ctesGCM, "fxyhypb", itau_dyn,
     202      CALL histwrite(nid_ctesGCM, "fxyhypb", itau_w,
    203203     .               zx_tmp_2d,iip1*jjp1,ndex2d)
    204204c
    205205      zx_tmp_2d(1:iip1,1:jjp1)=dzoomx
    206       CALL histwrite(nid_ctesGCM, "dzoomx", itau_dyn,
     206      CALL histwrite(nid_ctesGCM, "dzoomx", itau_w,
    207207     .               zx_tmp_2d,iip1*jjp1,ndex2d)
    208208c
    209209      zx_tmp_2d(1:iip1,1:jjp1)=dzoomy
    210       CALL histwrite(nid_ctesGCM, "dzoomy", itau_dyn,
     210      CALL histwrite(nid_ctesGCM, "dzoomy", itau_w,
    211211     .               zx_tmp_2d,iip1*jjp1,ndex2d)
    212212c
    213213      zx_tmp_2d(1:iip1,1:jjp1)=taux
    214       CALL histwrite(nid_ctesGCM, "taux", itau_dyn,
     214      CALL histwrite(nid_ctesGCM, "taux", itau_w,
    215215     .               zx_tmp_2d,iip1*jjp1,ndex2d)
    216216c
    217217      zx_tmp_2d(1:iip1,1:jjp1)=tauy
    218       CALL histwrite(nid_ctesGCM, "tauy", itau_dyn,
     218      CALL histwrite(nid_ctesGCM, "tauy", itau_w,
    219219     .               zx_tmp_2d,iip1*jjp1,ndex2d)
    220220c
    … …  
    224224       zx_tmp_2d(1:iip1,1:jjp1)=0.
    225225      ENDIF
    226       CALL histwrite(nid_ctesGCM, "ysinus", itau_dyn,
     226      CALL histwrite(nid_ctesGCM, "ysinus", itau_w,
    227227     .               zx_tmp_2d,iip1*jjp1,ndex2d)
    228228c
    229229      zx_tmp_2d(1:iip1,1:jjp1)=ip_ebil_dyn
    230       CALL histwrite(nid_ctesGCM, "ip_ebil_dyn", itau_dyn,
     230      CALL histwrite(nid_ctesGCM, "ip_ebil_dyn", itau_w,
    231231     .               zx_tmp_2d,iip1*jjp1,ndex2d)
    232232c
Note: See TracChangeset for help on using the changeset viewer.