Ignore:
Timestamp:
Aug 21, 2015, 9:23:13 AM (9 years ago)
Author:
Ehouarn Millour
Message:

Physics/dynamics separation: get rid of all the 'include "temps.h"' in the physics; variables in module time_phylmdz_mod must be used instead. Also added JD_cur, JH_cur and JD_ref in module phys_cal_mod, in preparation for having physics handle its calendar internally.
EM

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ5/trunk/libf/phylmd/write_histday_seri.h

    r2343 r2344  
    2424      zx_tmp_fi2d(1:klon)=moyglo
    2525!
    26       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     26      CALL gr_fi_ecrit(1, klon,nbp_lon,nbp_lat, zx_tmp_fi2d,zx_tmp_2d)
    2727      CALL histwrite(nid_day_seri,"bilTOA",itau_w, &
    28                      zx_tmp_2d,iim*jjmp1,ndex2d)
     28                     zx_tmp_2d,nbp_lon*nbp_lat,ndex2d)
    2929!
    3030      ok_msk=.FALSE.
     
    3333      zx_tmp_fi2d(1:klon)=moyglo
    3434!
    35       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     35      CALL gr_fi_ecrit(1, klon,nbp_lon,nbp_lat, zx_tmp_fi2d,zx_tmp_2d)
    3636      CALL histwrite(nid_day_seri,"bils",itau_w, &
    37                      zx_tmp_2d,iim*jjmp1,ndex2d)
     37                     zx_tmp_2d,nbp_lon*nbp_lat,ndex2d)
    3838!
    3939      DO k=1, klev
     
    4848      zx_tmp_fi2d(1:klon)=moyglo
    4949!
    50       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     50      CALL gr_fi_ecrit(1, klon,nbp_lon,nbp_lat, zx_tmp_fi2d,zx_tmp_2d)
    5151      CALL histwrite(nid_day_seri,"ecin",itau_w, &
    52                      zx_tmp_2d,iim*jjmp1,ndex2d)
     52                     zx_tmp_2d,nbp_lon*nbp_lat,ndex2d)
    5353!
    5454!IM 151004 BEG
     
    6565      zx_tmp_fi2d(1:klon)=moyglo
    6666!
    67       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     67      CALL gr_fi_ecrit(1, klon,nbp_lon,nbp_lat, zx_tmp_fi2d,zx_tmp_2d)
    6868      CALL histwrite(nid_day_seri,"momang",itau_w,zx_tmp_2d, &
    69                      iim*jjmp1,ndex2d)
     69                     nbp_lon*nbp_lat,ndex2d)
    7070!
    7171! friction torque
     
    8080      zx_tmp_fi2d(1:klon)=moyglo
    8181!
    82       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     82      CALL gr_fi_ecrit(1, klon,nbp_lon,nbp_lat, zx_tmp_fi2d,zx_tmp_2d)
    8383      CALL histwrite(nid_day_seri,"frictor",itau_w,zx_tmp_2d, &
    84                      iim*jjmp1,ndex2d)
     84                     nbp_lon*nbp_lat,ndex2d)
    8585!
    8686! mountain torque
    8787!
    8888!IM 190504 BEG
    89       CALL gr_fi_dyn(1,klon,iim+1,jjm+1,airephy,airedyn)
    90       CALL gr_fi_dyn(klev+1,klon,iim+1,jjm+1,paprs,padyn)
    91       CALL gr_fi_dyn(1,klon,iim+1,jjm+1,rlat,rlatdyn)
     89      CALL gr_fi_dyn(1,klon,nbp_lon+1,nbp_lat,airephy,airedyn)
     90      CALL gr_fi_dyn(klev+1,klon,nbp_lon+1,nbp_lat,paprs,padyn)
     91      CALL gr_fi_dyn(1,klon,nbp_lon+1,nbp_lat,rlat,rlatdyn)
    9292      mountor=0.
    9393      airetot=0.
    94       DO j = 1, jjmp1
    95        DO i = 1, iim+1
    96         ij=i+(iim+1)*(j-1)
     94      DO j = 1, nbp_lat
     95       DO i = 1, nbp_lon+1
     96        ij=i+(nbp_lon+1)*(j-1)
    9797        zx_tmp(ij)=0.
    9898        DO k = 1, klev
     
    113113!
    114114!IM 190504 END
    115       zx_tmp_2d(1:iim,1:jjmp1)=mountor
     115      zx_tmp_2d(1:nbp_lon,1:nbp_lat)=mountor
    116116      CALL histwrite(nid_day_seri,"mountor",itau_w,zx_tmp_2d, &
    117                      iim*jjmp1,ndex2d)
     117                     nbp_lon*nbp_lat,ndex2d)
    118118!
    119119      ENDIF !(1.EQ.0) THEN
    120120!
    121121!
    122       CALL gr_fi_dyn(1,klon,iim+1,jjm+1,airephy,airedyn)
    123       CALL gr_fi_ecrit(1,klon,iim,jjmp1,airephy,zx_tmp_2d)
     122      CALL gr_fi_dyn(1,klon,nbp_lon+1,nbp_lat,airephy,airedyn)
     123      CALL gr_fi_ecrit(1,klon,nbp_lon,nbp_lat,airephy,zx_tmp_2d)
    124124      airetot=0.
    125 !     DO j = 1, jjmp1
    126 !      DO i = 1, iim+1
    127 !       ij=i+(iim+1)*(j-1)
     125!     DO j = 1, nbp_lat
     126!      DO i = 1, nbp_lon+1
     127!       ij=i+(nbp_lon+1)*(j-1)
    128128!       DO k = 1, klev
    129129!        airetot=airetot+airedyn(i,j)
     
    139139!
    140140      airetot=0.
    141       DO j=1, jjmp1
    142        DO i=1, iim
     141      DO j=1, nbp_lat
     142       DO i=1, nbp_lon
    143143        airetot=airetot+zx_tmp_2d(i,j)
    144144       ENDDO
     
    149149!
    150150      zx_tmp_fi2d(1:klon)=aam/airetot
    151       CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
     151      CALL gr_fi_ecrit(1,klon,nbp_lon,nbp_lat,zx_tmp_fi2d,zx_tmp_2d)
    152152      CALL histwrite(nid_day_seri,"momang",itau_w,zx_tmp_2d, &
    153                      iim*jjmp1,ndex2d)
     153                     nbp_lon*nbp_lat,ndex2d)
    154154!
    155155      zx_tmp_fi2d(1:klon)=torsfc/airetot
    156       CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
     156      CALL gr_fi_ecrit(1,klon,nbp_lon,nbp_lat,zx_tmp_fi2d,zx_tmp_2d)
    157157      CALL histwrite(nid_day_seri,"torsfc",itau_w,zx_tmp_2d, &
    158                      iim*jjmp1,ndex2d)
     158                     nbp_lon*nbp_lat,ndex2d)
    159159!
    160160!IM 151004 END
     
    164164      zx_tmp_fi2d(1:klon)=moyglo
    165165!
    166       CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
     166      CALL gr_fi_ecrit(1,klon,nbp_lon,nbp_lat,zx_tmp_fi2d,zx_tmp_2d)
    167167      CALL histwrite(nid_day_seri,"tamv",itau_w, &
    168                      zx_tmp_2d,iim*jjmp1,ndex2d)
     168                     zx_tmp_2d,nbp_lon*nbp_lat,ndex2d)
    169169!
    170170      ok_msk=.FALSE.
     
    173173      zx_tmp_fi2d(1:klon)=moyglo
    174174!
    175       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     175      CALL gr_fi_ecrit(1, klon,nbp_lon,nbp_lat, zx_tmp_fi2d,zx_tmp_2d)
    176176      CALL histwrite(nid_day_seri,"psol",itau_w, &
    177                      zx_tmp_2d,iim*jjmp1,ndex2d)
     177                     zx_tmp_2d,nbp_lon*nbp_lat,ndex2d)
    178178!
    179179      ok_msk=.FALSE.
     
    182182      zx_tmp_fi2d(1:klon)=moyglo
    183183!
    184       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     184      CALL gr_fi_ecrit(1, klon,nbp_lon,nbp_lat, zx_tmp_fi2d,zx_tmp_2d)
    185185      CALL histwrite(nid_day_seri,"evap",itau_w, &
    186                      zx_tmp_2d,iim*jjmp1,ndex2d)
     186                     zx_tmp_2d,nbp_lon*nbp_lat,ndex2d)
    187187!
    188188!     DO i=1, klon
     
    196196!     zx_tmp_fi2d(1:klon)=moyglo
    197197!
    198 !     CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
     198!     CALL gr_fi_ecrit(1, klon,nbp_lon,nbp_lat,zx_tmp_fi2d,zx_tmp_2d)
    199199!     CALL histwrite(nid_day_seri,"SnowFrac",
    200 !    .               itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     200!    .               itau_w,zx_tmp_2d,nbp_lon*nbp_lat,ndex2d)
    201201!
    202202!     DO i=1, klon
     
    212212!     zx_tmp_fi2d(1:klon)=moyglo
    213213!
    214 !     CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
     214!     CALL gr_fi_ecrit(1, klon,nbp_lon,nbp_lat,zx_tmp_fi2d,zx_tmp_2d)
    215215!     CALL histwrite(nid_day_seri,"snow_depth",itau_w,
    216 !    .               zx_tmp_2d,iim*jjmp1,ndex2d)
     216!    .               zx_tmp_2d,nbp_lon*nbp_lat,ndex2d)
    217217!
    218218      DO i=1, klon
     
    226226      zx_tmp_fi2d(1:klon)=moyglo
    227227!
    228       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
     228      CALL gr_fi_ecrit(1, klon,nbp_lon,nbp_lat, zx_tmp_fi2d, zx_tmp_2d)
    229229      CALL histwrite(nid_day_seri,"tsol_"//clnsurf(is_oce), &
    230                      itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     230                     itau_w,zx_tmp_2d,nbp_lon*nbp_lat,ndex2d)
    231231!
    232232!=================================================================
Note: See TracChangeset for help on using the changeset viewer.