1 |
|
|
MODULE output_physiqex_mod |
2 |
|
|
|
3 |
|
|
|
4 |
|
|
|
5 |
|
|
CONTAINS |
6 |
|
|
|
7 |
|
|
SUBROUTINE output_physiqex(debut,zjulian,pdtphys,presnivs,paprs,u,v,t,qx,cf,zqr,zqs,zqg,ptke,theta) |
8 |
|
|
|
9 |
|
|
USE dimphy, only : klon,klev |
10 |
|
|
USE iophy, only : histbeg_phy,histwrite_phy |
11 |
|
|
USE ioipsl, only : histvert,histdef,histend,histsync |
12 |
|
|
USE mod_phys_lmdz_para, only : jj_nb |
13 |
|
|
USE ioipsl_getin_p_mod, ONLY : getin_p |
14 |
|
|
USE mod_grid_phy_lmdz, ONLY: nbp_lon,nbp_lat |
15 |
|
|
USE iophy, ONLY : init_iophy_new |
16 |
|
|
USE geometry_mod, ONLY: latitude_deg, longitude_deg |
17 |
|
|
USE infotrac_phy, only : nqtot |
18 |
|
|
|
19 |
|
|
|
20 |
|
|
|
21 |
|
|
implicit none |
22 |
|
|
logical, intent(in) :: debut |
23 |
|
|
real, intent(in) :: pdtphys,zjulian |
24 |
|
|
real,intent(in) :: presnivs(klev) ! pseudo-pressure (Pa) of mid-layers |
25 |
|
|
real,intent(in) :: u(klon,klev) ! eastward zonal wind (m/s) |
26 |
|
|
real,intent(in) :: v(klon,klev) ! northward meridional wind (m/s) |
27 |
|
|
real,intent(in) :: t(klon,klev) ! temperature (K) |
28 |
|
|
real,intent(in) :: theta(klon,klev) ! temperature (K) |
29 |
|
|
real,intent(in) :: paprs(klon,klev+1) ! interlayer pressure (Pa) |
30 |
|
|
real,intent(in) :: qx(klon,klev,nqtot) !tracers |
31 |
|
|
real,intent(in) :: cf(klon,klev) !cloud fraction |
32 |
|
|
real,intent(in) :: zqr(klon,klev) !rain specifiq content |
33 |
|
|
real,intent(in) :: zqs(klon,klev) !snow specifiq content |
34 |
|
|
real,intent(in) :: zqg(klon,klev) !graupel specifiq content |
35 |
|
|
real,intent(in) :: ptke(klon,klev) !tke |
36 |
|
|
|
37 |
|
|
real :: t_ops ! frequency of the IOIPSL operations (eg average over...) |
38 |
|
|
real :: t_wrt ! frequency of the IOIPSL outputs |
39 |
|
|
integer :: zvertid ! vertical coordinate ID |
40 |
|
|
real :: dtime |
41 |
|
|
|
42 |
|
|
integer,save :: iwrite_phys=1 ! output every iwrite_phys physics step |
43 |
|
|
!$OMP THREADPRIVATE(iwrite_phys) |
44 |
|
|
integer :: nhori ! horizontal coordinate ID |
45 |
|
|
integer,save :: nid_hist ! output file ID |
46 |
|
|
!$OMP THREADPRIVATE(nid_hist) |
47 |
|
|
integer, save :: itau=0 |
48 |
|
|
!$OMP THREADPRIVATE(itau) |
49 |
|
|
|
50 |
|
|
integer, save :: ioex=1 |
51 |
|
|
|
52 |
|
|
|
53 |
|
|
print*,'nnnnnnn ',nid_hist,debut,itau |
54 |
|
|
|
55 |
|
|
if(debut)then |
56 |
|
|
|
57 |
|
|
call getin_p("iwrite_phys",iwrite_phys) |
58 |
|
|
|
59 |
|
|
!$OMP MASTER |
60 |
|
|
CALL iophys_ini(pdtphys) |
61 |
|
|
!$OMP END MASTER |
62 |
|
|
!$OMP BARRIER |
63 |
|
|
|
64 |
|
|
endif |
65 |
|
|
|
66 |
|
|
|
67 |
|
|
itau=itau+1 |
68 |
|
|
|
69 |
|
|
if (modulo(itau,iwrite_phys)==0) then |
70 |
|
|
call iophys_ecrit('temp',klev,'Temperature','K',t) |
71 |
|
|
call iophys_ecrit('u',klev,'zonal wind','m/s',u) |
72 |
|
|
call iophys_ecrit('v',klev,'meridinal wind','m/s',v) |
73 |
|
|
call iophys_ecrit('ps',1,'Surface pressure','Pa',paprs(:,1)) |
74 |
|
|
call iophys_ecrit('qv',klev,'Water vapor specifiq content', 'kg/kg', qx(:,:,1)) |
75 |
|
|
call iophys_ecrit('qc',klev,'Cloud liquid water specifiq content', 'kg/kg', qx(:,:,2)) |
76 |
|
|
call iophys_ecrit('qi',klev,'Cloud solid water specifiq content', 'kg/kg', qx(:,:,3)) |
77 |
|
|
call iophys_ecrit('CF',klev,'Cloud fraction', '0-1', cf) |
78 |
|
|
call iophys_ecrit('qr',klev,'Rain specifiq content', 'kg/kg', zqr) |
79 |
|
|
call iophys_ecrit('qs',klev,'Snow specifiq content', 'kg/kg', zqs) |
80 |
|
|
call iophys_ecrit('qg',klev,'Graupel specifiq content', 'kg/kg', zqg) |
81 |
|
|
call iophys_ecrit('TKE',klev,'TKE', 'm2/s2', ptke) |
82 |
|
|
call iophys_ecrit('theta',klev,'Temperature potentielle', 'K', theta) |
83 |
|
|
endif |
84 |
|
|
|
85 |
|
|
|
86 |
|
|
END SUBROUTINE output_physiqex |
87 |
|
|
END MODULE output_physiqex_mod |