source: LMDZ4/branches/LMDZ4_AR5/libf/cosp/write_histdayCOSP.h @ 3536

Last change on this file since 3536 was 1928, checked in by idelkadi, 11 years ago
  • Correction du bug dans l'interface avec la physique : On ne veut pas que la distinction entre les nuages convectifs et stratiformes soit prise en compte dans les calculs Cosp. Dans l'interface avec la physique, sont passees en entree pour Cosp, les quantites totales (stratiforme + convective) du contenus en eau et d'autres variables. La fraction nuageuse convective calculee dans la physique est passee en entree pour Cosp dans la version buggee. Elle est remise a 0 dans cette version corrigee.
  • Mise a jour pour ISCCP : la fraction d'ensoleillement calculee par LMDZ est passee en argument pour Isccp
  • Rajouts d'autre variables diagnostiques : sunlit (=0 si nuit et =1 si jour) parasol_crefl (reflectance integree)
File size: 5.2 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_day_cosp,lCOSP,"cllcalipso", &
13     itau_wcosp,stlidar%cldlayer(:,1))
14         endif
15         if (cfg%Lclhcalipso) then
16          CALL histwrite_phy(nid_day_cosp,lCOSP,"clhcalipso", &
17     itau_wcosp,stlidar%cldlayer(:,3))
18         endif
19         if (cfg%Lclmcalipso) then
20          CALL histwrite_phy(nid_day_cosp,lCOSP,"clmcalipso", &
21     itau_wcosp,stlidar%cldlayer(:,2))
22         endif
23         if (cfg%Lcltcalipso) then
24          CALL histwrite_phy(nid_day_cosp,lCOSP,"cltcalipso", &
25     itau_wcosp,stlidar%cldlayer(:,4)) 
26         endif
27         if (cfg%Lclcalipso) then
28          CALL histwrite_phy(nid_day_cosp,lCOSP,"clcalipso", &
29     itau_wcosp,stlidar%lidarcld)
30         endif
31         if (cfg%Lcfad_lidarsr532) then
32           do ii=1,SR_BINS
33            CALL histwrite_phy(nid_day_cosp,lCOSP, &
34     "cfad_lidarsr532_"//chcol(ii),itau_wcosp,stlidar%cfad_sr(:,ii,:))
35           enddo
36         endif
37         if (cfg%Lparasol_refl) then
38           CALL histwrite_phy(nid_day_cosp,lCOSP,"parasol_refl", &
39     itau_wcosp,stlidar%parasolrefl)
40           CALL histwrite_phy(nid_day_cosp,lCOSP,"parasol_crefl", &
41     itau_wcosp,parasolcrefl)
42           CALL histwrite_phy(nid_day_cosp,lCOSP,"Ncrefl", &
43     itau_wcosp,Ncref)
44         endif
45         if (cfg%Latb532) then
46           do ii=1,Ncolumns
47            CALL histwrite_phy(nid_day_cosp,lCOSP,"atb532_"//chcol(ii), &
48     itau_wcosp,sglidar%beta_tot(:,ii,:))
49           enddo
50         endif
51         if (cfg%Lbeta_mol532) then
52           CALL histwrite_phy(nid_day_cosp,lCOSP,"beta_mol532", &
53     itau_wcosp,sglidar%beta_mol)
54         endif
55        endif ! Lidar
56
57! Sorties RADAR
58!Attention A FAIRE
59!        if (cfg%Lradar_sim) then
60!         print*,'Ecriture sorties Radar'
61!          if (cfg%Lcfad_dbze94) then
62!              print*,'Ecriture de cfad_dbze94.nc '
63!              A revoir l axe vertical Nlvgrid
64!               do ii=1,DBZE_BINS
65!                   dbze_ax(ii) = CFAD_ZE_MIN + CFAD_ZE_WIDTH*(ii - 0.5)
66!               enddo
67!               call write_netcdf4d('cfad_dbze94.nc',use_vgrid,nlon,nlat,Nlevout,DBZE_BINS, &
68!                                   x,y,out_levs,dbze_ax,i,ndays,time,stradar%cfad_ze)
69!          endif
70!          if (cfg%Lclcalipso2) then
71!               call write_netcdf3d('clcalipso2.nc',use_vgrid,'clcalipso2', &
72!                              nlon,nlat,Nlevout,x,y,out_levs,i,ndays,time,stradar%lidar_only_freq_cloud)
73!          endif
74!          if (cfg%Ldbze94) then
75!             do ii=1,Ncolumns
76!                xcol(ii)=float(i)
77!             enddo
78!             call write_netcdf4d('dbze94.nc',use_vgrid,nlon,nlat,Nlevout,Ncolumns, &
79!                                 x,y,out_levs,xcol,i,ndays,time,sgradar%Ze_tot)
80!          endif
81!          if (cfg%Lcltlidarradar) then
82!             call write_netcdf2d('cltlidarradar.nc','cltlidarradar', &
83!                                 nlon,nlat,x,y,i,ndays,time,stradar%radar_lidar_tcc)
84!          endif
85!        endif  ! Radar
86
87! Sorties MISR
88!Attention A FAIRE
89!        if (cfg%Lmisr_sim) then
90!         print*,'Ecriture sorties Misr'
91!            call write_netcdf4d('clMISR.nc',use_vgrid,nlon,nlat,MISR_N_CTH,7, &
92!                                x,y,MISR_CTH,ISCCP_TAU,i,ndays,time,misr%fq_MISR)
93!        endif
94
95! Sorties ISCCP
96        if (cfg%Lisccp_sim) then
97          CALL histwrite_phy(nid_day_cosp,lCOSP,"sunlit", &
98     itau_wcosp,sunlit)
99          if (cfg%Lclisccp2) then
100            do ii=1,7
101              CALL histwrite_phy(nid_day_cosp,lCOSP, &
102     "clisccp2_"//chcol(ii),itau_wcosp,isccp%fq_isccp(:,ii,:))
103            enddo
104          endif
105          if (cfg%Lboxtauisccp) then
106             CALL histwrite_phy(nid_day_cosp,lCOSP,"boxtauisccp", &
107     itau_wcosp,isccp%boxtau)
108          endif
109          if (cfg%Lboxptopisccp) then
110             CALL histwrite_phy(nid_day_cosp,lCOSP,"boxptopisccp", &
111     itau_wcosp,isccp%boxptop)
112          endif
113          if (cfg%Ltclisccp) then
114             CALL histwrite_phy(nid_day_cosp,lCOSP,"tclisccp", &
115     itau_wcosp,isccp%totalcldarea)
116          endif
117          if (cfg%Lctpisccp) then
118             CALL histwrite_phy(nid_day_cosp,lCOSP,"ctpisccp", &
119     itau_wcosp,isccp%meanptop)
120          endif
121          if (cfg%Ltauisccp) then
122             CALL histwrite_phy(nid_day_cosp,lCOSP,"tauisccp", &
123     itau_wcosp,isccp%meantaucld)
124          endif
125          if (cfg%Lalbisccp) then
126             CALL histwrite_phy(nid_day_cosp,lCOSP,"albisccp", &
127     itau_wcosp,isccp%meanalbedocld)
128          endif
129          if (cfg%Lmeantbisccp) then
130             CALL histwrite_phy(nid_day_cosp,lCOSP,"meantbisccp", &
131     itau_wcosp,isccp%meantb)
132          endif
133          if (cfg%Lmeantbclrisccp) then
134             CALL histwrite_phy(nid_day_cosp,lCOSP,"meantbclrisccp", &
135     itau_wcosp,isccp%meantbclr)
136          endif
137        endif ! Isccp
138
139!       if (ok_sync) then
140!$OMP MASTER
141        call histsync(nid_day_cosp)
142!$OMP END MASTER     
143!       endif
144
145      ENDIF ! if freq_COSP
Note: See TracBrowser for help on using the repository browser.