Ignore:
Timestamp:
May 25, 2005, 3:10:09 PM (20 years ago)
Author:
Laurent Fairhead
Message:

Synchronisation avec tous les diagnostiques de Ionela IM
Inclusion du slab ocean IM
LF

File:
1 edited

Legend:

Unmodified
Added
Removed

  • LMDZ4/trunk/libf/phylmd/write_histmthNMC.h

    r524 r644  
    55      IF (ok_mensuel) THEN
    66c
    7       ndex2d = 0
    8       ndex3d = 0
     7       ndex3d = 0
     8       itau_w = itau_phy + itap
     9ccc
     10c  Champs interpolles sur des niveaux de pression du NMC
    911c
    10 c Champs 2D:
     12c     PARAMETER(nout=3) !nout=1 : day; =2 : mth; =3 : NMC
     13ccc
     14      IF(type_run.EQ."CLIM".OR.type_run.EQ."ENSP") THEN
     15ccc
     16       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,tsumSTD(:,:,2),
     17     $      zx_tmp_3d)
     18       CALL histwrite(nid_nmc,"temp",itau_w,zx_tmp_3d,
     19     $       iim*jjmp1*nlevSTD,ndex3d)
    1120c
    12 
    13         zsto = dtime
    14         zout = dtime * ecrit_mth
    15         itau_w = itau_phy + itap
    16 
    17       i = NINT(zout/zsto)
     21       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,phisumSTD(:,:,2),
     22     $     zx_tmp_3d)
     23       CALL histwrite(nid_nmc,"phi",itau_w,zx_tmp_3d,
     24     $       iim*jjmp1*nlevSTD,ndex3d)
    1825c
    19 c  Champs interpolles sur des niveaux de pression du NMC
    20 c110304 BEG
    21       DO k=1, nlevSTD
     26       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,qsumSTD(:,:,2),
     27     $     zx_tmp_3d)
     28       CALL histwrite(nid_nmc,"q",itau_w,zx_tmp_3d,
     29     $       iim*jjmp1*nlevSTD,ndex3d)
    2230c
    23          bb=clevSTD(k)
     31       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,rhsumSTD(:,:,2),
     32     $     zx_tmp_3d)
     33       CALL histwrite(nid_nmc,"rh",itau_w,zx_tmp_3d,
     34     $       iim*jjmp1*nlevSTD,ndex3d)
    2435c
    25          IF(k.GE.2) THEN
    26           aa=clevSTD(k)
    27           bb=aa(1:lnblnk1(aa))
     36       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,usumSTD(:,:,2),
     37     $     zx_tmp_3d)
     38       CALL histwrite(nid_nmc,"u",itau_w,zx_tmp_3d,
     39     $       iim*jjmp1*nlevSTD,ndex3d)
     40c
     41       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,vsumSTD(:,:,2),
     42     $     zx_tmp_3d)
     43       CALL histwrite(nid_nmc,"v",itau_w,zx_tmp_3d,
     44     $       iim*jjmp1*nlevSTD,ndex3d)
     45ccc
     46      ELSE IF(type_run.EQ."AMIP".OR.type_run.EQ."CFMI") THEN
     47ccc
     48       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,tsumSTD(:,:,3),
     49     $     zx_tmp_3d)
     50       CALL histwrite(nid_nmc,"temp",itau_w,zx_tmp_3d,
     51     $       iim*jjmp1*nlevSTD,ndex3d)
     52c
     53       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,phisumSTD(:,:,3),
     54     $     zx_tmp_3d)
     55       CALL histwrite(nid_nmc,"phi",itau_w,zx_tmp_3d,
     56     $       iim*jjmp1*nlevSTD,ndex3d)
     57c
     58       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,qsumSTD(:,:,3),
     59     $     zx_tmp_3d)
     60       CALL histwrite(nid_nmc,"q",itau_w,zx_tmp_3d,
     61     $       iim*jjmp1*nlevSTD,ndex3d)
     62c
     63       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,rhsumSTD(:,:,3),
     64     $     zx_tmp_3d)
     65       CALL histwrite(nid_nmc,"rh",itau_w,zx_tmp_3d,
     66     $       iim*jjmp1*nlevSTD,ndex3d)
     67c
     68       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,usumSTD(:,:,3),
     69     $     zx_tmp_3d)
     70       CALL histwrite(nid_nmc,"u",itau_w,zx_tmp_3d,
     71     $       iim*jjmp1*nlevSTD,ndex3d)
     72c
     73       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,vsumSTD(:,:,3),
     74     $     zx_tmp_3d)
     75       CALL histwrite(nid_nmc,"v",itau_w,zx_tmp_3d,
     76     $       iim*jjmp1*nlevSTD,ndex3d)
     77c
     78       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,wsumSTD(:,:,3),
     79     $     zx_tmp_3d)
     80       CALL histwrite(nid_nmc,"w",itau_w,zx_tmp_3d,
     81     $       iim*jjmp1*nlevSTD,ndex3d)
     82c
     83       DO k=1, nlevSTD
     84        DO i=1, klon
     85         IF(tnondef(i,k,3).NE.1.E+20) THEN
     86          zx_tmp_fi3d(i,k) = (100.*tnondef(i,k,3))/FLOAT(ecrit_hf2mth)
     87         ELSE
     88          zx_tmp_fi3d(i,k) = 1.E+20
    2889         ENDIF
     90        ENDDO
     91       ENDDO !k=1, nlevSTD
    2992c
    30       CALL gr_fi_ecrit(1, klon,iim,jjmp1, tlevSTD(:,k),zx_tmp_2d)
    31       CALL histwrite(nid_nmc,"t"//bb,itau_w,zx_tmp_2d,
    32      $       iim*jjmp1,ndex2d)
     93       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
     94       CALL histwrite(nid_nmc,"psbg",itau_w,zx_tmp_3d,
     95     $       iim*jjmp1*nlevSTD,ndex3d)
    3396c
    34       CALL gr_fi_ecrit(1, klon,iim,jjmp1, philevSTD(:,k),zx_tmp_2d)
    35       CALL histwrite(nid_nmc,"phi"//bb,itau_w,zx_tmp_2d,
    36      $       iim*jjmp1,ndex2d)
     97       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,uvsumSTD(:,:,3),
     98     $     zx_tmp_3d)
     99       CALL histwrite(nid_nmc,"uv",itau_w,zx_tmp_3d,
     100     $       iim*jjmp1*nlevSTD,ndex3d)
    37101c
    38       CALL gr_fi_ecrit(1, klon,iim,jjmp1, qlevSTD(:,k),zx_tmp_2d)
    39       CALL histwrite(nid_nmc,"q"//bb,itau_w,zx_tmp_2d,
    40      $       iim*jjmp1,ndex2d)
     102       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,vqsumSTD(:,:,3),
     103     $     zx_tmp_3d)
     104       CALL histwrite(nid_nmc,"vq",itau_w,zx_tmp_3d,
     105     $       iim*jjmp1*nlevSTD,ndex3d)
    41106c
    42       CALL gr_fi_ecrit(1, klon,iim,jjmp1, rhlevSTD(:,k),zx_tmp_2d)
    43       CALL histwrite(nid_nmc,"rh"//bb,itau_w,zx_tmp_2d,
    44      $       iim*jjmp1,ndex2d)
     107       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,vTsumSTD(:,:,3),
     108     $     zx_tmp_3d)
     109       CALL histwrite(nid_nmc,"vT",itau_w,zx_tmp_3d,
     110     $       iim*jjmp1*nlevSTD,ndex3d)
    45111c
    46       CALL gr_fi_ecrit(1, klon,iim,jjmp1, ulevSTD(:,k),zx_tmp_2d)
    47       CALL histwrite(nid_nmc,"u"//bb,itau_w,zx_tmp_2d,
    48      $       iim*jjmp1,ndex2d)
     112       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1, wqsumSTD(:,:,3),
     113     $     zx_tmp_3d)
     114       CALL histwrite(nid_nmc,"wq",itau_w,zx_tmp_3d,
     115     $       iim*jjmp1*nlevSTD,ndex3d)
    49116c
    50       CALL gr_fi_ecrit(1, klon,iim,jjmp1, vlevSTD(:,k),zx_tmp_2d)
    51       CALL histwrite(nid_nmc,"v"//bb,itau_w,zx_tmp_2d,
    52      $       iim*jjmp1,ndex2d)
     117       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,vphisumSTD(:,:,3),
     118     $     zx_tmp_3d)
     119       CALL histwrite(nid_nmc,"vphi",itau_w,zx_tmp_3d,
     120     $       iim*jjmp1*nlevSTD,ndex3d)
    53121c
    54       ENDDO !k=1, nlevSTD
    55 ccc
    56 c110304 END
     122       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,wTsumSTD(:,:,3),
     123     $     zx_tmp_3d)
     124       CALL histwrite(nid_nmc,"wT",itau_w,zx_tmp_3d,
     125     $       iim*jjmp1*nlevSTD,ndex3d)
     126c
     127       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,u2sumSTD(:,:,3),
     128     $     zx_tmp_3d)
     129       CALL histwrite(nid_nmc,"uxu",itau_w,zx_tmp_3d,
     130     $       iim*jjmp1*nlevSTD,ndex3d)
     131c
     132       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,v2sumSTD(:,:,3),
     133     $     zx_tmp_3d)
     134       CALL histwrite(nid_nmc,"vxv",itau_w,zx_tmp_3d,
     135     $       iim*jjmp1*nlevSTD,ndex3d)
     136c
     137       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,T2sumSTD(:,:,3),
     138     $     zx_tmp_3d)
     139       CALL histwrite(nid_nmc,"TxT",itau_w,zx_tmp_3d,
     140     $       iim*jjmp1*nlevSTD,ndex3d)
     141c
     142      ENDIF !type_run
    57143c
    58144      if (ok_sync) then
Note: See TracChangeset for help on using the changeset viewer.