Changeset 902
- Timestamp:
- Mar 13, 2013, 11:06:21 AM (12 years ago)
- Location:
- trunk
- Files:
-
- 9 edited
-
LMDZ.TITAN/libf/phytitan/ini_histday.h (modified) (1 diff)
-
LMDZ.TITAN/libf/phytitan/ini_histins.h (modified) (2 diffs)
-
LMDZ.TITAN/libf/phytitan/ini_histmth.h (modified) (1 diff)
-
LMDZ.VENUS/libf/phyvenus/ini_histday.h (modified) (9 diffs)
-
LMDZ.VENUS/libf/phyvenus/ini_histins.h (modified) (4 diffs)
-
LMDZ.VENUS/libf/phyvenus/ini_histmth.h (modified) (6 diffs)
-
LMDZ.VENUS/libf/phyvenus/write_histday.h (modified) (5 diffs)
-
LMDZ.VENUS/libf/phyvenus/write_histins.h (modified) (2 diffs)
-
LMDZ.VENUS/libf/phyvenus/write_histmth.h (modified) (3 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/LMDZ.TITAN/libf/phytitan/ini_histday.h
r888 r902 16 16 CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlatd,zx_lat) 17 17 CALL histbeg("histday", iim,zx_lon(:,1), jjmp1,zx_lat(1,:), 18 . 1,iim,1,jjmp1, itau_phy, zjulian, zsto,18 . 1,iim,1,jjmp1, itau_phy, zjulian, dtime, 19 19 . nhori, nid_day) 20 20 write(*,*)'Journee ', itau_phy, zjulian -
trunk/LMDZ.TITAN/libf/phytitan/ini_histins.h
r888 r902 1 1 IF (ok_instan) THEN 2 2 c 3 zsto1= dtime4 3 zsto = dtime * FLOAT(ecrit_ins) 5 4 zout = dtime * FLOAT(ecrit_ins) … … 15 14 CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlatd,zx_lat) 16 15 CALL histbeg("histins", iim,zx_lon(:,1), jjmp1,zx_lat(1,:), 17 . 1,iim,1,jjmp1, itau_phy, zjulian, zsto1,16 . 1,iim,1,jjmp1, itau_phy, zjulian, dtime, 18 17 . nhori, nid_ins) 19 18 write(*,*)'Inst ', itau_phy, zjulian -
trunk/LMDZ.TITAN/libf/phytitan/ini_histmth.h
r888 r902 16 16 CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlatd,zx_lat) 17 17 CALL histbeg("histmth", iim,zx_lon(:,1), jjmp1,zx_lat(1,:), 18 . 1,iim,1,jjmp1, itau_phy, zjulian, zsto,18 . 1,iim,1,jjmp1, itau_phy, zjulian, dtime, 19 19 . nhori, nid_mth) 20 20 write(*,*)'Journee ', itau_phy, zjulian -
trunk/LMDZ.VENUS/libf/phyvenus/ini_histday.h
r892 r902 6 6 zsto = dtime 7 7 zout = dtime * REAL(ecrit_day) 8 zsto1= dtime * REAL(ecrit_day)9 8 10 9 idayref = day_ref … … 12 11 13 12 call histbeg_phy("histday.nc",itau_phy, 14 . zjulian, zsto,nhori,nid_day)13 . zjulian,dtime,nhori,nid_day) 15 14 16 15 !$OMP MASTER … … 39 38 . "ave(X)", zsto,zout) 40 39 c 41 CALL histdef(nid_day, "ue", "Zonal energy transport", "-", 42 . iim,jj_nb,nhori, 1,1,1, nvert, 32, 43 . "ave(X)", zsto,zout) 44 c 45 CALL histdef(nid_day, "ve", "Merid energy transport", "-", 46 . iim,jj_nb,nhori, 1,1,1, nvert, 32, 47 . "ave(X)", zsto,zout) 40 c CALL histdef(nid_day, "ue", "Zonal energy transport", "-", 41 c . iim,jj_nb,nhori, 1,1,1, nvert, 32, 42 c . "ave(X)", zsto,zout) 43 c 44 c CALL histdef(nid_day, "ve", "Merid energy transport", "-", 45 c . iim,jj_nb,nhori, 1,1,1, nvert, 32, 46 c . "ave(X)", zsto,zout) 47 c 48 c CALL histdef(nid_day, "cdragh", "Drag coef on T", "-", 49 c . iim,jj_nb,nhori, 1,1,1, nvert, 32, 50 c . "ave(X)", zsto,zout) 51 c 52 c CALL histdef(nid_day, "cdragm", "Drag coef on U", "-", 53 c . iim,jj_nb,nhori, 1,1,1, nvert, 32, 54 c . "ave(X)", zsto,zout) 48 55 c 49 56 ENDIF !lev_histday.GE.1 … … 78 85 . "ave(X)", zsto,zout) 79 86 c 87 CALL histdef(nid_day, "dudyn", "Dynamics dU", "m/s2", 88 . iim,jj_nb,nhori, klev,1,klev,nvert, 32, 89 . "ave(X)", zsto,zout) 90 c 91 CALL histdef(nid_day, "duvdf", "Boundary-layer dU", "m/s2", 92 . iim,jj_nb,nhori, klev,1,klev,nvert, 32, 93 . "ave(X)", zsto,zout) 94 c 95 c CALL histdef(nid_day, "mang", "Angular momentum", "kg m2/s", 96 c . iim,jj_nb,nhori, klev,1,klev,nvert, 32, 97 c . "ave(X)", zsto,zout) 98 c 99 c CALL histdef(nid_day, "Kz", "vertical diffusion coef", "m2/s", 100 c . iim,jj_nb,nhori, klev,1,klev,nvert, 32, 101 c . "ave(X)", zsto,zout) 102 c 103 c plusieurs traceurs 104 if (iflag_trac.eq.1) THEN 105 DO iq=1,nqmax 106 IF (iq.LE.99) THEN 107 WRITE(str2,'(i2.2)') iq 108 CALL histdef(nid_day, tname(iq), ttext(iq), "ppm", 109 . iim,jj_nb,nhori, klev,1,klev,nvert, 32, 110 . "ave(X)", zsto,zout) 111 ELSE 112 PRINT*, "Trop de traceurs" 113 CALL abort 114 ENDIF 115 ENDDO 116 endif 117 c 118 CALL histdef(nid_day, "tops", "Solar rad. at TOA", "W/m2", 119 . iim,jj_nb,nhori, 1,1,1, nvert, 32, 120 . "ave(X)", zsto,zout) 121 c 80 122 ENDIF !lev_histday.GE.2 81 123 c … … 87 129 c 2D 88 130 c 89 CALL histdef(nid_day, "tops", "Solar rad. at TOA", "W/m2",90 . iim,jj_nb,nhori, 1,1,1, nvert, 32,91 . "ave(X)", zsto1,zout)92 c93 131 CALL histdef(nid_day, "topl", "IR rad. at TOA", "W/m2", 94 132 . iim,jj_nb,nhori, 1,1,1, nvert, 32, 95 . "ave(X)", zsto 1,zout)133 . "ave(X)", zsto,zout) 96 134 c 97 135 CALL histdef(nid_day, "sols", "Solar rad. at surf.", "W/m2", 98 136 . iim,jj_nb,nhori, 1,1,1, nvert, 32, 99 . "ave(X)", zsto 1,zout)137 . "ave(X)", zsto,zout) 100 138 c 101 139 CALL histdef(nid_day, "soll", "IR rad. at surface", "W/m2", 102 140 . iim,jj_nb,nhori, 1,1,1, nvert, 32, 103 . "ave(X)", zsto 1,zout)141 . "ave(X)", zsto,zout) 104 142 c 105 143 c 3D … … 107 145 CALL histdef(nid_day, "SWnet", "Net SW flux","W/m2", 108 146 . iim,jj_nb,nhori, klev,1,klev,nvert, 109 . 32, "ave(X)", zsto 1,zout)147 . 32, "ave(X)", zsto,zout) 110 148 c 111 149 CALL histdef(nid_day, "LWnet", "Net LW flux","W/m2", 112 150 . iim,jj_nb,nhori, klev,1,klev,nvert, 113 . 32, "ave(X)", zsto1,zout) 151 . 32, "ave(X)", zsto,zout) 152 c 153 CALL histdef(nid_day, "fluxvdf", "PBL net flux","W/m2", 154 . iim,jj_nb,nhori, klev,1,klev,nvert, 155 . 32, "ave(X)", zsto,zout) 156 c 157 CALL histdef(nid_day, "fluxdyn", "Dyn. net flux","W/m2", 158 . iim,jj_nb,nhori, klev,1,klev,nvert, 159 . 32, "ave(X)", zsto,zout) 160 c 161 CALL histdef(nid_day, "fluxajs", "Dry adj. net flux","W/m2", 162 . iim,jj_nb,nhori, klev,1,klev,nvert, 163 . 32, "ave(X)", zsto,zout) 164 c 165 c CALL histdef(nid_day, "fluxec", "Cin. net flux","W/m2", 166 c . iim,jj_nb,nhori, klev,1,klev,nvert, 167 c . 32, "ave(X)", zsto,zout) 114 168 c 115 169 ENDIF !lev_histday.GE.3 … … 122 176 . "ave(X)", zsto,zout) 123 177 c 124 c 125 c 126 c 178 c CALL histdef(nid_day, "dtphy", "Physics dT", "K/s", 179 c . iim,jj_nb,nhori, klev,1,klev,nvert, 32, 180 c . "ave(X)", zsto,zout) 127 181 c 128 182 CALL histdef(nid_day, "dtvdf", "Boundary-layer dT", "K/s", … … 142 196 . "ave(X)", zsto,zout) 143 197 c 144 c CALL histdef(nid_day, "dtec", "Cinetic dissip dT", "K/s", 145 c . iim,jj_nb,nhori, klev,1,klev,nvert, 32, 146 c . "ave(X)", zsto,zout) 147 c 148 CALL histdef(nid_mth, "dudyn", "Dynamics dU", "m/s2", 149 . iim,jj_nb,nhori, klev,1,klev,nvert, 32, 150 . "ave(X)", zsto,zout) 151 c 152 CALL histdef(nid_mth, "duvdf", "Boundary-layer dU", "m/s2", 153 . iim,jj_nb,nhori, klev,1,klev,nvert, 32, 154 . "ave(X)", zsto,zout) 198 c CALL histdef(nid_day, "dtec", "Cinetic dissip dT", "K/s", 199 c . iim,jj_nb,nhori, klev,1,klev,nvert, 32, 200 c . "ave(X)", zsto,zout) 155 201 c 156 202 CALL histdef(nid_mth, "duajs", "Dry convection dU", "m/s2", … … 174 220 c------------------------------------------------------- 175 221 IF(lev_histday.GE.5) THEN 176 c177 if (nqmax.GE.3) THEN178 DO iq=1,nqmax-2179 IF (iq.LE.99) THEN180 WRITE(str2,'(i2.2)') iq181 CALL histdef(nid_day, "trac"//str2, "Tracer No."//str2, "-",182 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,183 . "ave(X)", zsto,zout)184 ELSE185 PRINT*, "Trop de traceurs"186 CALL abort187 ENDIF188 ENDDO189 endif190 222 c 191 223 c call histdef(nid_day, "taux", -
trunk/LMDZ.VENUS/libf/phyvenus/ini_histins.h
r892 r902 4 4 IF (ok_instan) THEN 5 5 6 zsto1= dtime 7 zsto = dtime * ecrit_ins 8 zout = dtime * ecrit_ins 6 zsto = dtime * ecrit_ins 7 zout = dtime * ecrit_ins 9 8 10 9 idayref = day_ref … … 12 11 13 12 call histbeg_phy("histins.nc",itau_phy, 14 . zjulian, zsto,nhori,nid_ins)13 . zjulian,dtime,nhori,nid_ins) 15 14 16 15 !$OMP MASTER … … 47 46 c . "ins(X)", zsto,zout) 48 47 c 49 50 51 52 c 53 54 55 48 c CALL histdef(nid_ins, "cdragh", "Drag coef on T", "-", 49 c . iim,jj_nb,nhori, 1,1,1, nvert, 32, 50 c . "ins(X)", zsto,zout) 51 c 52 c CALL histdef(nid_ins, "cdragm", "Drag coef on U", "-", 53 c . iim,jj_nb,nhori, 1,1,1, nvert, 32, 54 c . "ins(X)", zsto,zout) 56 55 c 57 56 ENDIF !lev_histday.GE.1 … … 98 97 c . "ins(X)", zsto,zout) 99 98 c 100 101 102 99 c CALL histdef(nid_ins, "Kz", "vertical diffusion coef", "m2/s", 100 c . iim,jj_nb,nhori, klev,1,klev,nvert, 32, 101 c . "ins(X)", zsto,zout) 103 102 c 104 103 c plusieurs traceurs -
trunk/LMDZ.VENUS/libf/phyvenus/ini_histmth.h
r892 r902 7 7 zout = dtime * ecrit_mth 8 8 9 c zsto1: pour des valeurs "instantanees" mensuelles10 zsto1 = dtime * ecrit_mth11 c zsto2: pour des flux radiatifs calcules avec un autre pas de temps12 zsto2 = dtime * radpas13 PRINT*,' zsto,zsto1,zsto2,zout=',zsto, zsto1, zsto2,zout14 15 9 idayref = day_ref 16 10 CALL ymds2ju(annee_ref, 1, idayref, zero, zjulian) 17 11 18 12 call histbeg_phy("histmth.nc",itau_phy, 19 . zjulian, zsto,nhori,nid_mth)13 . zjulian,dtime,nhori,nid_mth) 20 14 21 15 !$OMP MASTER … … 52 46 c . "ave(X)", zsto,zout) 53 47 c 54 55 56 57 c 58 59 60 48 c CALL histdef(nid_mth, "cdragh", "Drag coef on T", "-", 49 c . iim,jj_nb,nhori, 1,1,1, nvert, 32, 50 c . "ave(X)", zsto,zout) 51 c 52 c CALL histdef(nid_mth, "cdragm", "Drag coef on U", "-", 53 c . iim,jj_nb,nhori, 1,1,1, nvert, 32, 54 c . "ave(X)", zsto,zout) 61 55 c 62 56 ENDIF !lev_histmth.GE.1 … … 103 97 c . "ave(X)", zsto,zout) 104 98 c 105 106 107 99 c CALL histdef(nid_mth, "Kz", "vertical diffusion coef", "m2/s", 100 c . iim,jj_nb,nhori, klev,1,klev,nvert, 32, 101 c . "ave(X)", zsto,zout) 108 102 c 109 103 c plusieurs traceurs … … 124 118 CALL histdef(nid_mth, "tops", "Solar rad. at TOA", "W/m2", 125 119 . iim,jj_nb,nhori, 1,1,1, nvert, 32, 126 . "ave(X)", zsto 2,zout)120 . "ave(X)", zsto,zout) 127 121 c 128 122 ENDIF !lev_histmth.GE.2 … … 137 131 CALL histdef(nid_mth, "topl", "IR rad. at TOA", "W/m2", 138 132 . iim,jj_nb,nhori, 1,1,1, nvert, 32, 139 . "ave(X)", zsto 2,zout)133 . "ave(X)", zsto,zout) 140 134 c 141 135 CALL histdef(nid_mth, "sols", "Solar rad. at surf.", "W/m2", 142 136 . iim,jj_nb,nhori, 1,1,1, nvert, 32, 143 . "ave(X)", zsto 2,zout)137 . "ave(X)", zsto,zout) 144 138 c 145 139 CALL histdef(nid_mth, "soll", "IR rad. at surface", "W/m2", 146 140 . iim,jj_nb,nhori, 1,1,1, nvert, 32, 147 . "ave(X)", zsto 2,zout)141 . "ave(X)", zsto,zout) 148 142 c 149 143 c 3D … … 151 145 CALL histdef(nid_mth, "SWnet", "Net SW flux","W/m2", 152 146 . iim,jj_nb,nhori, klev,1,klev,nvert, 153 . 32, "ave(X)", zsto 2,zout)147 . 32, "ave(X)", zsto,zout) 154 148 c 155 149 CALL histdef(nid_mth, "LWnet", "Net LW flux","W/m2", 156 150 . iim,jj_nb,nhori, klev,1,klev,nvert, 157 . 32, "ave(X)", zsto 2,zout)151 . 32, "ave(X)", zsto,zout) 158 152 c 159 153 CALL histdef(nid_mth, "fluxvdf", "PBL net flux","W/m2", -
trunk/LMDZ.VENUS/libf/phyvenus/write_histday.h
r892 r902 4 4 IF (ok_journe) THEN 5 5 6 zsto = dtime7 zout = dtime * REAL(ecrit_day)8 6 itau_w = itau_phy + itap 9 7 … … 54 52 55 53 ENDIF !lev_histday.GE.2 56 c 54 57 55 c------------------------------------------------------- 58 56 IF(lev_histday.GE.3) THEN 59 c 57 60 58 cccccccccccccccccc Radiative transfer 61 59 … … 76 74 77 75 ENDIF !lev_histday.GE.3 78 c 76 79 77 c------------------------------------------------------- 80 78 IF(lev_histday.GE.4) THEN … … 105 103 106 104 ENDIF !lev_histday.GE.4 107 c 105 108 106 c------------------------------------------------------- 109 107 IF(lev_histday.GE.5) THEN … … 116 114 ENDIF !lev_histday.GE.5 117 115 c------------------------------------------------------- 118 c 116 119 117 if (ok_sync) then 120 118 call histsync(nid_day) -
trunk/LMDZ.VENUS/libf/phyvenus/write_histins.h
r892 r902 4 4 IF (ok_instan) THEN 5 5 6 zsto = dtime * ecrit_ins7 zout = dtime * ecrit_ins8 6 itau_w = itau_phy + itap 9 7 … … 107 105 108 106 c------------------------------------------------------- 109 107 IF(lev_histday.GE.5) THEN 110 108 111 109 c call histwrite_phy(nid_ins,.false.,"taux_",itau_w,fluxu) -
trunk/LMDZ.VENUS/libf/phyvenus/write_histmth.h
r892 r902 4 4 IF (ok_mensuel) THEN 5 5 6 zsto = dtime7 zout = dtime * ecrit_mth8 6 itau_w = itau_phy + itap 9 7 … … 105 103 106 104 ENDIF !lev_histmth.GE.4 107 c 105 108 106 c------------------------------------------------------- 109 107 IF(lev_histmth.GE.5) THEN 110 108 111 109 c call histwrite_phy(nid_mth,.false.,"taux_",itau_w,fluxu) … … 116 114 ENDIF !lev_histmth.GE.5 117 115 c------------------------------------------------------- 118 c 116 119 117 if (ok_sync) then 120 118 call histsync(nid_mth)
Note: See TracChangeset
for help on using the changeset viewer.
