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