Changeset 696 for LMDZ4/trunk/libf/dyn3d
- Timestamp:
- Apr 21, 2006, 10:41:15 AM (19 years ago)
- File:
-
- 1 edited
-
LMDZ4/trunk/libf/dyn3d/write_paramLMDZ_dyn.h (modified) (1 diff)
Legend:
- Unmodified
- Added
- Removed
-
LMDZ4/trunk/libf/dyn3d/write_paramLMDZ_dyn.h
r692 r696 1 1 c 2 dt_cum = dtvr*day_step 3 4 zan = annee_ref 5 dayref = day_ref 6 CALL ymds2ju(zan, 1, dayref, 0.0, zjulian) 7 tau0 = itau_dyn 8 c 9 pi = 4.0 * ATAN(1.0) 10 degres = 180./pi 11 rlong = rlonu * degres 12 rlatg = rlatu * degres 13 c 14 CALL histbeg("paramLMDZ_dyn.nc", 15 . iip1,rlong, jjp1,rlatg, 16 . 1,1,1,1, 17 . tau0, zjulian, dt_cum, 18 . thoriid, nid_ctesGCM) 19 c 20 CALL histdef(nid_ctesGCM, "prt_level", 21 . "Niveau impression debuggage dynamique", 22 . "-",iip1,jjp1,thoriid, 1,1,1, -99, 32, 23 . "once", dt_cum,dt_cum) 24 c 25 CALL histdef(nid_ctesGCM, "dayref", 26 . "Jour de l etat initial ( = 350 si 20 Decembre par ex.)", 27 . "-",iip1,jjp1,thoriid, 1,1,1, -99, 32, 28 . "once", dt_cum,dt_cum) 29 c 30 CALL histdef(nid_ctesGCM, "anneeref", 31 . "Annee de l etat initial", 32 . "-",iip1,jjp1,thoriid, 1,1,1, -99, 32, 33 . "once", dt_cum,dt_cum) 34 c 35 CALL histdef(nid_ctesGCM, "anneelim", 36 . "Annee du fichier limitxxxx.nc si ok_limitvrai =y", 37 . "-",iip1,jjp1,thoriid, 1,1,1, -99, 32, 38 . "once", dt_cum,dt_cum) 39 c 40 CALL histdef(nid_ctesGCM, "raz_date", 41 . "Remise a zero (raz) date init.: 0 pas de raz;1=date gcm.def", 42 . "-",iip1,jjp1,thoriid, 1,1,1, -99, 32, 43 . "once", dt_cum,dt_cum) 44 c 45 CALL histdef(nid_ctesGCM, "nday", 46 . "Nombre de jours d integration", 47 . "-",iip1,jjp1,thoriid, 1,1,1, -99, 32, 48 . "once", dt_cum,dt_cum) 49 c 50 CALL histdef(nid_ctesGCM, "day_step", 51 . "nombre de pas par jour pour dt = 1 min", 52 . "-",iip1,jjp1,thoriid, 1,1,1, -99, 32, 53 . "once", dt_cum,dt_cum) 54 c 55 CALL histdef(nid_ctesGCM, "iperiod", 56 . "periode pour le pas Matsuno (en pas de temps)", 57 . "-",iip1,jjp1,thoriid, 1,1,1, -99, 32, 58 . "once", dt_cum,dt_cum) 59 c 60 CALL histdef(nid_ctesGCM, "iapp_tracvl", 61 . "frequence du groupement des flux (en pas de temps)", 62 . "-",iip1,jjp1,thoriid, 1,1,1, -99, 32, 63 . "once", dt_cum,dt_cum) 64 c 65 CALL histdef(nid_ctesGCM, "iconser", 66 . "periode de sortie des variables de controle (en pas de temps)", 67 . "-",iip1,jjp1,thoriid, 1,1,1, -99, 32, 68 . "once", dt_cum,dt_cum) 69 c 70 CALL histdef(nid_ctesGCM, "iecri", 71 . "periode d ecriture du fichier histoire (en jour)", 72 . "-",iip1,jjp1,thoriid, 1,1,1, -99, 32, 73 . "once", dt_cum,dt_cum) 74 c 75 CALL histdef(nid_ctesGCM, "periodav", 76 . "periode de stockage fichier histmoy (en jour)", 77 . "-",iip1,jjp1,thoriid, 1,1,1, -99, 32, 78 . "once", dt_cum,dt_cum) 79 c 80 CALL histdef(nid_ctesGCM, "idissip", 81 . "periode de la dissipation (en pas) ... a completer", 82 . "-",iip1,jjp1,thoriid, 1,1,1, -99, 32, 83 . "once", dt_cum,dt_cum) 84 c 85 CALL histdef(nid_ctesGCM, "lstardis", 86 . "choix de l operateur de dissipation: 1= star,0=non-star ??", 87 . "-",iip1,jjp1,thoriid, 1,1,1, -99, 32, 88 . "once", dt_cum,dt_cum) 89 c 90 CALL histdef(nid_ctesGCM, "nitergdiv", 91 . "nombre d iterations de l operateur de dissipation gradiv", 92 . "-",iip1,jjp1,thoriid, 1,1,1, -99, 32, 93 . "once", dt_cum,dt_cum) 94 c 95 CALL histdef(nid_ctesGCM, "nitergrot", 96 . "nombre d iterations de l operateur de dissipation nxgradrot", 97 . "-",iip1,jjp1,thoriid, 1,1,1, -99, 32, 98 . "once", dt_cum,dt_cum) 99 c 100 CALL histdef(nid_ctesGCM, "niterh", 101 . "nombre d iterations de l operateur de dissipation divgrad", 102 . "-",iip1,jjp1,thoriid, 1,1,1, -99, 32, 103 . "once", dt_cum,dt_cum) 104 c 105 CALL histdef(nid_ctesGCM, "tetagdiv", 106 ."temps dissipation des + petites long. d ondes pour u,v (gradiv)", 107 . "-",iip1,jjp1,thoriid, 1,1,1, -99, 32, 108 . "once", dt_cum,dt_cum) 109 c 110 CALL histdef(nid_ctesGCM, "tetagrot", 111 ."temps diss. des + petites long. d ondes pour u,v (nxgradrot)", 112 . "-",iip1,jjp1,thoriid, 1,1,1, -99, 32, 113 . "once", dt_cum,dt_cum) 114 c 115 CALL histdef(nid_ctesGCM, "tetatemp", 116 ."temps diss. des + petites long. d ondes pour h (divgrad)", 117 . "-",iip1,jjp1,thoriid, 1,1,1, -99, 32, 118 . "once", dt_cum,dt_cum) 119 c 120 CALL histdef(nid_ctesGCM, "coefdis", 121 ."coefficient pour gamdissip", 122 . "-",iip1,jjp1,thoriid, 1,1,1, -99, 32, 123 . "once", dt_cum,dt_cum) 124 c 125 CALL histdef(nid_ctesGCM, "purmats", 126 ."Choix schema integration temporel: 1=Matsuno,0=Matsuno-leapfrog", 127 . "-",iip1,jjp1,thoriid, 1,1,1, -99, 32, 128 . "once", dt_cum,dt_cum) 129 c 130 CALL histdef(nid_ctesGCM, "ok_guide", 131 ."Guidage: 1=true ,0=false", 132 . "-",iip1,jjp1,thoriid, 1,1,1, -99, 32, 133 . "once", dt_cum,dt_cum) 134 c 135 CALL histdef(nid_ctesGCM, "true_calendar", 136 ."Choix du calendrier: 1=gregorien ,0=calen. a 360 j", 137 . "-",iip1,jjp1,thoriid, 1,1,1, -99, 32, 138 . "once", dt_cum,dt_cum) 139 c 140 CALL histdef(nid_ctesGCM, "guide_calend", 141 ."Guidage calendrier gregorien: 1=oui ,0=non", 142 . "-",iip1,jjp1,thoriid, 1,1,1, -99, 32, 143 . "once", dt_cum,dt_cum) 144 c 145 CALL histdef(nid_ctesGCM, "iflag_phys", 146 ."Permet de faire tourner le modele sans physique: 1=avec ,0=sans", 147 . "-",iip1,jjp1,thoriid, 1,1,1, -99, 32, 148 . "once", dt_cum,dt_cum) 149 c 150 CALL histdef(nid_ctesGCM, "iphysiq", 151 ."Periode de la physique en pas de temps de la dynamique", 152 . "-",iip1,jjp1,thoriid, 1,1,1, -99, 32, 153 . "once", dt_cum,dt_cum) 154 c 155 CALL histdef(nid_ctesGCM, "ecritphy", 156 ."frequence de l ecriture du fichier histphy en jours", 157 . "-",iip1,jjp1,thoriid, 1,1,1, -99, 32, 158 . "once", dt_cum,dt_cum) 159 c 160 CALL histdef(nid_ctesGCM, "cycle_diurne", 161 ."Avoir un cycle diurne: 1=true, 0=false", 162 . "-",iip1,jjp1,thoriid, 1,1,1, -99, 32, 163 . "once", dt_cum,dt_cum) 164 c 165 CALL histdef(nid_ctesGCM, "soil_model", 166 ."Choix du modele de sol (Thermique ?): 1=true, 0=false", 167 . "-",iip1,jjp1,thoriid, 1,1,1, -99, 32, 168 . "once", dt_cum,dt_cum) 169 c 170 CALL histdef(nid_ctesGCM, "new_oliq", 171 ."Nouvelle eau liquide: 1=true, 0=false", 172 . "-",iip1,jjp1,thoriid, 1,1,1, -99, 32, 173 . "once", dt_cum,dt_cum) 174 c 175 CALL histdef(nid_ctesGCM, "ok_orodr", 176 ."Orography drag : 1=true, 0=false", 177 . "-",iip1,jjp1,thoriid, 1,1,1, -99, 32, 178 . "once", dt_cum,dt_cum) 179 c 180 CALL histdef(nid_ctesGCM, "ok_orolf", 181 ."Orography lift : 1=true, 0=false", 182 . "-",iip1,jjp1,thoriid, 1,1,1, -99, 32, 183 . "once", dt_cum,dt_cum) 184 c 185 CALL histdef(nid_ctesGCM, "ok_limitvrai", 186 ."Force la lecture de la bonne annee: 1=true, 0=false", 187 . "-",iip1,jjp1,thoriid, 1,1,1, -99, 32, 188 . "once", dt_cum,dt_cum) 189 c 190 CALL histdef(nid_ctesGCM, "nbapp_rad", 191 ."Nombre d appels des routines de rayonnements par jour", 192 . "-",iip1,jjp1,thoriid, 1,1,1, -99, 32, 193 . "once", dt_cum,dt_cum) 194 c 195 CALL histdef(nid_ctesGCM, "iflag_con", 196 ."Flag pour la convection: 1=LMD, 2=Tiedtke, 3=KE, 4=KE vect.", 197 . "-",iip1,jjp1,thoriid, 1,1,1, -99, 32, 198 . "once", dt_cum,dt_cum) 199 c 200 CALL histdef(nid_ctesGCM, "clon", 201 ."longitude en degres du centre du zoom", 202 . "deg",iip1,jjp1,thoriid, 1,1,1, -99, 32, 203 . "once", dt_cum,dt_cum) 204 c 205 CALL histdef(nid_ctesGCM, "clat", 206 ."latitude en degres du centre du zoom", 207 . "deg",iip1,jjp1,thoriid, 1,1,1, -99, 32, 208 . "once", dt_cum,dt_cum) 209 c 210 CALL histdef(nid_ctesGCM, "grossismx", 211 ."facteur de grossissement du zoom, selon la longitude", 212 . "-",iip1,jjp1,thoriid, 1,1,1, -99, 32, 213 . "once", dt_cum,dt_cum) 214 c 215 CALL histdef(nid_ctesGCM, "grossismy", 216 ."facteur de grossissement du zoom, selon la latitude", 217 . "-",iip1,jjp1,thoriid, 1,1,1, -99, 32, 218 . "once", dt_cum,dt_cum) 219 c 220 CALL histdef(nid_ctesGCM, "fxyhypb", 221 ."Fonction f(y) hyperbolique si true=1, sinusoidale si false=0", 222 . "-",iip1,jjp1,thoriid, 1,1,1, -99, 32, 223 . "once", dt_cum,dt_cum) 224 c 225 CALL histdef(nid_ctesGCM, "dzoomx", 226 ."extension en longitude de la zone du zoom (fraction zone totale)" 227 . ,"-",iip1,jjp1,thoriid, 1,1,1, -99, 32, 228 . "once", dt_cum,dt_cum) 229 c 230 CALL histdef(nid_ctesGCM, "dzoomy", 231 ."extension en latitude de la zone du zoom (fraction zone totale)" 232 . ,"-",iip1,jjp1,thoriid, 1,1,1, -99, 32, 233 . "once", dt_cum,dt_cum) 234 c 235 CALL histdef(nid_ctesGCM, "taux", 236 ."raideur du zoom en X" 237 . ,"-",iip1,jjp1,thoriid, 1,1,1, -99, 32, 238 . "once", dt_cum,dt_cum) 239 c 240 CALL histdef(nid_ctesGCM, "tauy", 241 ."raideur du zoom en Y" 242 . ,"-",iip1,jjp1,thoriid, 1,1,1, -99, 32, 243 . "once", dt_cum,dt_cum) 244 c 245 CALL histdef(nid_ctesGCM, "ysinus", 246 ."ysinus=1: Ftion f(y) avec y=Sin(latit.)/ ysinus=0: y = latit" 247 . ,"-",iip1,jjp1,thoriid, 1,1,1, -99, 32, 248 . "once", dt_cum,dt_cum) 249 c 250 CALL histdef(nid_ctesGCM, "ip_ebil_dyn", 251 ."PRINTlevel for energy conservation diag.; 0/1= pas de print, 252 . 2= print","-",iip1,jjp1,thoriid, 1,1,1, -99, 32, 253 . "once", dt_cum,dt_cum) 2 ndex2d = 0 3 itau_dyn=itau_dyn+itau 4 print*,'write_paramLMDZ_dyn :itau_dyn=',itau_dyn 5 c 6 zx_tmp_2d(1:iip1,1:jjp1)=FLOAT(prt_level) 7 CALL histwrite(nid_ctesGCM, "prt_level", itau_dyn, 8 . zx_tmp_2d,iip1*jjp1,ndex2d) 9 c 10 zx_tmp_2d(1:iip1,1:jjp1)=FLOAT(dayref) 11 CALL histwrite(nid_ctesGCM, "dayref", itau_dyn, 12 . zx_tmp_2d,iip1*jjp1,ndex2d) 13 c 14 zx_tmp_2d(1:iip1,1:jjp1)=FLOAT(anneeref) 15 CALL histwrite(nid_ctesGCM, "anneeref", itau_dyn, 16 . zx_tmp_2d,iip1*jjp1,ndex2d) 17 c 18 zx_tmp_2d(1:iip1,1:jjp1)=FLOAT(raz_date) 19 CALL histwrite(nid_ctesGCM, "raz_date", itau_dyn, 20 . zx_tmp_2d,iip1*jjp1,ndex2d) 21 c 22 zx_tmp_2d(1:iip1,1:jjp1)=FLOAT(nday) 23 CALL histwrite(nid_ctesGCM, "nday", itau_dyn, 24 . zx_tmp_2d,iip1*jjp1,ndex2d) 25 c 26 zx_tmp_2d(1:iip1,1:jjp1)=FLOAT(day_step) 27 CALL histwrite(nid_ctesGCM, "day_step", itau_dyn, 28 . zx_tmp_2d,iip1*jjp1,ndex2d) 29 c 30 zx_tmp_2d(1:iip1,1:jjp1)=FLOAT(iperiod) 31 CALL histwrite(nid_ctesGCM, "iperiod", itau_dyn, 32 . zx_tmp_2d,iip1*jjp1,ndex2d) 33 c 34 zx_tmp_2d(1:iip1,1:jjp1)=FLOAT(iapp_tracvl) 35 CALL histwrite(nid_ctesGCM, "iapp_tracvl", itau_dyn, 36 . zx_tmp_2d,iip1*jjp1,ndex2d) 37 c 38 zx_tmp_2d(1:iip1,1:jjp1)=FLOAT(iconser) 39 CALL histwrite(nid_ctesGCM, "iconser", itau_dyn, 40 . zx_tmp_2d,iip1*jjp1,ndex2d) 41 c 42 zx_tmp_2d(1:iip1,1:jjp1)=FLOAT(iecri) 43 CALL histwrite(nid_ctesGCM, "iecri", itau_dyn, 44 . zx_tmp_2d,iip1*jjp1,ndex2d) 45 c 46 zx_tmp_2d(1:iip1,1:jjp1)=periodav 47 CALL histwrite(nid_ctesGCM, "periodav", itau_dyn, 48 . zx_tmp_2d,iip1*jjp1,ndex2d) 49 c 50 zx_tmp_2d(1:iip1,1:jjp1)=FLOAT(idissip) 51 CALL histwrite(nid_ctesGCM, "idissip", itau_dyn, 52 . zx_tmp_2d,iip1*jjp1,ndex2d) 53 c 54 IF(lstardis) THEN 55 zx_tmp_2d(1:iip1,1:jjp1)=1. 56 ELSE 57 zx_tmp_2d(1:iip1,1:jjp1)=0. 58 ENDIF 59 CALL histwrite(nid_ctesGCM, "lstardis", itau_dyn, 60 . zx_tmp_2d,iip1*jjp1,ndex2d) 61 c 62 zx_tmp_2d(1:iip1,1:jjp1)=FLOAT(nitergdiv) 63 CALL histwrite(nid_ctesGCM, "nitergdiv", itau_dyn, 64 . zx_tmp_2d,iip1*jjp1,ndex2d) 65 c 66 zx_tmp_2d(1:iip1,1:jjp1)=FLOAT(nitergrot) 67 CALL histwrite(nid_ctesGCM, "nitergrot", itau_dyn, 68 . zx_tmp_2d,iip1*jjp1,ndex2d) 69 c 70 zx_tmp_2d(1:iip1,1:jjp1)=FLOAT(niterh) 71 CALL histwrite(nid_ctesGCM, "niterh", itau_dyn, 72 . zx_tmp_2d,iip1*jjp1,ndex2d) 73 c 74 zx_tmp_2d(1:iip1,1:jjp1)=tetagdiv 75 CALL histwrite(nid_ctesGCM, "tetagdiv", itau_dyn, 76 . zx_tmp_2d,iip1*jjp1,ndex2d) 77 c 78 zx_tmp_2d(1:iip1,1:jjp1)=tetagrot 79 CALL histwrite(nid_ctesGCM, "tetagrot", itau_dyn, 80 . zx_tmp_2d,iip1*jjp1,ndex2d) 81 c 82 zx_tmp_2d(1:iip1,1:jjp1)=tetatemp 83 CALL histwrite(nid_ctesGCM, "tetatemp", itau_dyn, 84 . zx_tmp_2d,iip1*jjp1,ndex2d) 85 c 86 zx_tmp_2d(1:iip1,1:jjp1)=coefdis 87 CALL histwrite(nid_ctesGCM, "coefdis", itau_dyn, 88 . zx_tmp_2d,iip1*jjp1,ndex2d) 89 c 90 IF(purmats) THEN 91 zx_tmp_2d(1:iip1,1:jjp1)=1. 92 ELSE 93 zx_tmp_2d(1:iip1,1:jjp1)=0. 94 ENDIF 95 CALL histwrite(nid_ctesGCM, "purmats", itau_dyn, 96 . zx_tmp_2d,iip1*jjp1,ndex2d) 97 c 98 IF(ok_guide) THEN 99 zx_tmp_2d(1:iip1,1:jjp1)=1. 100 ELSE 101 zx_tmp_2d(1:iip1,1:jjp1)=0. 102 ENDIF 103 CALL histwrite(nid_ctesGCM, "ok_guide", itau_dyn, 104 . zx_tmp_2d,iip1*jjp1,ndex2d) 105 c 106 IF(true_calendar) THEN 107 zx_tmp_2d(1:iip1,1:jjp1)=1. 108 ELSE 109 zx_tmp_2d(1:iip1,1:jjp1)=0. 110 ENDIF 111 CALL histwrite(nid_ctesGCM, "true_calendar", itau_dyn, 112 . zx_tmp_2d,iip1*jjp1,ndex2d) 113 c 114 zx_tmp_2d(1:iip1,1:jjp1)=FLOAT(iflag_phys) 115 CALL histwrite(nid_ctesGCM, "iflag_phys", itau_dyn, 116 . zx_tmp_2d,iip1*jjp1,ndex2d) 117 c 118 zx_tmp_2d(1:iip1,1:jjp1)=FLOAT(iphysiq) 119 CALL histwrite(nid_ctesGCM, "iphysiq", itau_dyn, 120 . zx_tmp_2d,iip1*jjp1,ndex2d) 121 c 122 zx_tmp_2d(1:iip1,1:jjp1)=FLOAT(ecritphy) 123 CALL histwrite(nid_ctesGCM, "ecritphy", itau_dyn, 124 . zx_tmp_2d,iip1*jjp1,ndex2d) 125 c 126 IF(cycle_diurne) THEN 127 zx_tmp_2d(1:iip1,1:jjp1)=1. 128 ELSE 129 zx_tmp_2d(1:iip1,1:jjp1)=0. 130 ENDIF 131 CALL histwrite(nid_ctesGCM, "cycle_diurne", itau_dyn, 132 . zx_tmp_2d,iip1*jjp1,ndex2d) 133 c 134 IF(soil_model) THEN 135 zx_tmp_2d(1:iip1,1:jjp1)=1. 136 ELSE 137 zx_tmp_2d(1:iip1,1:jjp1)=0. 138 ENDIF 139 CALL histwrite(nid_ctesGCM, "soil_model", itau_dyn, 140 . zx_tmp_2d,iip1*jjp1,ndex2d) 141 c 142 IF(new_oliq) THEN 143 zx_tmp_2d(1:iip1,1:jjp1)=1. 144 ELSE 145 zx_tmp_2d(1:iip1,1:jjp1)=0. 146 ENDIF 147 CALL histwrite(nid_ctesGCM, "new_oliq", itau_dyn, 148 . zx_tmp_2d,iip1*jjp1,ndex2d) 149 c 150 IF(ok_orodr) THEN 151 zx_tmp_2d(1:iip1,1:jjp1)=1. 152 ELSE 153 zx_tmp_2d(1:iip1,1:jjp1)=0. 154 ENDIF 155 CALL histwrite(nid_ctesGCM, "ok_orodr", itau_dyn, 156 . zx_tmp_2d,iip1*jjp1,ndex2d) 157 c 158 IF(ok_orolf) THEN 159 zx_tmp_2d(1:iip1,1:jjp1)=1. 160 ELSE 161 zx_tmp_2d(1:iip1,1:jjp1)=0. 162 ENDIF 163 CALL histwrite(nid_ctesGCM, "ok_orolf", itau_dyn, 164 . zx_tmp_2d,iip1*jjp1,ndex2d) 165 c 166 IF(ok_limitvrai) THEN 167 zx_tmp_2d(1:iip1,1:jjp1)=1. 168 ELSE 169 zx_tmp_2d(1:iip1,1:jjp1)=0. 170 ENDIF 171 CALL histwrite(nid_ctesGCM, "ok_limitvrai", itau_dyn, 172 . zx_tmp_2d,iip1*jjp1,ndex2d) 173 c 174 zx_tmp_2d(1:iip1,1:jjp1)=FLOAT(nbapp_rad) 175 CALL histwrite(nid_ctesGCM, "nbapp_rad", itau_dyn, 176 . zx_tmp_2d,iip1*jjp1,ndex2d) 177 c 178 zx_tmp_2d(1:iip1,1:jjp1)=FLOAT(iflag_con) 179 CALL histwrite(nid_ctesGCM, "iflag_con", itau_dyn, 180 . zx_tmp_2d,iip1*jjp1,ndex2d) 181 c 182 zx_tmp_2d(1:iip1,1:jjp1)=clon 183 CALL histwrite(nid_ctesGCM, "clon", itau_dyn, 184 . zx_tmp_2d,iip1*jjp1,ndex2d) 185 c 186 zx_tmp_2d(1:iip1,1:jjp1)=clat 187 CALL histwrite(nid_ctesGCM, "clat", itau_dyn, 188 . zx_tmp_2d,iip1*jjp1,ndex2d) 189 c 190 zx_tmp_2d(1:iip1,1:jjp1)=grossismx 191 CALL histwrite(nid_ctesGCM, "grossismx", itau_dyn, 192 . zx_tmp_2d,iip1*jjp1,ndex2d) 193 c 194 zx_tmp_2d(1:iip1,1:jjp1)=grossismy 195 CALL histwrite(nid_ctesGCM, "grossismy", itau_dyn, 196 . zx_tmp_2d,iip1*jjp1,ndex2d) 197 c 198 IF(fxyhypb) THEN 199 zx_tmp_2d(1:iip1,1:jjp1)=1. 200 ELSE 201 zx_tmp_2d(1:iip1,1:jjp1)=0. 202 ENDIF 203 CALL histwrite(nid_ctesGCM, "fxyhypb", itau_dyn, 204 . zx_tmp_2d,iip1*jjp1,ndex2d) 205 c 206 zx_tmp_2d(1:iip1,1:jjp1)=dzoomx 207 CALL histwrite(nid_ctesGCM, "dzoomx", itau_dyn, 208 . zx_tmp_2d,iip1*jjp1,ndex2d) 209 c 210 zx_tmp_2d(1:iip1,1:jjp1)=dzoomy 211 CALL histwrite(nid_ctesGCM, "dzoomy", itau_dyn, 212 . zx_tmp_2d,iip1*jjp1,ndex2d) 213 c 214 zx_tmp_2d(1:iip1,1:jjp1)=taux 215 CALL histwrite(nid_ctesGCM, "taux", itau_dyn, 216 . zx_tmp_2d,iip1*jjp1,ndex2d) 217 c 218 zx_tmp_2d(1:iip1,1:jjp1)=tauy 219 CALL histwrite(nid_ctesGCM, "tauy", itau_dyn, 220 . zx_tmp_2d,iip1*jjp1,ndex2d) 221 c 222 IF(ysinus) THEN 223 zx_tmp_2d(1:iip1,1:jjp1)=1. 224 ELSE 225 zx_tmp_2d(1:iip1,1:jjp1)=0. 226 ENDIF 227 CALL histwrite(nid_ctesGCM, "ysinus", itau_dyn, 228 . zx_tmp_2d,iip1*jjp1,ndex2d) 229 c 230 zx_tmp_2d(1:iip1,1:jjp1)=FLOAT(ip_ebil_dyn) 231 CALL histwrite(nid_ctesGCM, "ip_ebil_dyn", itau_dyn, 232 . zx_tmp_2d,iip1*jjp1,ndex2d) 254 233 c 255 234 c================================================================= 256 235 c 257 CALL histend(nid_ctesGCM) 236 if (ok_sync) then 237 call histsync(nid_ctesGCM) 238 endif 258 239 c 259 240 c=================================================================
Note: See TracChangeset
for help on using the changeset viewer.
