source: LMDZ4/trunk/libf/cosp/write_histmthCOSP.h @ 1414

Last change on this file since 1414 was 1279, checked in by Laurent Fairhead, 15 years ago

Merged LMDZ4-dev branch changes r1241:1278 into the trunk
Running trunk and LMDZ4-dev in LMDZOR configuration on local
machine (sequential) and SX8 (4-proc) yields identical results
(restart and restartphy are identical binarily)
Log history from r1241 to r1278 is available by switching to
source:LMDZ4/branches/LMDZ4-dev-20091210

File size: 4.7 KB
Line 
1! Ecriture des fichiers de sorties COSP
2! Sorties journalierres
3! Abderrahmane Idelkadi Septembre 2009
4
5      IF (MOD(itap,NINT(freq_COSP/dtime)).EQ.0) THEN
6
7       itau_wcosp = itau_phy + itap
8
9! Sorties LIDAR
10       if (cfg%Llidar_sim) then
11         if (cfg%Lcllcalipso) then
12          CALL histwrite_phy(nid_mth_cosp,"cllcalipso",itau_wcosp,stlidar%cldlayer(:,1))
13         endif
14         if (cfg%Lclhcalipso) then
15          CALL histwrite_phy(nid_mth_cosp,"clhcalipso",itau_wcosp,stlidar%cldlayer(:,3))
16         endif
17         if (cfg%Lclmcalipso) then
18          CALL histwrite_phy(nid_mth_cosp,"clmcalipso",itau_wcosp,stlidar%cldlayer(:,2))
19         endif
20         if (cfg%Lcltcalipso) then
21          CALL histwrite_phy(nid_mth_cosp,"cltcalipso",itau_wcosp,stlidar%cldlayer(:,4)) 
22         endif
23         if (cfg%Lclcalipso) then
24          CALL histwrite_phy(nid_mth_cosp,"clcalipso",itau_wcosp,stlidar%lidarcld)
25         endif
26         if (cfg%Lcfad_lidarsr532) then
27           do ii=1,SR_BINS
28            CALL histwrite_phy(nid_mth_cosp,"cfad_lidarsr532_"//chcol(ii),itau_wcosp,stlidar%cfad_sr(:,ii,:))
29           enddo
30         endif
31         if (cfg%Lparasol_refl) then
32           CALL histwrite_phy(nid_mth_cosp,"parasol_refl",itau_wcosp,stlidar%parasolrefl)
33         endif
34         if (cfg%Latb532) then
35           do ii=1,Ncolumns
36            CALL histwrite_phy(nid_mth_cosp,"atb532_"//chcol(ii),itau_wcosp,sglidar%beta_tot(:,ii,:))
37           enddo
38         endif
39         if (cfg%Lbeta_mol532) then
40           CALL histwrite_phy(nid_mth_cosp,"beta_mol532",itau_wcosp,sglidar%beta_mol)
41         endif
42        endif ! Lidar
43
44! Sorties RADAR
45!Attention A FAIRE
46!        if (cfg%Lradar_sim) then
47!         print*,'Ecriture sorties Radar'
48!          if (cfg%Lcfad_dbze94) then
49!              print*,'Ecriture de cfad_dbze94.nc '
50!              A revoir l axe vertical Nlvgrid
51!               do ii=1,DBZE_BINS
52!                   dbze_ax(ii) = CFAD_ZE_MIN + CFAD_ZE_WIDTH*(ii - 0.5)
53!               enddo
54!               call write_netcdf4d('cfad_dbze94.nc',use_vgrid,nlon,nlat,Nlevout,DBZE_BINS, &
55!                                   x,y,out_levs,dbze_ax,i,ndays,time,stradar%cfad_ze)
56!          endif
57!          if (cfg%Lclcalipso2) then
58!               call write_netcdf3d('clcalipso2.nc',use_vgrid,'clcalipso2', &
59!                              nlon,nlat,Nlevout,x,y,out_levs,i,ndays,time,stradar%lidar_only_freq_cloud)
60!          endif
61!          if (cfg%Ldbze94) then
62!             do ii=1,Ncolumns
63!                xcol(ii)=float(i)
64!             enddo
65!             call write_netcdf4d('dbze94.nc',use_vgrid,nlon,nlat,Nlevout,Ncolumns, &
66!                                 x,y,out_levs,xcol,i,ndays,time,sgradar%Ze_tot)
67!          endif
68!          if (cfg%Lcltlidarradar) then
69!             call write_netcdf2d('cltlidarradar.nc','cltlidarradar', &
70!                                 nlon,nlat,x,y,i,ndays,time,stradar%radar_lidar_tcc)
71!          endif
72!        endif  ! Radar
73
74! Sorties MISR
75!Attention A FAIRE
76!        if (cfg%Lmisr_sim) then
77!         print*,'Ecriture sorties Misr'
78!            call write_netcdf4d('clMISR.nc',use_vgrid,nlon,nlat,MISR_N_CTH,7, &
79!                                x,y,MISR_CTH,ISCCP_TAU,i,ndays,time,misr%fq_MISR)
80!        endif
81
82! Sorties ISCCP
83        if (cfg%Lisccp_sim) then
84          if (cfg%Lclisccp2) then
85            do ii=1,7
86              CALL histwrite_phy(nid_mth_cosp,"clisccp2_"//chcol(ii),itau_wcosp,isccp%fq_isccp(:,ii,:))
87            enddo
88          endif
89          if (cfg%Lboxtauisccp) then
90             CALL histwrite_phy(nid_mth_cosp,"boxtauisccp",itau_wcosp,isccp%boxtau)
91          endif
92          if (cfg%Lboxptopisccp) then
93             CALL histwrite_phy(nid_mth_cosp,"boxptopisccp",itau_wcosp,isccp%boxptop)
94          endif
95          if (cfg%Ltclisccp) then
96             CALL histwrite_phy(nid_mth_cosp,"tclisccp",itau_wcosp,isccp%totalcldarea)
97          endif
98          if (cfg%Lctpisccp) then
99             CALL histwrite_phy(nid_mth_cosp,"ctpisccp",itau_wcosp,isccp%meanptop)
100          endif
101          if (cfg%Ltauisccp) then
102             CALL histwrite_phy(nid_mth_cosp,"tauisccp",itau_wcosp,isccp%meantaucld)
103          endif
104          if (cfg%Lalbisccp) then
105             CALL histwrite_phy(nid_mth_cosp,"albisccp",itau_wcosp,isccp%meanalbedocld)
106          endif
107          if (cfg%Lmeantbisccp) then
108             CALL histwrite_phy(nid_mth_cosp,"meantbisccp",itau_wcosp,isccp%meantb)
109          endif
110          if (cfg%Lmeantbclrisccp) then
111             CALL histwrite_phy(nid_mth_cosp,"meantbclrisccp",itau_wcosp,isccp%meantbclr)
112          endif
113        endif ! Isccp
114
115!       if (ok_sync) then
116!$OMP MASTER
117        call histsync(nid_mth_cosp)
118!$OMP END MASTER     
119!       endif
120
121      ENDIF ! if freq_COSP
Note: See TracBrowser for help on using the repository browser.