| 1 | ! |
|---|
| 2 | ! $Header: /home/cvsroot/LMDZ4/libf/phylmd/write_histmth.h,v 1.2 2005/03/09 12:30:16 fairhead Exp $ |
|---|
| 3 | ! |
|---|
| 4 | IF (ok_mensuel) THEN |
|---|
| 5 | |
|---|
| 6 | itau_w = itau_phy + itap |
|---|
| 7 | |
|---|
| 8 | c------------------------------------------------------- |
|---|
| 9 | IF(lev_histmth.GE.1) THEN |
|---|
| 10 | |
|---|
| 11 | ccccccccccccc 2D fields, basics |
|---|
| 12 | |
|---|
| 13 | call histwrite_phy(nid_mth,.false.,"phis",itau_w,pphis) |
|---|
| 14 | call histwrite_phy(nid_mth,.false.,"aire",itau_w,airephy) |
|---|
| 15 | call histwrite_phy(nid_mth,.false.,"tsol",itau_w,ftsol) |
|---|
| 16 | call histwrite_phy(nid_mth,.false.,"psol",itau_w,paprs(:,1)) |
|---|
| 17 | c call histwrite_phy(nid_mth,.false.,"ue",itau_w,ue) |
|---|
| 18 | c VENUS: regardee a l'envers!!!!!!!!!!!!!!! |
|---|
| 19 | c call histwrite_phy(nid_mth,.false.,"ve",itau_w,-1.*ve) |
|---|
| 20 | c call histwrite_phy(nid_mth,.false.,"cdragh",itau_w,cdragh) |
|---|
| 21 | c call histwrite_phy(nid_mth,.false.,"cdragm",itau_w,cdragm) |
|---|
| 22 | |
|---|
| 23 | ENDIF !lev_histmth.GE.1 |
|---|
| 24 | |
|---|
| 25 | c------------------------------------------------------- |
|---|
| 26 | IF(lev_histmth.GE.2) THEN |
|---|
| 27 | |
|---|
| 28 | ccccccccccccc 3D fields, basics |
|---|
| 29 | |
|---|
| 30 | call histwrite_phy(nid_mth,.false.,"temp",itau_w,t_seri) |
|---|
| 31 | call histwrite_phy(nid_mth,.false.,"pres",itau_w,pplay) |
|---|
| 32 | call histwrite_phy(nid_mth,.false.,"geop",itau_w,zphi) |
|---|
| 33 | call histwrite_phy(nid_mth,.false.,"vitu",itau_w,u_seri) |
|---|
| 34 | c VENUS: regardee a l'envers!!!!!!!!!!!!!!! |
|---|
| 35 | call histwrite_phy(nid_mth,.false.,"vitv",itau_w,-1.*v_seri) |
|---|
| 36 | call histwrite_phy(nid_mth,.false.,"vitw",itau_w,omega) |
|---|
| 37 | c en (m/s)/s |
|---|
| 38 | call histwrite_phy(nid_mth,.false.,"dudyn",itau_w,d_u_dyn) |
|---|
| 39 | c en (m/s)/s |
|---|
| 40 | call histwrite_phy(nid_mth,.false.,"duvdf",itau_w,d_u_vdf) |
|---|
| 41 | c call histwrite_phy(nid_mth,.false.,"mang",itau_w,mang) |
|---|
| 42 | c call histwrite_phy(nid_mth,.false.,"Kz",itau_w,ycoefh) |
|---|
| 43 | |
|---|
| 44 | c plusieurs traceurs |
|---|
| 45 | IF (iflag_trac.eq.1) THEN |
|---|
| 46 | DO iq=1,nqmax |
|---|
| 47 | call histwrite_phy(nid_mth,.false.,tname(iq),itau_w,qx(:,:,iq)) |
|---|
| 48 | ENDDO |
|---|
| 49 | ENDIF |
|---|
| 50 | |
|---|
| 51 | call histwrite_phy(nid_mth,.false.,"tops",itau_w,topsw) |
|---|
| 52 | |
|---|
| 53 | ENDIF !lev_histmth.GE.2 |
|---|
| 54 | |
|---|
| 55 | c------------------------------------------------------- |
|---|
| 56 | IF(lev_histmth.GE.3) THEN |
|---|
| 57 | |
|---|
| 58 | cccccccccccccccccc Radiative transfer |
|---|
| 59 | |
|---|
| 60 | c 2D |
|---|
| 61 | |
|---|
| 62 | call histwrite_phy(nid_mth,.false.,"topl",itau_w,toplw) |
|---|
| 63 | call histwrite_phy(nid_mth,.false.,"sols",itau_w,solsw) |
|---|
| 64 | call histwrite_phy(nid_mth,.false.,"soll",itau_w,sollw) |
|---|
| 65 | |
|---|
| 66 | c 3D |
|---|
| 67 | |
|---|
| 68 | call histwrite_phy(nid_mth,.false.,"SWnet",itau_w,swnet) |
|---|
| 69 | call histwrite_phy(nid_mth,.false.,"LWnet",itau_w,lwnet) |
|---|
| 70 | c call histwrite_phy(nid_mth,.false.,"fluxvdf",itau_w,fluxt) |
|---|
| 71 | c call histwrite_phy(nid_mth,.false.,"fluxdyn",itau_w,flux_dyn) |
|---|
| 72 | c call histwrite_phy(nid_mth,.false.,"fluxajs",itau_w,flux_ajs) |
|---|
| 73 | c call histwrite_phy(nid_mth,.false.,"fluxec",itau_w,flux_ec) |
|---|
| 74 | |
|---|
| 75 | ENDIF !lev_histmth.GE.3 |
|---|
| 76 | |
|---|
| 77 | c------------------------------------------------------- |
|---|
| 78 | IF(lev_histmth.GE.4) THEN |
|---|
| 79 | |
|---|
| 80 | c en K/s |
|---|
| 81 | call histwrite_phy(nid_mth,.false.,"dtdyn",itau_w,d_t_dyn) |
|---|
| 82 | c en K/s |
|---|
| 83 | c call histwrite_phy(nid_mth,.false.,"dtphy",itau_w,d_t) |
|---|
| 84 | c en K/s |
|---|
| 85 | call histwrite_phy(nid_mth,.false.,"dtvdf",itau_w,d_t_vdf) |
|---|
| 86 | c en K/s |
|---|
| 87 | call histwrite_phy(nid_mth,.false.,"dtajs",itau_w,d_t_ajs) |
|---|
| 88 | c K/day ==> K/s |
|---|
| 89 | call histwrite_phy(nid_mth,.false.,"dtswr",itau_w,heat/RDAY) |
|---|
| 90 | c K/day ==> K/s |
|---|
| 91 | call histwrite_phy(nid_mth,.false.,"dtlwr",itau_w,-1.*cool/RDAY) |
|---|
| 92 | c en K/s |
|---|
| 93 | c call histwrite_phy(nid_mth,.false.,"dtec",itau_w,d_t_ec) |
|---|
| 94 | c en (m/s)/s |
|---|
| 95 | call histwrite_phy(nid_mth,.false.,"duajs",itau_w,d_u_ajs) |
|---|
| 96 | c en (m/s)/s |
|---|
| 97 | call histwrite_phy(nid_mth,.false.,"dugwo",itau_w,d_u_oro) |
|---|
| 98 | c en (m/s)/s |
|---|
| 99 | call histwrite_phy(nid_mth,.false.,"dugwno",itau_w,d_u_hin) |
|---|
| 100 | c en (m/s)/s |
|---|
| 101 | c VENUS: regardee a l'envers!!!!!!!!!!!!!!! |
|---|
| 102 | c call histwrite_phy(nid_mth,.false.,"dvvdf",itau_w,-1.*d_v_vdf) |
|---|
| 103 | |
|---|
| 104 | ENDIF !lev_histmth.GE.4 |
|---|
| 105 | |
|---|
| 106 | c------------------------------------------------------- |
|---|
| 107 | IF(lev_histmth.GE.5) THEN |
|---|
| 108 | |
|---|
| 109 | c call histwrite_phy(nid_mth,.false.,"taux_",itau_w,fluxu) |
|---|
| 110 | c call histwrite_phy(nid_mth,.false.,"tauy_",itau_w,fluxv) |
|---|
| 111 | c call histwrite_phy(nid_mth,.false.,"cdrm",itau_w,cdragm) |
|---|
| 112 | c call histwrite_phy(nid_mth,.false.,"cdrh",itau_w,cdragh) |
|---|
| 113 | |
|---|
| 114 | ENDIF !lev_histmth.GE.5 |
|---|
| 115 | c------------------------------------------------------- |
|---|
| 116 | |
|---|
| 117 | if (ok_sync) then |
|---|
| 118 | call histsync(nid_mth) |
|---|
| 119 | endif |
|---|
| 120 | |
|---|
| 121 | ENDIF |
|---|
