source: LMDZ.3.3/branches/rel-LF/libf/phylmd/ini_histISCCP.h @ 5392

Last change on this file since 5392 was 511, checked in by lmdzadmin, 21 years ago

Enleve region "1d" des fichiers ini*.h
IM

  • Property svn:eol-style set to native
  • Property svn:executable set to *
  • Property svn:keywords set to Author Date Id Revision
File size: 1.9 KB
Line 
1      IF (ok_isccp) THEN
2c
3c pour les champs instantannes, il faut mettre la meme valeur pour
4c zout et zsto.
5c dtime est passe par ailleurs a histbeg
6c zsto = frequence de stockage des champs
7c zout = frequence d'ecriture des champs
8         zsto = dtime
9c
10c ecriture 8 fois par jour
11c       zout = dtime * REAL(NINT(86400./dtime*ecrit_isccp))
12c ecriture toutes les 2h (12 fois par jour)
13c       zout = dtime * 4.
14c ecriture toutes les 1/h (48 fois par jour)
15c       zout = dtime
16c
17c ecriture mensuelle
18        zout = dtime * ecrit_mth
19c
20        print*,'ISCCP zout,zsto=',zout,zsto
21c
22c       PRINT*, 'La frequence de sortie ISCCP est de ', ecrit_isccp
23c
24        idayref = day_ref
25        CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian)
26        write(*,*)'ISCCP ', itau_phy, zjulian
27c
28c
29c definition coordonnees lon,lat en globale
30c
31        CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon)
32        DO i = 1, iim
33          zx_lon(i,1) = rlon(i+1)
34          zx_lon(i,jjmp1) = rlon(i+1)
35        ENDDO
36
37        CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat)
38c
39         CALL histbeg("histISCCP.nc", iim,zx_lon(:,1),jjmp1,zx_lat(1,:),
40     .                 1, iim, 1, jjmp1,
41     .                 itau_phy, zjulian, dtime,
42     .                 nhori, nid_isccp)
43c
44         CALL histvert(nid_isccp, "cldtopres","Cloud Top Pressure","mb",
45     .                 lmaxm1, cldtopres, nvert,'down')
46c
47c variables a ecrire
48c
49         DO k=1, kmaxm1
50          CALL histdef(nid_isccp, "cldISCCP_"//taulev(k),
51     .                "LMDZ ISCCP cld", "%",
52     .                iim, jjmp1,nhori,lmaxm1,1,lmaxm1,nvert,32,
53     .                "ave(X)", zsto,zout)
54         ENDDO
55c
56         CALL histdef(nid_isccp, "nsunlit",
57     .                "Nb of calls with sunlit ", "%",
58     .                iim, jjmp1,nhori,1,1,1,-99,32,
59     .                "ave(X)", zsto,zout)
60c
61        CALL histend(nid_isccp)
62c
63      ENDIF ! ok_isccp
Note: See TracBrowser for help on using the repository browser.