source: LMDZ5/branches/IPSLCM5A2.1_ISO/libf/phyiso/write_histREGDYN.h @ 5452

Last change on this file since 5452 was 3331, checked in by acozic, 7 years ago

Add modification for isotopes

  • Property svn:executable set to *
File size: 2.0 KB
Line 
1!
2! $Header$
3!
4      if (ok_regdyn) then
5     
6      if (is_sequential) then
7
8
9      ndex3d = 0
10      itau_w = itau_phy + itap + start_time * day_step / iphysiq
11c
12       CALL histwrite(nid_regdyn,"hw1",itau_w,histoW(:,:,:,1),
13     &               kmaxm1*lmaxm1*iwmax,ndex3d)
14c
15       CALL histwrite(nid_regdyn,"nh1",itau_w,nhistoW(:,:,:,1),
16     &               kmaxm1*lmaxm1*iwmax,ndex3d)
17c
18       CALL histwrite(nid_regdyn,"nht1",itau_w,nhistoWt(:,:,:,1),
19     &               kmaxm1*lmaxm1*iwmax,ndex3d)
20c
21       CALL histwrite(nid_regdyn,"hw2",itau_w,histoW(:,:,:,2),
22     &               kmaxm1*lmaxm1*iwmax,ndex3d)
23c
24       CALL histwrite(nid_regdyn,"nh2",itau_w,nhistoW(:,:,:,2),
25     &               kmaxm1*lmaxm1*iwmax,ndex3d)
26c
27       CALL histwrite(nid_regdyn,"nht2",itau_w,nhistoWt(:,:,:,2),
28     &               kmaxm1*lmaxm1*iwmax,ndex3d)
29c
30       CALL histwrite(nid_regdyn,"hw3",itau_w,histoW(:,:,:,3),
31     &               kmaxm1*lmaxm1*iwmax,ndex3d)
32c
33       CALL histwrite(nid_regdyn,"nh3",itau_w,nhistoW(:,:,:,3),
34     &               kmaxm1*lmaxm1*iwmax,ndex3d)
35c
36       CALL histwrite(nid_regdyn,"nht3",itau_w,nhistoWt(:,:,:,3),
37     &               kmaxm1*lmaxm1*iwmax,ndex3d)
38c
39       CALL histwrite(nid_regdyn,"hw4",itau_w,histoW(:,:,:,4),
40     &               kmaxm1*lmaxm1*iwmax,ndex3d)
41c
42       CALL histwrite(nid_regdyn,"nh4",itau_w,nhistoW(:,:,:,4),
43     &               kmaxm1*lmaxm1*iwmax,ndex3d)
44c
45       CALL histwrite(nid_regdyn,"nht4",itau_w,nhistoWt(:,:,:,4),
46     &               kmaxm1*lmaxm1*iwmax,ndex3d)
47c
48       CALL histwrite(nid_regdyn,"hw5",itau_w,histoW(:,:,:,5),
49     &               kmaxm1*lmaxm1*iwmax,ndex3d)
50c
51       CALL histwrite(nid_regdyn,"nh5",itau_w,nhistoW(:,:,:,5),
52     &               kmaxm1*lmaxm1*iwmax,ndex3d)
53c
54       CALL histwrite(nid_regdyn,"nht5",itau_w,nhistoWt(:,:,:,5),
55     &               kmaxm1*lmaxm1*iwmax,ndex3d)
56
57      if (ok_sync) then
58        call histsync(nid_regdyn)
59      endif
60
61      endif ! is_sequential
62
63      endif
Note: See TracBrowser for help on using the repository browser.