source: LMDZ.3.3/branches/rel-LF/libf/phylmd/write_histhf.h @ 486

Last change on this file since 486 was 486, checked in by lmdzadmin, 20 years ago

Phasage avec la version de Ionela
IM/LF

  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 2.5 KB
Line 
1      if (ok_hf) then
2
3c   Comprendre comment marche el i=nint(zout/zsto)
4c
5c$$$      print*,'ACRITURE HF !!! ACRITURE HF !!! ACRITURE HF !!! '
6      ndex2d = 0
7      ndex3d = 0
8
9      itau_w =  itau_phy + itap
10c
11c
12c     i = NINT(zout/zsto)
13c     CALL gr_fi_ecrit(1,klon,iim,jjmp1,pphis,zx_tmp_2d)
14c     CALL histwrite(nid_hf,"phis",i,zx_tmp_2d,iim*jjmp1,ndex2d)
15c
16c     i = NINT(zout/zsto)
17c     CALL gr_fi_ecrit(1,klon,iim,jjmp1,paire,zx_tmp_2d)
18c     CALL histwrite(nid_hf,"aire",i,zx_tmp_2d,iim*jjmp1,ndex2d)
19C
20      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxtsol,zx_tmp_2d)
21      CALL histwrite(nid_hf,"tsol",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
22c
23      DO i = 1, klon
24         zx_tmp_fi2d(i) = paprs(i,1)
25      ENDDO
26      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
27      CALL histwrite(nid_hf,"psol",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
28c
29      DO i = 1, klon
30         zx_tmp_fi2d(i) = rain_fall(i) + snow_fall(i)
31      ENDDO
32      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
33      CALL histwrite(nid_hf,"rain",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
34
35      CALL gr_fi_ecrit(1, klon,iim,jjmp1, u850,zx_tmp_2d)
36      CALL histwrite(nid_hf,"u850",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
37
38      CALL gr_fi_ecrit(1, klon,iim,jjmp1, v850,zx_tmp_2d)
39      CALL histwrite(nid_hf,"v850",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
40
41      CALL gr_fi_ecrit(1, klon,iim,jjmp1, u500,zx_tmp_2d)
42      CALL histwrite(nid_hf,"u500",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
43
44      CALL gr_fi_ecrit(1, klon,iim,jjmp1, v500,zx_tmp_2d)
45      CALL histwrite(nid_hf,"v500",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
46
47      CALL gr_fi_ecrit(1, klon,iim,jjmp1, u200,zx_tmp_2d)
48      CALL histwrite(nid_hf,"u200",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
49
50      CALL gr_fi_ecrit(1, klon,iim,jjmp1, v200,zx_tmp_2d)
51      CALL histwrite(nid_hf,"v200",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
52
53      CALL gr_fi_ecrit(1, klon,iim,jjmp1, phi500,zx_tmp_2d)
54      CALL histwrite(nid_hf,"phi500",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
55
56cIM cf FH
57      CALL gr_fi_ecrit(1, klon,iim,jjmp1, u_seri(:,1),zx_tmp_2d)
58      CALL histwrite(nid_hf,"u1",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
59c
60      CALL gr_fi_ecrit(1, klon,iim,jjmp1, v_seri(:,1),zx_tmp_2d)
61      CALL histwrite(nid_hf,"v1",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
62c
63      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragm,zx_tmp_2d)
64      CALL histwrite(nid_hf,"cdrm",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
65c
66      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragh,zx_tmp_2d)
67      CALL histwrite(nid_hf,"cdrh",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
68
69      if (ok_sync) then
70        call histsync(nid_hf)
71      endif
72
73      endif
Note: See TracBrowser for help on using the repository browser.