Changeset 1862 for LMDZ5/trunk/libf/phylmd/ini_histday_seri.h
- Timestamp:
- Sep 10, 2013, 2:18:14 PM (11 years ago)
- File:
-
- 1 edited
-
LMDZ5/trunk/libf/phylmd/ini_histday_seri.h (modified) (2 diffs)
Legend:
- Unmodified
- Added
- Removed
-
LMDZ5/trunk/libf/phylmd/ini_histday_seri.h
r1403 r1862 1 c 2 c$Id$3 c 4 cym Ne fonctionnera pas en mode parallele1 ! 2 ! $Id$ 3 ! 4 !ym Ne fonctionnera pas en mode parallele 5 5 IF (is_sequential) THEN 6 6 7 7 IF (type_run.EQ."AMIP") THEN 8 c 8 ! 9 9 zstophy = dtime 10 10 zout = ecrit_day 11 c 11 ! 12 12 idayref = day_ref 13 13 CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian) 14 c 14 ! 15 15 CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon) 16 16 DO i = 1, iim … … 22 22 ENDDO 23 23 CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat) 24 c 24 ! 25 25 imin_debut=1 26 26 nbpti=1 27 27 jmin_debut=1 28 28 nbptj=1 29 c 30 CALL histbeg("histday_seri.nc", 31 . iim,zx_lon(:,1), jjmp1,zx_lat(1,:),32 . imin_debut,nbpti,jmin_debut,nbptj,33 . itau_phy, zjulian, dtime,34 .nhori, nid_day_seri)35 c 36 CALL histvert(nid_day_seri, "presnivs", 37 . "Vertical levels","mb",38 .klev, presnivs/100., nvert)39 c 40 CALL histdef(nid_day_seri, "bilTOA", 41 . "Net radiation at model top", "W/m2",42 . iim,jjmp1,nhori, 1,1,1, -99, 32,43 ."ave(X)", zstophy,zout)44 c 45 CALL histdef(nid_day_seri, "bils", 46 . "Net downward energy flux at surface","W/m2",47 . iim,jjmp1,nhori, 1,1,1, -99, 32,48 ."ave(X)", zstophy,zout)49 c 50 CALL histdef(nid_day_seri, "ecin", 51 . "Total kinetic energy (per unit area)","J/m2",52 . iim,jjmp1,nhori, 1,1,1, -99, 32,53 ."ave(X)", zstophy,zout)54 c 55 cIM 151004 BEG29 ! 30 CALL histbeg("histday_seri.nc", & 31 iim,zx_lon(:,1), jjmp1,zx_lat(1,:), & 32 imin_debut,nbpti,jmin_debut,nbptj, & 33 itau_phy, zjulian, dtime, & 34 nhori, nid_day_seri) 35 ! 36 CALL histvert(nid_day_seri, "presnivs", & 37 "Vertical levels","mb", & 38 klev, presnivs/100., nvert) 39 ! 40 CALL histdef(nid_day_seri, "bilTOA", & 41 "Net radiation at model top", "W/m2", & 42 iim,jjmp1,nhori, 1,1,1, -99, 32, & 43 "ave(X)", zstophy,zout) 44 ! 45 CALL histdef(nid_day_seri, "bils", & 46 "Net downward energy flux at surface","W/m2", & 47 iim,jjmp1,nhori, 1,1,1, -99, 32, & 48 "ave(X)", zstophy,zout) 49 ! 50 CALL histdef(nid_day_seri, "ecin", & 51 "Total kinetic energy (per unit area)","J/m2", & 52 iim,jjmp1,nhori, 1,1,1, -99, 32, & 53 "ave(X)", zstophy,zout) 54 ! 55 !IM 151004 BEG 56 56 IF(1.EQ.0) THEN 57 c 58 CALL histdef(nid_day_seri, "momang", 59 . "Total relative angular momentum (per unit area)",60 . "kg/s",61 . iim,jjmp1,nhori, 1,1,1, -99, 32,62 ."ave(X)", zstophy,zout)63 c 64 CALL histdef(nid_day_seri, "frictor", 65 . "Friction torque (per unit area)", "N/m",66 . iim,jjmp1,nhori, 1,1,1, -99, 32,67 ."ave(X)", zstophy,zout)68 c 69 CALL histdef(nid_day_seri, "mountor", 70 . "Mountain torque (per unit area)", "N/m",71 . iim,jjmp1,nhori, 1,1,1, -99, 32,72 ."ave(X)", zstophy,zout)73 c 57 ! 58 CALL histdef(nid_day_seri, "momang", & 59 "Total relative angular momentum (per unit area)", & 60 "kg/s", & 61 iim,jjmp1,nhori, 1,1,1, -99, 32, & 62 "ave(X)", zstophy,zout) 63 ! 64 CALL histdef(nid_day_seri, "frictor", & 65 "Friction torque (per unit area)", "N/m", & 66 iim,jjmp1,nhori, 1,1,1, -99, 32, & 67 "ave(X)", zstophy,zout) 68 ! 69 CALL histdef(nid_day_seri, "mountor", & 70 "Mountain torque (per unit area)", "N/m", & 71 iim,jjmp1,nhori, 1,1,1, -99, 32, & 72 "ave(X)", zstophy,zout) 73 ! 74 74 ENDIF !(1.EQ.0) THEN 75 c 76 CALL histdef(nid_day_seri, "momang", 77 . "Axial angular momentum (per unit area)",78 . "kg/s",79 . iim,jjmp1,nhori, 1,1,1, -99, 32,80 ."ave(X)", zstophy,zout)81 c 82 CALL histdef(nid_day_seri, "torsfc", 83 . "Total surface torque (including mountain torque)", "N/m",84 . iim,jjmp1,nhori, 1,1,1, -99, 32,85 ."ave(X)", zstophy,zout)86 c 87 cIM 151004 END88 c 89 CALL histdef(nid_day_seri, "tamv", 90 . "Temperature (mass-weighted vert. ave)", "K",91 . iim,jjmp1,nhori, 1,1,1, -99, 32,92 ."ave(X)", zstophy,zout)93 c 94 CALL histdef(nid_day_seri, "psol", 95 . "Surface pressure", "Pa",96 . iim,jjmp1,nhori, 1,1,1, -99, 32,97 ."ave(X)", zstophy,zout)98 c 99 CALL histdef(nid_day_seri, "evap", 100 . "Evaporation and sublimation (per unit area)",101 . "kg/(m2*s)",102 . iim,jjmp1,nhori, 1,1,1, -99, 32,103 ."ave(X)", zstophy,zout)104 c 105 ccall histdef(nid_day_seri,106 c. "SnowFrac",107 c. "Snow-covered area ", "%",108 c. iim,jjmp1,nhori, 1,1,1, -99, 32,109 c. "ave(X)", zstophy,zout)110 c 111 cCALL histdef(nid_day_seri, "snow_depth",112 cIM 080904 . "Snow Depth (water equivalent)", "m",113 cIM 191104 . "Snow Depth (water equivalent)", "kg/m2",114 c. "Snow Mass", "kg/m2",115 c. iim,jjmp1,nhori, 1,1,1, -99, 32,116 c. "ave(X)", zstophy,zout)117 c 118 call histdef(nid_day_seri, 119 . "tsol_"//clnsurf(is_oce),120 . "SST over open (ice-free) ocean ", "K",121 . iim,jjmp1,nhori, 1,1,1, -99, 32,122 ."ave(X)", zstophy,zout)123 c 124 c=================================================================125 c 75 ! 76 CALL histdef(nid_day_seri, "momang", & 77 "Axial angular momentum (per unit area)", & 78 "kg/s", & 79 iim,jjmp1,nhori, 1,1,1, -99, 32, & 80 "ave(X)", zstophy,zout) 81 ! 82 CALL histdef(nid_day_seri, "torsfc", & 83 "Total surface torque (including mountain torque)", "N/m", & 84 iim,jjmp1,nhori, 1,1,1, -99, 32, & 85 "ave(X)", zstophy,zout) 86 ! 87 !IM 151004 END 88 ! 89 CALL histdef(nid_day_seri, "tamv", & 90 "Temperature (mass-weighted vert. ave)", "K", & 91 iim,jjmp1,nhori, 1,1,1, -99, 32, & 92 "ave(X)", zstophy,zout) 93 ! 94 CALL histdef(nid_day_seri, "psol", & 95 "Surface pressure", "Pa", & 96 iim,jjmp1,nhori, 1,1,1, -99, 32, & 97 "ave(X)", zstophy,zout) 98 ! 99 CALL histdef(nid_day_seri, "evap", & 100 "Evaporation and sublimation (per unit area)", & 101 "kg/(m2*s)", & 102 iim,jjmp1,nhori, 1,1,1, -99, 32, & 103 "ave(X)", zstophy,zout) 104 ! 105 ! call histdef(nid_day_seri, 106 ! . "SnowFrac", 107 ! . "Snow-covered area ", "%", 108 ! . iim,jjmp1,nhori, 1,1,1, -99, 32, 109 ! . "ave(X)", zstophy,zout) 110 ! 111 ! CALL histdef(nid_day_seri, "snow_depth", 112 !IM 080904 . "Snow Depth (water equivalent)", "m", 113 !IM 191104 . "Snow Depth (water equivalent)", "kg/m2", 114 ! . "Snow Mass", "kg/m2", 115 ! . iim,jjmp1,nhori, 1,1,1, -99, 32, 116 ! . "ave(X)", zstophy,zout) 117 ! 118 call histdef(nid_day_seri, & 119 "tsol_"//clnsurf(is_oce), & 120 "SST over open (ice-free) ocean ", "K", & 121 iim,jjmp1,nhori, 1,1,1, -99, 32, & 122 "ave(X)", zstophy,zout) 123 ! 124 !================================================================= 125 ! 126 126 CALL histend(nid_day_seri) 127 c 128 c=================================================================127 ! 128 !================================================================= 129 129 ENDIF ! fin de test sur type_run.EQ.AMIP 130 130
Note: See TracChangeset
for help on using the changeset viewer.
