1 | c IF (ok_mensuelNMC) THEN |
---|
2 | IF (ok_mensuel) THEN |
---|
3 | c |
---|
4 | ndex2d = 0 |
---|
5 | ndex3d = 0 |
---|
6 | c |
---|
7 | c Champs 2D: |
---|
8 | c |
---|
9 | |
---|
10 | zsto = dtime |
---|
11 | zout = dtime * ecrit_mth |
---|
12 | itau_w = itau_phy + itap |
---|
13 | |
---|
14 | i = NINT(zout/zsto) |
---|
15 | c |
---|
16 | c Champs interpolles sur des niveaux de pression du NMC |
---|
17 | c110304 BEG |
---|
18 | DO k=1, nlevSTD |
---|
19 | c |
---|
20 | bb=clevSTD(k) |
---|
21 | c |
---|
22 | IF(k.GE.2) THEN |
---|
23 | aa=clevSTD(k) |
---|
24 | bb=aa(1:lnblnk1(aa)) |
---|
25 | ENDIF |
---|
26 | c |
---|
27 | CALL gr_fi_ecrit(1, klon,iim,jjmp1, tlevSTD(:,k),zx_tmp_2d) |
---|
28 | CALL histwrite(nid_nmc,"t"//bb,itau_w,zx_tmp_2d, |
---|
29 | $ iim*jjmp1,ndex2d) |
---|
30 | c |
---|
31 | CALL gr_fi_ecrit(1, klon,iim,jjmp1, philevSTD(:,k),zx_tmp_2d) |
---|
32 | CALL histwrite(nid_nmc,"phi"//bb,itau_w,zx_tmp_2d, |
---|
33 | $ iim*jjmp1,ndex2d) |
---|
34 | c |
---|
35 | CALL gr_fi_ecrit(1, klon,iim,jjmp1, qlevSTD(:,k),zx_tmp_2d) |
---|
36 | CALL histwrite(nid_nmc,"q"//bb,itau_w,zx_tmp_2d, |
---|
37 | $ iim*jjmp1,ndex2d) |
---|
38 | c |
---|
39 | CALL gr_fi_ecrit(1, klon,iim,jjmp1, rhlevSTD(:,k),zx_tmp_2d) |
---|
40 | CALL histwrite(nid_nmc,"rh"//bb,itau_w,zx_tmp_2d, |
---|
41 | $ iim*jjmp1,ndex2d) |
---|
42 | c |
---|
43 | CALL gr_fi_ecrit(1, klon,iim,jjmp1, ulevSTD(:,k),zx_tmp_2d) |
---|
44 | CALL histwrite(nid_nmc,"u"//bb,itau_w,zx_tmp_2d, |
---|
45 | $ iim*jjmp1,ndex2d) |
---|
46 | c |
---|
47 | CALL gr_fi_ecrit(1, klon,iim,jjmp1, vlevSTD(:,k),zx_tmp_2d) |
---|
48 | CALL histwrite(nid_nmc,"v"//bb,itau_w,zx_tmp_2d, |
---|
49 | $ iim*jjmp1,ndex2d) |
---|
50 | c |
---|
51 | ENDDO !k=1, nlevSTD |
---|
52 | ccc |
---|
53 | c110304 END |
---|
54 | c |
---|
55 | if (ok_sync) then |
---|
56 | call histsync(nid_nmc) |
---|
57 | endif |
---|
58 | |
---|
59 | ENDIF |
---|