Changeset 4514 for branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM
- Timestamp:
- 2014-02-25T12:41:51+01:00 (10 years ago)
- Location:
- branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts
- Files:
-
- 3 added
- 5 deleted
- 12 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_IceSal.pro
r4507 r4514 1 pro std_plot_IceSal, Isal1, Isal2, ARC = arc, ANT = ant, APRIL = april, SEPT = sept, POSTSCRIPT = postscript, _extra = ex1 pro std_plot_IceSal, Isal1, Isal2, ARC = arc, ANT = ant, FEBR = febr, MARCH = march, SEPT = sept, POSTSCRIPT = postscript, _extra = ex 2 2 3 3 compile_opt idl2, strictarrsubs … … 12 12 IF keyword_set(arc) THEN var = var+'_Arc_' 13 13 IF keyword_set(ant) THEN var = var+'_Ant_' 14 IF keyword_set(april) THEN var = var+'April' 14 IF keyword_set(febr) THEN var = var+'Febr' 15 IF keyword_set(march) THEN var = var+'March' 15 16 IF keyword_set(sept) THEN var = var+'Sept' 16 17 … … 31 32 ; 32 33 title = var+'!C'+std_file1_I 33 plt, Isal1.arr, MIN = .1, MAX = 20.1, INTER = 1., /STRICTFILL, CELL_FILL = 2, format = '(f4.1)' $34 plt, Isal1.arr, MIN = 0., MAX = 12., INTER = 1., /STRICTFILL, CELL_FILL = 2, format = '(f4.1)' $ 34 35 , small = [1, 2, 1], COAST_THICK = 2, TITLE = title $ 35 36 , CHARSIZE = 1.05, GLINETHICK = 2., /ORTHO, MAP = map, /PORTRAIT, _extra = ex … … 37 38 if std_file1_I NE std_file2_I then begin 38 39 title = title + std_file2_I 39 plt, Isal1.arr - Isal2.arr, MIN = -4, MAX = 4, INTER = .2, STYLE = 'so0so', format = '(f4.1)' $40 plt, Isal1.arr - Isal2.arr, MIN = -4, MAX = 4, INTER = 1., STYLE = 'so0so', format = '(f4.1)' $ 40 41 , small = [1, 2, 2], COAST_THICK = 2, CELL_FILL = 2, TITLE = title $ 41 42 , CHARSIZE = 1.05, GLINETHICK = 2., /ORTHO, MAP = map, /NOERASE, _extra = ex -
branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_SnowThick.pro
r4501 r4514 1 pro std_plot_SnowThick, Isnow1, Isnow2, ARC = arc, ANT = ant, APRIL = april, SEPT = sept, POSTSCRIPT = postscript, _extra = ex1 pro std_plot_SnowThick, Isnow1, Isnow2, ARC = arc, ANT = ant, MARCH = march, SEPT = sept, POSTSCRIPT = postscript, _extra = ex 2 2 3 3 compile_opt idl2, strictarrsubs … … 12 12 IF keyword_set(arc) THEN var = var+'_Arc_' 13 13 IF keyword_set(ant) THEN var = var+'_Ant_' 14 IF keyword_set( april) THEN var = var+'April'14 IF keyword_set(march) THEN var = var+'March' 15 15 IF keyword_set(sept) THEN var = var+'Sept' 16 16 … … 22 22 domdef, 20, 380, 50, 90 23 23 vmin = 0. 24 vmax = 1.24 vmax = .6 25 25 vint = 0.1 26 26 ;fmt = '(i2)' -
branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_all.pro
r4507 r4514 187 187 ENDELSE 188 188 ; 189 Iage_1 = read_ncdf(getenv('VAR1_Iage'), allrecords = allrec, filename = std_file1_I ) 190 caldat, time, mm 191 march = where(mm EQ 3, cnt) 192 Iage_march_1 = {arr:1./float(cnt) * total(reform(Iage_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Iage_1.unit} 193 sept = where(mm EQ 9, cnt) 194 Iage_sept_1 = {arr:1./float(cnt) * total(reform(Iage_1.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Iage_1.unit} 195 undefine, Iage_1 196 ; 197 IF std_file2_I NE std_file1_I THEN BEGIN 198 Iage_2 = read_ncdf(getenv('VAR2_Iage'), allrecords = allrec, filename = std_file2_I ) 199 caldat, time, mm 200 march = where(mm EQ 3, cnt) 201 Iage_march_2 = {arr:1./float(cnt) * total(reform(Iage_2.arr[*, *, temporary(april)],nxt,nyt,cnt), 3), unit:Iage_2.unit} 202 sept = where(mm EQ 9, cnt) 203 Iage_sept_2 = {arr:1./float(cnt) * total(reform(Iage_2.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Iage_2.unit} 204 undefine, Iage_2 205 ENDIF ELSE BEGIN 206 Iage_march_2 = {arr:-1} 207 Iage_sept_2 = {arr:-1} 208 ENDELSE 209 ; 189 210 Ifra_1 = read_ncdf(getenv('VAR1_Ifrac'), allrecords = allrec, filename = std_file1_I ) 190 ;help, Ifra_1.arr, jpi,jpj191 211 caldat, time, mm 192 212 march = where(mm EQ 3, cnt) … … 210 230 ; 211 231 Isnow_1 = read_ncdf(getenv('VAR1_Isnow'), allrecords = allrec, filename = std_file1_I ) 212 ;help, Isnow_1.arr, jpi,jpj213 232 caldat, time, mm 214 april = where(mm EQ 4, cnt)215 Isnow_ april_1 = {arr:1./float(cnt) * total(reform(Isnow_1.arr[*, *, temporary(april)],nxt,nyt,cnt), 3), unit:Isnow_1.unit}233 march = where(mm EQ 3, cnt) 234 Isnow_march_1 = {arr:1./float(cnt) * total(reform(Isnow_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Isnow_1.unit} 216 235 sept = where(mm EQ 9, cnt) 217 236 Isnow_sept_1 = {arr:1./float(cnt) * total(reform(Isnow_1.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Isnow_1.unit} … … 221 240 Isnow_2 = read_ncdf(getenv('VAR2_Isnow'), allrecords = allrec, filename = std_file2_I ) 222 241 caldat, time, mm 223 april = where(mm EQ 4, cnt)224 Isnow_ april_2 = {arr:1./float(cnt) * total(reform(Isnow_2.arr[*, *, temporary(april)],nxt,nyt,cnt), 3), unit:Isnow_2.unit}242 march = where(mm EQ 3, cnt) 243 Isnow_march_2 = {arr:1./float(cnt) * total(reform(Isnow_2.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Isnow_2.unit} 225 244 sept = where(mm EQ 9, cnt) 226 245 Isnow_sept_2 = {arr:1./float(cnt) * total(reform(Isnow_2.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Isnow_2.unit} 227 246 ; undefine, Isnow_2 228 247 ENDIF ELSE BEGIN 229 Isnow_ april_2 = {arr:-1}248 Isnow_march_2 = {arr:-1} 230 249 Isnow_sept_2 = {arr:-1} 231 250 ENDELSE 232 251 ; 233 252 Isal_1 = read_ncdf(getenv('VAR1_Isal'), allrecords = allrec, filename = std_file1_I ) 234 ;help, Isal_1.arr, jpi,jpj 253 Ifra_1 = read_ncdf(getenv('VAR1_Ifrac'), allrecords = allrec, filename = std_file1_I ) 254 msk = Ifra_1.arr gt 0.15 ; remove 0.15% for observations 255 help, Ifra_1 256 help, mask 235 257 caldat, time, mm 236 april = where(mm EQ 4, cnt) 237 Isal_april_1 = {arr:1./float(cnt) * total(reform(Isal_1.arr[*, *, temporary(april)],nxt,nyt,cnt), 3), unit:Isal_1.unit} 258 march = where(mm EQ 3, cnt) 259 febr = where(mm EQ 2, cnt) 260 ;SF Isal_march_1 = {arr:1./float(cnt) * total(reform(Isal_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Isal_1.unit} 261 Isal_1.arr = Isal_1.arr * msk 262 help, Isal_1 263 Isal_march_1 = {arr:1./float(cnt) * total(reform(Isal_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Isal_1.unit} 264 help, Isal_march_1 265 Isal_febr_1 = {arr:1./float(cnt) * total(reform(Isal_1.arr[*, *, temporary(febr)],nxt,nyt,cnt), 3), unit:Isal_1.unit} 238 266 sept = where(mm EQ 9, cnt) 239 267 Isal_sept_1 = {arr:1./float(cnt) * total(reform(Isal_1.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Isal_1.unit} 240 268 ;undefine, Isal_1 241 269 ; 270 Isal_2 = read_ncdf(getenv('VAR2_Isal'), allrecords = allrec, filename = std_file2_I ) 271 Ifra_2 = read_ncdf(getenv('VAR2_Ifrac'), allrecords = allrec, filename = std_file2_I ) 272 msk = Ifra_2.arr gt 0.15 ; remove 0.15% for observations 242 273 IF std_file2_I NE std_file1_I THEN BEGIN 243 274 Isal_2 = read_ncdf(getenv('VAR2_Isal'), allrecords = allrec, filename = std_file2_I ) 275 Isal_2.arr = Isal_2.arr * msk 244 276 caldat, time, mm 245 april = where(mm EQ 4, cnt) 246 Isal_april_2 = {arr:1./float(cnt) * total(reform(Isal_2.arr[*, *, temporary(april)],nxt,nyt,cnt), 3), unit:Isal_2.unit} 277 march = where(mm EQ 3, cnt) 278 febr = where(mm EQ 2, cnt) 279 Isal_march_2 = {arr:1./float(cnt) * total(reform(Isal_2.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Isal_2.unit} 280 Isal_febr_2 = {arr:1./float(cnt) * total(reform(Isal_2.arr[*, *, temporary(febr)],nxt,nyt,cnt), 3), unit:Isal_2.unit} 247 281 sept = where(mm EQ 9, cnt) 248 282 Isal_sept_2 = {arr:1./float(cnt) * total(reform(Isal_2.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Isal_2.unit} 249 283 ; undefine, Isal_2 250 284 ENDIF ELSE BEGIN 251 Isal_april_2 = {arr:-1} 285 Isal_febr_2 = {arr:-1} 286 Isal_march_2 = {arr:-1} 252 287 Isal_sept_2 = {arr:-1} 253 288 ENDELSE … … 375 410 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_sept_1, Ithi_sept_2, /ANT, /SEPT, _extra = ex 376 411 ; 412 cnt = cnt+1 & blabla = 'Arctic Ice Age: MARCH' 413 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceAge, Iage_march_1, Iage_march_2, /ARC, /MARCH, _extra = ex 414 ; 415 cnt = cnt+1 & blabla = 'Arctic Ice Age: SEPT' 416 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceAge, Iage_sept_1, Iage_sept_2, /ARC, /SEPT, _extra = ex 417 ; 418 cnt = cnt+1 & blabla = 'Antarctic Ice Age: MARCH' 419 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceAge, Iage_march_1, Iage_march_2, /ANT, /MARCH, _extra = ex 420 ; 421 cnt = cnt+1 & blabla = 'Antarctic Ice Age: SEPT' 422 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceAge, Iage_sept_1, Iage_sept_2, /ANT, /SEPT, _extra = ex 423 ; 377 424 cnt = cnt+1 & blabla = 'Arctic Ice Fraction: MARCH' 378 425 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceFrac, Ifra_march_1, Ifra_march_2, /ARC, /MARCH, _extra = ex … … 387 434 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceFrac, Ifra_sept_1, Ifra_sept_2, /ANT, /SEPT, _extra = ex 388 435 ; 389 cnt = cnt+1 & blabla = 'Arctic SNOW Thickness: APRIL'390 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_ april_1, Isnow_april_2, /ARC, /APRIL, _extra = ex436 cnt = cnt+1 & blabla = 'Arctic SNOW Thickness: MARCH' 437 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_march_1, Isnow_march_2, /ARC, /MARCH, _extra = ex 391 438 ; 392 439 cnt = cnt+1 & blabla = 'Arctic SNOW Thickness: SEPT' 393 440 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_sept_1, Isnow_sept_2, /ARC, /SEPT, _extra = ex 394 441 ; 395 cnt = cnt+1 & blabla = 'Antarctic SNOW Thickness: APRIL'396 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_ april_1, Isnow_april_2, /ANT, /APRIL, _extra = ex442 cnt = cnt+1 & blabla = 'Antarctic SNOW Thickness: MARCH' 443 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_march_1, Isnow_march_2, /ANT, /MARCH, _extra = ex 397 444 ; 398 445 cnt = cnt+1 & blabla = 'Antarctic SNOW Thickness: SEPT' 399 446 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_sept_1, Isnow_sept_2, /ANT, /SEPT, _extra = ex 400 447 ; 401 cnt = cnt+1 & blabla = 'Arctic Ice Salinity: APRIL'402 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceSal, Isal_ april_1, Isal_april_2, /ARC, /APRIL, _extra = ex448 cnt = cnt+1 & blabla = 'Arctic Ice Salinity: MARCH' 449 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceSal, Isal_march_1, Isal_march_2, /ARC, /MARCH, _extra = ex 403 450 ; 404 451 cnt = cnt+1 & blabla = 'Arctic Ice Salinity: SEPT' 405 452 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceSal, Isal_sept_1, Isal_sept_2, /ARC, /SEPT, _extra = ex 406 453 ; 407 cnt = cnt+1 & blabla = 'Antarctic Ice Salinity: APRIL'408 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceSal, Isal_ april_1, Isal_april_2, /ANT, /APRIL, _extra = ex454 cnt = cnt+1 & blabla = 'Antarctic Ice Salinity: FEBRUARY' 455 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceSal, Isal_febr_1, Isal_febr_2, /ANT, /FEBR, _extra = ex 409 456 ; 410 457 cnt = cnt+1 & blabla = 'Antarctic Ice Salinity: SEPT' -
branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_all_1Y.pro
r4320 r4514 115 115 EMP1 = read_ncdf(getenv('VAR1_EMP'), allrecords = allrec, direc = 't', filename = std_file1_T ) 116 116 EMP1 = {arr:EMP1.arr * 86400., unit:'mm/day', grid:'T'} 117 IF std_file2_T NE std_file1_T THEN BEGIN 118 EMP2 = read_ncdf(getenv('VAR2_EMP'), allrecords = allrec, direc = 't', filename = std_file2_T ) 119 EMP2 = {arr:EMP2.arr * 86400., unit:'mm/day', grid:'T'} 120 ENDIF ELSE EMP2 = {arr:-1} 121 122 ;mixed layer depth 123 MLD1 = read_ncdf(getenv('VAR1_MLD'), allrecords = allrec, direc = 't', filename = std_file1_T ) ; 10 m 124 IF std_file2_T NE std_file1_T THEN BEGIN 125 MLD2 = read_ncdf(getenv('VAR2_MLD'), allrecords = allrec, direc = 't', filename = std_file2_T ) ; 10 m 126 ENDIF ELSE MLD2 = {arr:-1} 127 ;climatology 128 MLD = read_ncdf(getenv('VAR_MLD'), filename = std_file_mld ) 129 130 ; velocities 131 U1 = read_ncdf(getenv('VAR1_U'), allrecords = allrec, direc = 't', filename = std_file1_U ) 132 IF strlowcase(getenv('VAR1_U')) EQ 'uocetr_eff' THEN BEGIN 133 U1.arr = U1.arr / e3u_3d(/e2) * umask() 134 U1.unit = 'm/s' 135 ENDIF 136 IF std_file2_U NE std_file1_U THEN BEGIN 137 U2 = read_ncdf(getenv('VAR2_U'), allrecords = allrec, direc = 't', filename = std_file2_U ) 138 IF strlowcase(getenv('VAR2_U')) EQ 'uocetr_eff' THEN BEGIN 139 U2.arr = U2.arr / e3u_3d(/e2) * umask() 140 U2.unit = 'm/s' 141 ENDIF 142 ENDIF ELSE U2 = {arr:-1} 143 ; 144 V1 = read_ncdf(getenv('VAR1_V'), allrecords = allrec, direc = 't', filename = std_file1_V ) 145 IF strlowcase(getenv('VAR1_V')) EQ 'vocetr_eff' THEN BEGIN 146 V1.arr = V1.arr / e3v_3d(/e1) * vmask() 147 V1.unit = 'm/s' 148 ENDIF 149 IF std_file2_V NE std_file1_V THEN BEGIN 150 V2 = read_ncdf(getenv('VAR2_V'), allrecords = allrec, direc = 't', filename = std_file2_V ) 151 IF strlowcase(getenv('VAR2_V')) EQ 'vocetr_eff' THEN BEGIN 152 V2.arr = V2.arr / e3v_3d(/e1) * vmask() 153 V2.unit = 'm/s' 154 ENDIF 155 ENDIF ELSE V2 = {arr:-1} 117 IF std_file2_T NE std_file1_T THEN BEGIN 118 EMP2 = read_ncdf(getenv('VAR2_EMP'), allrecords = allrec, direc = 't', filename = std_file2_T ) 119 EMP2 = {arr:EMP2.arr * 86400., unit:'mm/day', grid:'T'} 120 ENDIF ELSE EMP2 = {arr:-1} 121 122 ;mixed layer depth 123 MLD1 = read_ncdf(getenv('VAR1_MLD'), allrecords = allrec, direc = 't', filename = std_file1_T ) ; 10 m 124 IF std_file2_T NE std_file1_T THEN BEGIN 125 MLD2 = read_ncdf(getenv('VAR2_MLD'), allrecords = allrec, direc = 't', filename = std_file2_T ) ; 10 m 126 ENDIF ELSE MLD2 = {arr:-1} 127 ;climatology 128 MLD = read_ncdf(getenv('VAR_MLD'), filename = std_file_mld ) 129 130 ; velocities 131 U1 = read_ncdf(getenv('VAR1_U'), allrecords = allrec, direc = 't', filename = std_file1_U ) 132 ; old formulation: we tested variable name 133 ; IF strlowcase(getenv('VAR1_U')) EQ 'uocetr_eff' OR strlowcase(getenv('VAR1_U')) EQ 'vozoeftr' THEN BEGIN 134 IF strlowcase(U1.unit) EQ 'm3/s' THEN BEGIN 135 ;IF it is a transport it is transofrmed in velocity 136 U1.arr = U1.arr / e3u_3d(/e2) * umask() 137 U1.unit = 'm/s' 138 ENDIF 139 IF std_file2_U NE std_file1_U THEN BEGIN 140 U2 = read_ncdf(getenv('VAR2_U'), allrecords = allrec, direc = 't', filename = std_file2_U ) 141 ; old formulation: we tested variable name 142 ; IF strlowcase(getenv('VAR2_U')) EQ 'uocetr_eff' OR strlowcase(getenv('VAR2_U')) EQ 'vozoeftr' THEN BEGIN 143 IF strlowcase(U2.unit) EQ 'm3/s' THEN BEGIN 144 U2.arr = U2.arr / e3u_3d(/e2) * umask() 145 U2.unit = 'm/s' 146 ENDIF 147 ENDIF ELSE U2 = {arr:-1} 148 ; 149 V1 = read_ncdf(getenv('VAR1_V'), allrecords = allrec, direc = 't', filename = std_file1_V ) 150 ; old formulation: we tested variable name 151 ; IF strlowcase(getenv('VAR1_V')) EQ 'vocetr_eff' OR strlowcase(getenv('VAR1_V')) EQ 'vomeeftr' THEN BEGIN 152 IF strlowcase(V1.unit) EQ 'm3/s' THEN BEGIN ; test on unit to understand if it is a transport or velocity 153 ;IF it is a transport it is transofrmed in velocity 154 V1.arr = V1.arr / e3v_3d(/e1) * vmask() 155 V1.unit = 'm/s' 156 ENDIF 157 IF std_file2_V NE std_file1_V THEN BEGIN 158 V2 = read_ncdf(getenv('VAR2_V'), allrecords = allrec, direc = 't', filename = std_file2_V ) 159 ; old formulation 160 ; IF strlowcase(getenv('VAR2_V')) EQ 'vocetr_eff' OR strlowcase(getenv('VAR2_V')) EQ 'vozoeftr' THEN BEGIN 161 IF strlowcase(V2.unit) EQ 'm3/s' THEN BEGIN 162 V2.arr = V2.arr / e3v_3d(/e1) * vmask() 163 V2.unit = 'm/s' 164 ENDIF 165 ENDIF ELSE V2 = {arr:-1} 156 166 157 167 ; ice … … 177 187 ENDELSE 178 188 ; 189 Iage_1 = read_ncdf(getenv('VAR1_Iage'), allrecords = allrec, filename = std_file1_I ) 190 caldat, time, mm 191 march = where(mm EQ 3, cnt) 192 Iage_march_1 = {arr:1./float(cnt) * total(reform(Iage_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Iage_1.unit} 193 sept = where(mm EQ 9, cnt) 194 Iage_sept_1 = {arr:1./float(cnt) * total(reform(Iage_1.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Iage_1.unit} 195 undefine, Iage_1 196 ; 197 IF std_file2_I NE std_file1_I THEN BEGIN 198 Iage_2 = read_ncdf(getenv('VAR2_Iage'), allrecords = allrec, filename = std_file2_I ) 199 caldat, time, mm 200 march = where(mm EQ 3, cnt) 201 Iage_march_2 = {arr:1./float(cnt) * total(reform(Iage_2.arr[*, *, temporary(april)],nxt,nyt,cnt), 3), unit:Iage_2.unit} 202 sept = where(mm EQ 9, cnt) 203 Iage_sept_2 = {arr:1./float(cnt) * total(reform(Iage_2.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Iage_2.unit} 204 undefine, Iage_2 205 ENDIF ELSE BEGIN 206 Iage_march_2 = {arr:-1} 207 Iage_sept_2 = {arr:-1} 208 ENDELSE 209 ; 179 210 Ifra_1 = read_ncdf(getenv('VAR1_Ifrac'), allrecords = allrec, filename = std_file1_I ) 180 help, Ifra_1.arr, jpi,jpj181 211 caldat, time, mm 182 212 march = where(mm EQ 3, cnt) … … 199 229 ENDELSE 200 230 ; 201 ;; Sf commentato perché non funziona snow 202 ;; Sf Isnow_1 = read_ncdf(getenv('VAR1_Isnow'), allrecords = allrec, filename = std_file1_I ) 203 ;; Sf help, Isnow_1.arr, jpi,jpj 204 ;; Sf caldat, time, mm 205 ;; Sf april = where(mm EQ 4, cnt) 206 ;; Sf Isnow_april_1 = {arr:1./float(cnt) * total(reform(Isnow_1.arr[*, *, temporary(april)],nxt,nyt,cnt), 3), unit:Isnow_1.unit} 207 ;; Sf sept = where(mm EQ 9, cnt) 208 ;; Sf Isnow_sept_1 = {arr:1./float(cnt) * total(reform(Isnow_1.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Isnow_1.unit} 209 ;; Sf ;undefine, Isnow_1 210 ;; Sf ; 211 ;; Sf IF std_file2_I NE std_file1_I THEN BEGIN 212 ;; Sf Isnow_2 = read_ncdf(getenv('VAR2_Isnow'), allrecords = allrec, filename = std_file2_I ) 213 ;; Sf caldat, time, mm 214 ;; Sf april = where(mm EQ 4, cnt) 215 ;; Sf Isnow_april_2 = {arr:1./float(cnt) * total(reform(Isnow_2.arr[*, *, temporary(april)],nxt,nyt,cnt), 3), unit:Isnow_2.unit} 216 ;; Sf sept = where(mm EQ 9, cnt) 217 ;; Sf Isnow_sept_2 = {arr:1./float(cnt) * total(reform(Isnow_2.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Isnow_2.unit} 218 ;; Sf ; undefine, Isnow_2 219 ;; Sf ENDIF ELSE BEGIN 220 ;; Sf Isnow_april_2 = {arr:-1} 221 ;; Sf Isnow_sept_2 = {arr:-1} 222 ;; Sf ENDELSE 231 Isnow_1 = read_ncdf(getenv('VAR1_Isnow'), allrecords = allrec, filename = std_file1_I ) 232 caldat, time, mm 233 march = where(mm EQ 3, cnt) 234 Isnow_march_1 = {arr:1./float(cnt) * total(reform(Isnow_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Isnow_1.unit} 235 sept = where(mm EQ 9, cnt) 236 Isnow_sept_1 = {arr:1./float(cnt) * total(reform(Isnow_1.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Isnow_1.unit} 237 ;undefine, Isnow_1 238 ; 239 IF std_file2_I NE std_file1_I THEN BEGIN 240 Isnow_2 = read_ncdf(getenv('VAR2_Isnow'), allrecords = allrec, filename = std_file2_I ) 241 caldat, time, mm 242 march = where(mm EQ 3, cnt) 243 Isnow_march_2 = {arr:1./float(cnt) * total(reform(Isnow_2.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Isnow_2.unit} 244 sept = where(mm EQ 9, cnt) 245 Isnow_sept_2 = {arr:1./float(cnt) * total(reform(Isnow_2.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Isnow_2.unit} 246 ; undefine, Isnow_2 247 ENDIF ELSE BEGIN 248 Isnow_march_2 = {arr:-1} 249 Isnow_sept_2 = {arr:-1} 250 ENDELSE 251 ; 252 Isal_1 = read_ncdf(getenv('VAR1_Isal'), allrecords = allrec, filename = std_file1_I ) 253 Ifra_1 = read_ncdf(getenv('VAR1_Ifrac'), allrecords = allrec, filename = std_file1_I ) 254 msk = Ifra_1.arr gt 0.15 ; remove 0.15% for observations 255 help, Ifra_1 256 help, mask 257 caldat, time, mm 258 march = where(mm EQ 3, cnt) 259 febr = where(mm EQ 2, cnt) 260 ;SF Isal_march_1 = {arr:1./float(cnt) * total(reform(Isal_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Isal_1.unit} 261 Isal_1.arr = Isal_1.arr * msk 262 help, Isal_1 263 Isal_march_1 = {arr:1./float(cnt) * total(reform(Isal_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Isal_1.unit} 264 help, Isal_march_1 265 Isal_febr_1 = {arr:1./float(cnt) * total(reform(Isal_1.arr[*, *, temporary(febr)],nxt,nyt,cnt), 3), unit:Isal_1.unit} 266 sept = where(mm EQ 9, cnt) 267 Isal_sept_1 = {arr:1./float(cnt) * total(reform(Isal_1.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Isal_1.unit} 268 ;undefine, Isal_1 269 ; 270 Isal_2 = read_ncdf(getenv('VAR2_Isal'), allrecords = allrec, filename = std_file2_I ) 271 Ifra_2 = read_ncdf(getenv('VAR2_Ifrac'), allrecords = allrec, filename = std_file2_I ) 272 msk = Ifra_2.arr gt 0.15 ; remove 0.15% for observations 273 IF std_file2_I NE std_file1_I THEN BEGIN 274 Isal_2 = read_ncdf(getenv('VAR2_Isal'), allrecords = allrec, filename = std_file2_I ) 275 Isal_2.arr = Isal_2.arr * msk 276 caldat, time, mm 277 march = where(mm EQ 3, cnt) 278 febr = where(mm EQ 2, cnt) 279 Isal_march_2 = {arr:1./float(cnt) * total(reform(Isal_2.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Isal_2.unit} 280 Isal_febr_2 = {arr:1./float(cnt) * total(reform(Isal_2.arr[*, *, temporary(febr)],nxt,nyt,cnt), 3), unit:Isal_2.unit} 281 sept = where(mm EQ 9, cnt) 282 Isal_sept_2 = {arr:1./float(cnt) * total(reform(Isal_2.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Isal_2.unit} 283 ; undefine, Isal_2 284 ENDIF ELSE BEGIN 285 Isal_febr_2 = {arr:-1} 286 Isal_march_2 = {arr:-1} 287 Isal_sept_2 = {arr:-1} 288 ENDELSE 223 289 ; 224 290 jpt = 1 … … 256 322 cnt = cnt+1 & blabla = 'Erp salinity damping term' 257 323 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_erp, ERP1, ERP2, _extra = ex 258 324 ; 259 325 cnt = cnt+1 & blabla = 'Evaporation - Precipitation - Runoff term' 260 326 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_emp, EMP1, EMP2, _extra = ex 261 327 ; 262 328 cnt = cnt+1 & blabla = 'Net heat flux' 263 329 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_qnet, Q1, Q2, QNET, _extra = ex 264 330 ; 265 331 cnt = cnt+1 & blabla = 'Meridionnal Heat Transport' 266 332 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_mht, Q1.arr+geo, Q2.arr+geo, masknp, std_file_msksub, _extra = ex 267 333 ; 268 334 cnt = cnt+1 & blabla = 'Global Barotropic stream Function' 269 335 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_bsf, U1, U2, _extra = ex 270 336 ; 271 337 cnt = cnt+1 & blabla = 'mean Temperature diff with New Reynolds' 272 338 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_sst, T1, T2, TRey, _extra = ex 273 339 ; 274 340 cnt = cnt+1 & blabla = 'mean Salinity diff with Levitus' 275 341 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_sss, S1, S2, SLev, _extra = ex 276 342 ; 277 343 cnt = cnt+1 & blabla = 'Arctic mean Salinity diff with Levitus' 278 344 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_ArcSal, S1, SLev, _extra = ex 279 345 ; 280 346 cnt = cnt+1 & blabla = 'Arctic mean Salinity diff with Levitus and exp2' 281 347 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_ArcSal, S1, S2, SLev, _extra = ex 282 348 ; 283 349 cnt = cnt+1 & blabla = 'Arctic mean Salinity diff with Levitus at z=100 meters' 284 350 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_ArcSal, S1, SLev, /z100, _extra = ex 285 351 ; 286 352 cnt = cnt+1 & blabla = 'Arctic mean Salinity diff with Levitus and exp2 at z=100 meters' 287 353 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_ArcSal, S1, S2, SLev, /z100, _extra = ex 288 354 ; 289 355 cnt = cnt+1 & blabla = 'mean Temperature diff with Levitus at z=100 meters' 290 356 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_T100m, T1, T2, Tlev, _extra = ex 291 357 ; 292 358 cnt = cnt+1 & blabla = 'mean Salinity diff with Levitus at z=100 meters' 293 359 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_S100m, S1, S2, SLev, _extra = ex 294 360 ; 295 361 cnt = cnt+1 & blabla = 'Mixed layer depth' 296 362 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_mld, MLD1, MLD, _extra = ex 297 363 ; 298 364 cnt = cnt+1 & blabla = 'Mixed layer depth differences' 299 365 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_mld, MLD1, MLD2, MLD, _extra = ex 300 366 ; 301 367 cnt = cnt+1 & blabla = 'Zonal mean Mixed layer depth' 302 368 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_ZonMld, MLD1, MLD2, MLD, _extra = ex 303 369 ; 304 370 cnt = cnt+1 & blabla = 'Zonal mean Temperature diff with Levitus: Global' 305 371 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_zonal_T, T1, T2, TLev, _extra = ex 306 372 ; 307 373 cnt = cnt+1 & blabla = 'Zonal mean Temperature diff with Levitus: Atlantic' 308 374 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_zonal_T, T1, T2, TLev, SUBBASIN = 'Atl', _extra = ex 309 375 ; 310 376 cnt = cnt+1 & blabla = 'Zonal mean Temperature diff with Levitus: Indian' 311 377 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_zonal_T, T1, T2, TLev, SUBBASIN = 'Ind', _extra = ex 312 378 ; 313 379 cnt = cnt+1 & blabla = 'Zonal mean Temperature diff with Levitus: Pacific' 314 380 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_zonal_T, T1, T2, TLev, SUBBASIN = 'Pac', _extra = ex 315 381 ; 316 382 cnt = cnt+1 & blabla = 'Zonal mean Salinity diff with Levitus: Global' 317 383 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_zonal_S, S1, S2, SLev, _extra = ex 318 384 ; 319 385 cnt = cnt+1 & blabla = 'Zonal mean Salinity diff with Levitus: Atlantic' 320 386 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_zonal_S, S1, S2, SLev, SUBBASIN = 'Atl', _extra = ex 321 387 ; 322 388 cnt = cnt+1 & blabla = 'Zonal mean Salinity diff with Levitus: Indian' 323 389 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_zonal_S, S1, S2, SLev, SUBBASIN = 'Ind', _extra = ex 324 390 ; 325 391 cnt = cnt+1 & blabla = 'Zonal mean Salinity diff with Levitus: Pacific' 326 392 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_zonal_S, S1, S2, SLev, SUBBASIN = 'Pac', _extra = ex 327 393 ; 328 394 ; cnt = cnt+1 & blabla = 'Arctic Ice Thickness: MARCH' 329 395 ; IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_march_1, Ithi_march_2, /ARC, /MARCH, _extra = ex 330 396 ; 331 397 cnt = cnt+1 & blabla = 'Arctic Ice Thickness: APRIL' 332 398 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_april_1, Ithi_april_2, /ARC, /APRIL, _extra = ex 333 399 ; 334 400 cnt = cnt+1 & blabla = 'Arctic Ice Thickness: SEPT' 335 401 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_sept_1, Ithi_sept_2, /ARC, /SEPT, _extra = ex 336 402 ; 337 403 ; cnt = cnt+1 & blabla = 'Antarctic Ice Thickness: MARCH' 338 404 ; IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_march_1, Ithi_march_2, /ANT, /MARCH, _extra = ex 339 405 ; 340 406 cnt = cnt+1 & blabla = 'Antarctic Ice Thickness: APRIL' 341 407 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_april_1, Ithi_april_2, /ANT, /APRIL, _extra = ex 342 408 ; 343 409 cnt = cnt+1 & blabla = 'Antarctic Ice Thickness: SEPT' 344 410 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_sept_1, Ithi_sept_2, /ANT, /SEPT, _extra = ex 345 411 ; 412 cnt = cnt+1 & blabla = 'Arctic Ice Age: MARCH' 413 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceAge, Iage_march_1, Iage_march_2, /ARC, /MARCH, _extra = ex 414 ; 415 cnt = cnt+1 & blabla = 'Arctic Ice Age: SEPT' 416 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceAge, Iage_sept_1, Iage_sept_2, /ARC, /SEPT, _extra = ex 417 ; 418 cnt = cnt+1 & blabla = 'Antarctic Ice Age: MARCH' 419 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceAge, Iage_march_1, Iage_march_2, /ANT, /MARCH, _extra = ex 420 ; 421 cnt = cnt+1 & blabla = 'Antarctic Ice Age: SEPT' 422 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceAge, Iage_sept_1, Iage_sept_2, /ANT, /SEPT, _extra = ex 423 ; 346 424 cnt = cnt+1 & blabla = 'Arctic Ice Fraction: MARCH' 347 425 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceFrac, Ifra_march_1, Ifra_march_2, /ARC, /MARCH, _extra = ex 348 426 ; 349 427 cnt = cnt+1 & blabla = 'Arctic Ice Fraction: SEPT' 350 428 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceFrac, Ifra_sept_1, Ifra_sept_2, /ARC, /SEPT, _extra = ex 351 429 ; 352 430 cnt = cnt+1 & blabla = 'Antarctic Ice Fraction: MARCH' 353 431 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceFrac, Ifra_march_1, Ifra_march_2, /ANT, /MARCH, _extra = ex 354 432 ; 355 433 cnt = cnt+1 & blabla = 'Antarctic Ice Fraction: SEPT' 356 434 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceFrac, Ifra_sept_1, Ifra_sept_2, /ANT, /SEPT, _extra = ex 357 ;; SF non funziona, dice che deve essere una struttura e non e!!!! 358 ;;SF cnt = cnt+1 & blabla = 'Arctic SNOW Thickness: APRIL' 359 ;;SF IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_april_1, Isnow_april_2, /ARC, /APRIL, _extra = ex 360 ;; SF 361 ;; SF cnt = cnt+1 & blabla = 'Arctic SNOW Thickness: SEPT' 362 ;; SF IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_sept_1, Isnow_sept_2, /ARC, /SEPT, _extra = ex 363 ;; SF 364 ;; SF cnt = cnt+1 & blabla = 'Antarctic SNOW Thickness: APRIL' 365 ;; SF IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_april_1, Isnow_april_2, /ANT, /APRIL, _extra = ex 366 ;; SF 367 ;; SF cnt = cnt+1 & blabla = 'Antarctic SNOW Thickness: SEPT' 368 ;; SF IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_sept_1, Isnow_sept_2, /ANT, /SEPT, _extra = ex 369 ; 370 ; cnt = cnt+1 & blabla = 'Arctic Ice Extend: ' 371 ; IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceExt, Ifra_1, Ifra_2, masknp, /ARC, _extra = ex 372 373 ; cnt = cnt+1 & blabla = 'Antarctic Ice Extend: ' 374 ; IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceExt, Ifra_1, Ifra_2, masknp, /ANT, _extra = ex 375 435 ; 436 cnt = cnt+1 & blabla = 'Arctic SNOW Thickness: MARCH' 437 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_march_1, Isnow_march_2, /ARC, /MARCH, _extra = ex 438 ; 439 cnt = cnt+1 & blabla = 'Arctic SNOW Thickness: SEPT' 440 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_sept_1, Isnow_sept_2, /ARC, /SEPT, _extra = ex 441 ; 442 cnt = cnt+1 & blabla = 'Antarctic SNOW Thickness: MARCH' 443 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_march_1, Isnow_march_2, /ANT, /MARCH, _extra = ex 444 ; 445 cnt = cnt+1 & blabla = 'Antarctic SNOW Thickness: SEPT' 446 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_sept_1, Isnow_sept_2, /ANT, /SEPT, _extra = ex 447 ; 448 cnt = cnt+1 & blabla = 'Arctic Ice Salinity: MARCH' 449 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceSal, Isal_march_1, Isal_march_2, /ARC, /MARCH, _extra = ex 450 ; 451 cnt = cnt+1 & blabla = 'Arctic Ice Salinity: SEPT' 452 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceSal, Isal_sept_1, Isal_sept_2, /ARC, /SEPT, _extra = ex 453 ; 454 cnt = cnt+1 & blabla = 'Antarctic Ice Salinity: FEBRUARY' 455 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceSal, Isal_febr_1, Isal_febr_2, /ANT, /FEBR, _extra = ex 456 ; 457 cnt = cnt+1 & blabla = 'Antarctic Ice Salinity: SEPT' 458 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceSal, Isal_sept_1, Isal_sept_2, /ANT, /SEPT, _extra = ex 459 ; 376 460 cnt = cnt+1 & blabla = 'Meridional stream Function: Global (no Med)' 377 461 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_msf, V1, V2, SUBBASIN = 'GloNoMed', _extra = ex 378 462 ; 379 463 cnt = cnt+1 & blabla = 'Meridional stream Function: Atlantic' 380 464 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_msf, V1, V2, SUBBASIN = 'Atl', _extra = ex 381 465 ; 382 466 cnt = cnt+1 & blabla = 'Meridional stream Function: Indian' 383 467 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_msf, V1, V2, SUBBASIN = 'Ind', _extra = ex 384 468 ; 385 469 cnt = cnt+1 & blabla = 'Meridional stream Function: Indo-Pacific' 386 470 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_msf, V1, V2, SUBBASIN = 'IndoPac', _extra = ex 387 471 ; 388 472 cnt = cnt+1 & blabla = 'Equatorial Temperature' 389 473 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_EqT, T1, T2, Tlev, _extra = ex 390 474 ; 391 475 cnt = cnt+1 & blabla = 'Equatorial Salinity' 392 476 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_EqS, S1, S2, SLev, _extra = ex 393 477 ; 394 478 cnt = cnt+1 & blabla = 'Equatorial zonal velocity' 395 479 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_EqU, U1, U2, _extra = ex 396 480 ; 397 481 cnt = cnt+1 & blabla = 'Mediterranean salt tongue at depth=700' 398 482 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_Med_Sspread, S1, S2, SLev, 700, _extra = ex 399 483 ; 400 484 cnt = cnt+1 & blabla = 'Mediterranean salt tongue at depth=1000' 401 485 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_Med_Sspread, S1, S2, SLev, 1000, _extra = ex 402 486 ; 403 487 cnt = cnt+1 & blabla = 'Mediterranean water at lat=40°N' 404 488 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_Med_Sdepth, S1, S2, SLev, 40, _extra = ex 405 489 ; 406 490 cnt = cnt+1 & blabla = 'Mediterranean water at lat=38°N' 407 491 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_Med_Sdepth, S1, S2, SLev, 38, _extra = ex 408 492 ; 409 493 cnt = cnt+1 & blabla = 'Vertical Global mean T & S' 410 494 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_GlobMeanTS, T1, T2, TLev, S1, S2, SLev, _extra = ex 411 495 ; 412 496 IF n_elements(htmltxt) GT 1 THEN putfile, psdir+'std_plot_html_body.txt', htmltxt[1:*] 413 497 -
branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_all_last_year_5D.pro
r4320 r4514 203 203 ENDELSE 204 204 ; 205 Iage_1 = read_ncdf(getenv('VAR1_Iage'), allrecords = allrec, filename = std_file1_I ) 206 caldat, time, mm 207 march = where(mm EQ 3, cnt) 208 Iage_march_1 = {arr:1./float(cnt) * total(reform(Iage_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Iage_1.unit} 209 sept = where(mm EQ 9, cnt) 210 Iage_sept_1 = {arr:1./float(cnt) * total(reform(Iage_1.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Iage_1.unit} 211 undefine, Iage_1 212 ; 213 IF std_file2_I NE std_file1_I THEN BEGIN 214 Iage_2 = read_ncdf(getenv('VAR2_Iage'), allrecords = allrec, filename = std_file2_I ) 215 caldat, time, mm 216 march = where(mm EQ 3, cnt) 217 Iage_march_2 = {arr:1./float(cnt) * total(reform(Iage_2.arr[*, *, temporary(april)],nxt,nyt,cnt), 3), unit:Iage_2.unit} 218 sept = where(mm EQ 9, cnt) 219 Iage_sept_2 = {arr:1./float(cnt) * total(reform(Iage_2.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Iage_2.unit} 220 undefine, Iage_2 221 ENDIF ELSE BEGIN 222 Iage_march_2 = {arr:-1} 223 Iage_sept_2 = {arr:-1} 224 ENDELSE 225 ; 205 226 Ifra_1 = read_ncdf(getenv('VAR1_Ifrac'), allrecords = allrec, filename = std_file1_I ) 206 227 ;; help, Ifra_1.arr, jpi,jpj … … 225 246 ENDELSE 226 247 ; 227 ;; SF commentato perché non funziona snow 228 ;; SF Isnow_1 = read_ncdf(getenv('VAR1_Isnow'), allrecords = allrec, filename = std_file1_I ) 229 ;; SF help, Isnow_1.arr, jpi,jpj 230 ;; SF caldat, time, mm 231 ;; SF april = where(mm EQ 4, cnt) 232 ;; SF Isnow_april_1 = {arr:1./float(cnt) * total(reform(Isnow_1.arr[*, *, temporary(april)],nxt,nyt,cnt), 3), unit:Isnow_1.unit} 233 ;; SF sept = where(mm EQ 9, cnt) 234 ;; SF Isnow_sept_1 = {arr:1./float(cnt) * total(reform(Isnow_1.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Isnow_1.unit} 235 ;; SF ;undefine, Isnow_1 236 ;; SF ; 237 ;; SF IF std_file2_I NE std_file1_I THEN BEGIN 238 ;; SF Isnow_2 = read_ncdf(getenv('VAR2_Isnow'), allrecords = allrec, filename = std_file2_I ) 239 ;; SF caldat, time, mm 240 ;; SF april = where(mm EQ 4, cnt) 241 ;; SF Isnow_april_2 = {arr:1./float(cnt) * total(reform(Isnow_2.arr[*, *, temporary(april)],nxt,nyt,cnt), 3), unit:Isnow_2.unit} 242 ;; SF sept = where(mm EQ 9, cnt) 243 ;; SF Isnow_sept_2 = {arr:1./float(cnt) * total(reform(Isnow_2.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Isnow_2.unit} 244 ;; SF ; undefine, Isnow_2 245 ;; SF ENDIF ELSE BEGIN 246 ;; SF Isnow_april_2 = {arr:-1} 247 ;; SF Isnow_sept_2 = {arr:-1} 248 ;; SF ENDELSE 248 Isnow_1 = read_ncdf(getenv('VAR1_Isnow'), allrecords = allrec, filename = std_file1_I ) 249 caldat, time, mm 250 march = where(mm EQ 3, cnt) 251 Isnow_march_1 = {arr:1./float(cnt) * total(reform(Isnow_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Isnow_1.unit} 252 sept = where(mm EQ 9, cnt) 253 Isnow_sept_1 = {arr:1./float(cnt) * total(reform(Isnow_1.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Isnow_1.unit} 254 ;undefine, Isnow_1 255 ; 256 IF std_file2_I NE std_file1_I THEN BEGIN 257 Isnow_2 = read_ncdf(getenv('VAR2_Isnow'), allrecords = allrec, filename = std_file2_I ) 258 caldat, time, mm 259 march = where(mm EQ 3, cnt) 260 Isnow_march_2 = {arr:1./float(cnt) * total(reform(Isnow_2.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Isnow_2.unit} 261 sept = where(mm EQ 9, cnt) 262 Isnow_sept_2 = {arr:1./float(cnt) * total(reform(Isnow_2.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Isnow_2.unit} 263 ; undefine, Isnow_2 264 ENDIF ELSE BEGIN 265 Isnow_march_2 = {arr:-1} 266 Isnow_sept_2 = {arr:-1} 267 ENDELSE 268 ; 269 Isal_1 = read_ncdf(getenv('VAR1_Isal'), allrecords = allrec, filename = std_file1_I ) 270 Ifra_1 = read_ncdf(getenv('VAR1_Ifrac'), allrecords = allrec, filename = std_file1_I ) 271 msk = Ifra_1.arr gt 0.15 ; remove 0.15% for observations 272 help, Ifra_1 273 help, mask 274 caldat, time, mm 275 march = where(mm EQ 3, cnt) 276 febr = where(mm EQ 2, cnt) 277 ;SF Isal_march_1 = {arr:1./float(cnt) * total(reform(Isal_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Isal_1.unit} 278 Isal_1.arr = Isal_1.arr * msk 279 help, Isal_1 280 Isal_march_1 = {arr:1./float(cnt) * total(reform(Isal_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Isal_1.unit} 281 help, Isal_march_1 282 Isal_febr_1 = {arr:1./float(cnt) * total(reform(Isal_1.arr[*, *, temporary(febr)],nxt,nyt,cnt), 3), unit:Isal_1.unit} 283 sept = where(mm EQ 9, cnt) 284 Isal_sept_1 = {arr:1./float(cnt) * total(reform(Isal_1.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Isal_1.unit} 285 ;undefine, Isal_1 286 ; 287 Isal_2 = read_ncdf(getenv('VAR2_Isal'), allrecords = allrec, filename = std_file2_I ) 288 Ifra_2 = read_ncdf(getenv('VAR2_Ifrac'), allrecords = allrec, filename = std_file2_I ) 289 msk = Ifra_2.arr gt 0.15 ; remove 0.15% for observations 290 IF std_file2_I NE std_file1_I THEN BEGIN 291 Isal_2 = read_ncdf(getenv('VAR2_Isal'), allrecords = allrec, filename = std_file2_I ) 292 Isal_2.arr = Isal_2.arr * msk 293 caldat, time, mm 294 march = where(mm EQ 3, cnt) 295 febr = where(mm EQ 2, cnt) 296 Isal_march_2 = {arr:1./float(cnt) * total(reform(Isal_2.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Isal_2.unit} 297 Isal_febr_2 = {arr:1./float(cnt) * total(reform(Isal_2.arr[*, *, temporary(febr)],nxt,nyt,cnt), 3), unit:Isal_2.unit} 298 sept = where(mm EQ 9, cnt) 299 Isal_sept_2 = {arr:1./float(cnt) * total(reform(Isal_2.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Isal_2.unit} 300 ; undefine, Isal_2 301 ENDIF ELSE BEGIN 302 Isal_febr_2 = {arr:-1} 303 Isal_march_2 = {arr:-1} 304 Isal_sept_2 = {arr:-1} 305 ENDELSE 249 306 ; 250 307 jpt = 1 … … 282 339 cnt = cnt+1 & blabla = 'Erp salinity damping term' 283 340 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_erp, ERP1, ERP2, _extra = ex 284 341 ; 285 342 cnt = cnt+1 & blabla = 'Evaporation - Precipitation - Runoff term' 286 343 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_emp, EMP1, EMP2, _extra = ex 287 344 ; 288 345 cnt = cnt+1 & blabla = 'Net heat flux' 289 346 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_qnet, Q1, Q2, QNET, _extra = ex 290 347 ; 291 348 cnt = cnt+1 & blabla = 'Meridionnal Heat Transport' 292 349 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_mht, Q1.arr+geo, Q2.arr+geo, masknp, std_file_msksub, _extra = ex 293 350 ; 294 351 cnt = cnt+1 & blabla = 'Global Barotropic stream Function' 295 352 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_bsf, U1, U2, _extra = ex 296 353 ; 297 354 cnt = cnt+1 & blabla = 'mean Temperature diff with New Reynolds' 298 355 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_sst, T1, T2, TRey, _extra = ex 299 356 ; 300 357 cnt = cnt+1 & blabla = 'mean Salinity diff with Levitus' 301 358 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_sss, S1, S2, SLev, _extra = ex 302 359 ; 303 360 cnt = cnt+1 & blabla = 'Arctic mean Salinity diff with Levitus' 304 361 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_ArcSal, S1, SLev, _extra = ex 305 362 ; 306 363 cnt = cnt+1 & blabla = 'Arctic mean Salinity diff with Levitus and exp2' 307 364 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_ArcSal, S1, S2, SLev, _extra = ex 308 365 ; 309 366 cnt = cnt+1 & blabla = 'Arctic mean Salinity diff with Levitus at z=100 meters' 310 367 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_ArcSal, S1, SLev, /z100, _extra = ex 311 368 ; 312 369 cnt = cnt+1 & blabla = 'Arctic mean Salinity diff with Levitus and exp2 at z=100 meters' 313 370 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_ArcSal, S1, S2, SLev, /z100, _extra = ex 314 371 ; 315 372 cnt = cnt+1 & blabla = 'mean Temperature diff with Levitus at z=100 meters' 316 373 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_T100m, T1, T2, Tlev, _extra = ex 317 374 ; 318 375 cnt = cnt+1 & blabla = 'mean Salinity diff with Levitus at z=100 meters' 319 376 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_S100m, S1, S2, SLev, _extra = ex 320 321 ;;SF need monthly climatology to plot difference of MAX of MLD 322 cnt = cnt+1 & blabla = 'Mixed layer depth & difference with DeBoyer' 323 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_mld, MLD_1, MLD, _extra = ex 324 377 ; 378 cnt = cnt+1 & blabla = 'Mixed layer depth' 379 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_mld, MLD1, MLD, _extra = ex 380 ; 325 381 cnt = cnt+1 & blabla = 'Mixed layer depth differences' 326 382 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_mld, MLD_1, MLD_2, MLD, _extra = ex … … 339 395 cnt = cnt+1 & blabla = 'Zonal mean Temperature diff with Levitus: Global' 340 396 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_zonal_T, T1, T2, TLev, _extra = ex 341 397 ; 342 398 cnt = cnt+1 & blabla = 'Zonal mean Temperature diff with Levitus: Atlantic' 343 399 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_zonal_T, T1, T2, TLev, SUBBASIN = 'Atl', _extra = ex 344 400 ; 345 401 cnt = cnt+1 & blabla = 'Zonal mean Temperature diff with Levitus: Indian' 346 402 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_zonal_T, T1, T2, TLev, SUBBASIN = 'Ind', _extra = ex 347 403 ; 348 404 cnt = cnt+1 & blabla = 'Zonal mean Temperature diff with Levitus: Pacific' 349 405 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_zonal_T, T1, T2, TLev, SUBBASIN = 'Pac', _extra = ex 350 406 ; 351 407 cnt = cnt+1 & blabla = 'Zonal mean Salinity diff with Levitus: Global' 352 408 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_zonal_S, S1, S2, SLev, _extra = ex 353 409 ; 354 410 cnt = cnt+1 & blabla = 'Zonal mean Salinity diff with Levitus: Atlantic' 355 411 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_zonal_S, S1, S2, SLev, SUBBASIN = 'Atl', _extra = ex 356 412 ; 357 413 cnt = cnt+1 & blabla = 'Zonal mean Salinity diff with Levitus: Indian' 358 414 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_zonal_S, S1, S2, SLev, SUBBASIN = 'Ind', _extra = ex 359 415 ; 360 416 cnt = cnt+1 & blabla = 'Zonal mean Salinity diff with Levitus: Pacific' 361 417 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_zonal_S, S1, S2, SLev, SUBBASIN = 'Pac', _extra = ex 362 418 ; 363 419 ; cnt = cnt+1 & blabla = 'Arctic Ice Thickness: MARCH' 364 420 ; IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_march_1, Ithi_march_2, /ARC, /MARCH, _extra = ex 365 421 ; 366 422 cnt = cnt+1 & blabla = 'Arctic Ice Thickness: APRIL' 367 423 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_april_1, Ithi_april_2, /ARC, /APRIL, _extra = ex 368 424 ; 369 425 cnt = cnt+1 & blabla = 'Arctic Ice Thickness: SEPT' 370 426 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_sept_1, Ithi_sept_2, /ARC, /SEPT, _extra = ex 371 427 ; 372 428 ; cnt = cnt+1 & blabla = 'Antarctic Ice Thickness: MARCH' 373 429 ; IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_march_1, Ithi_march_2, /ANT, /MARCH, _extra = ex 374 430 ; 375 431 cnt = cnt+1 & blabla = 'Antarctic Ice Thickness: APRIL' 376 432 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_april_1, Ithi_april_2, /ANT, /APRIL, _extra = ex 377 433 ; 378 434 cnt = cnt+1 & blabla = 'Antarctic Ice Thickness: SEPT' 379 435 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_sept_1, Ithi_sept_2, /ANT, /SEPT, _extra = ex 380 436 ; 437 cnt = cnt+1 & blabla = 'Arctic Ice Age: MARCH' 438 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceAge, Iage_march_1, Iage_march_2, /ARC, /MARCH, _extra = ex 439 ; 440 cnt = cnt+1 & blabla = 'Arctic Ice Age: SEPT' 441 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceAge, Iage_sept_1, Iage_sept_2, /ARC, /SEPT, _extra = ex 442 ; 443 cnt = cnt+1 & blabla = 'Antarctic Ice Age: MARCH' 444 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceAge, Iage_march_1, Iage_march_2, /ANT, /MARCH, _extra = ex 445 ; 446 cnt = cnt+1 & blabla = 'Antarctic Ice Age: SEPT' 447 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceAge, Iage_sept_1, Iage_sept_2, /ANT, /SEPT, _extra = ex 448 ; 381 449 cnt = cnt+1 & blabla = 'Arctic Ice Fraction: MARCH' 382 450 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceFrac, Ifra_march_1, Ifra_march_2, /ARC, /MARCH, _extra = ex 383 451 ; 384 452 cnt = cnt+1 & blabla = 'Arctic Ice Fraction: SEPT' 385 453 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceFrac, Ifra_sept_1, Ifra_sept_2, /ARC, /SEPT, _extra = ex 386 454 ; 387 455 cnt = cnt+1 & blabla = 'Antarctic Ice Fraction: MARCH' 388 456 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceFrac, Ifra_march_1, Ifra_march_2, /ANT, /MARCH, _extra = ex 389 457 ; 390 458 cnt = cnt+1 & blabla = 'Antarctic Ice Fraction: SEPT' 391 459 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceFrac, Ifra_sept_1, Ifra_sept_2, /ANT, /SEPT, _extra = ex 392 393 ;; SF non funziona, dice che deve essere una struttura e non e!!!! 394 ;;SF cnt = cnt+1 & blabla = 'Arctic SNOW Thickness: APRIL' 395 ;;SF IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_april_1, Isnow_april_2, /ARC, /APRIL, _extra = ex 396 ;; SF 397 ;; SF cnt = cnt+1 & blabla = 'Arctic SNOW Thickness: SEPT' 398 ;; SF IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_sept_1, Isnow_sept_2, /ARC, /SEPT, _extra = ex 399 ;; SF 400 ;; SF cnt = cnt+1 & blabla = 'Antarctic SNOW Thickness: APRIL' 401 ;; SF IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_april_1, Isnow_april_2, /ANT, /APRIL, _extra = ex 402 ;; SF 403 ;; SF cnt = cnt+1 & blabla = 'Antarctic SNOW Thickness: SEPT' 404 ;; SF IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_sept_1, Isnow_sept_2, /ANT, /SEPT, _extra = ex 405 ; 406 ; cnt = cnt+1 & blabla = 'Arctic Ice Extend: ' 407 ; IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceExt, Ifra_1, Ifra_2, masknp, /ARC, _extra = ex 408 409 ; cnt = cnt+1 & blabla = 'Antarctic Ice Extend: ' 410 ; IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceExt, Ifra_1, Ifra_2, masknp, /ANT, _extra = ex 411 460 ; 461 cnt = cnt+1 & blabla = 'Arctic SNOW Thickness: MARCH' 462 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_march_1, Isnow_march_2, /ARC, /MARCH, _extra = ex 463 ; 464 cnt = cnt+1 & blabla = 'Arctic SNOW Thickness: SEPT' 465 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_sept_1, Isnow_sept_2, /ARC, /SEPT, _extra = ex 466 ; 467 cnt = cnt+1 & blabla = 'Antarctic SNOW Thickness: MARCH' 468 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_march_1, Isnow_march_2, /ANT, /MARCH, _extra = ex 469 ; 470 cnt = cnt+1 & blabla = 'Antarctic SNOW Thickness: SEPT' 471 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_sept_1, Isnow_sept_2, /ANT, /SEPT, _extra = ex 472 ; 473 cnt = cnt+1 & blabla = 'Arctic Ice Salinity: MARCH' 474 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceSal, Isal_march_1, Isal_march_2, /ARC, /MARCH, _extra = ex 475 ; 476 cnt = cnt+1 & blabla = 'Arctic Ice Salinity: SEPT' 477 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceSal, Isal_sept_1, Isal_sept_2, /ARC, /SEPT, _extra = ex 478 ; 479 cnt = cnt+1 & blabla = 'Antarctic Ice Salinity: FEBRUARY' 480 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceSal, Isal_febr_1, Isal_febr_2, /ANT, /FEBR, _extra = ex 481 ; 482 cnt = cnt+1 & blabla = 'Antarctic Ice Salinity: SEPT' 483 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceSal, Isal_sept_1, Isal_sept_2, /ANT, /SEPT, _extra = ex 484 ; 412 485 cnt = cnt+1 & blabla = 'Meridional stream Function: Global (no Med)' 413 486 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_msf, V1, V2, SUBBASIN = 'GloNoMed', _extra = ex 414 487 ; 415 488 cnt = cnt+1 & blabla = 'Meridional stream Function: Atlantic' 416 489 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_msf, V1, V2, SUBBASIN = 'Atl', _extra = ex 417 490 ; 418 491 cnt = cnt+1 & blabla = 'Meridional stream Function: Indian' 419 492 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_msf, V1, V2, SUBBASIN = 'Ind', _extra = ex 420 493 ; 421 494 cnt = cnt+1 & blabla = 'Meridional stream Function: Indo-Pacific' 422 495 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_msf, V1, V2, SUBBASIN = 'IndoPac', _extra = ex 423 496 ; 424 497 cnt = cnt+1 & blabla = 'Equatorial Temperature' 425 498 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_EqT, T1, T2, Tlev, _extra = ex 426 499 ; 427 500 cnt = cnt+1 & blabla = 'Equatorial Salinity' 428 501 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_EqS, S1, S2, SLev, _extra = ex 429 502 ; 430 503 cnt = cnt+1 & blabla = 'Equatorial zonal velocity' 431 504 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_EqU, U1, U2, _extra = ex 432 505 ; 433 506 cnt = cnt+1 & blabla = 'Mediterranean salt tongue at depth=700' 434 507 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_Med_Sspread, S1, S2, SLev, 700, _extra = ex 435 508 ; 436 509 cnt = cnt+1 & blabla = 'Mediterranean salt tongue at depth=1000' 437 510 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_Med_Sspread, S1, S2, SLev, 1000, _extra = ex 438 511 ; 439 512 cnt = cnt+1 & blabla = 'Mediterranean water at lat=40°N' 440 513 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_Med_Sdepth, S1, S2, SLev, 40, _extra = ex 441 514 ; 442 515 cnt = cnt+1 & blabla = 'Mediterranean water at lat=38°N' 443 516 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_Med_Sdepth, S1, S2, SLev, 38, _extra = ex 444 517 ; 445 518 cnt = cnt+1 & blabla = 'Vertical Global mean T & S' 446 519 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_GlobMeanTS, T1, T2, TLev, S1, S2, SLev, _extra = ex 447 520 ; 448 521 IF n_elements(htmltxt) GT 1 THEN putfile, psdir+'std_plot_html_body.txt', htmltxt[1:*] 449 522 -
branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_vardef.sh
r4507 r4514 69 69 VAR1_Ithick=sithic 70 70 VAR1_Ifrac=siconc 71 VAR1_Isnow=sn dept71 VAR1_Isnow=snvolu 72 72 VAR1_Isal=sisali 73 VAR1_Iage=siages 73 74 # 74 75 #===================== EXP2 ===================== … … 88 89 VAR2_Ithick=sithic 89 90 VAR2_Ifrac=siconc 90 VAR2_Isnow=sn dept91 VAR2_Isnow=snvolu 91 92 VAR2_Isal=sisali 93 VAR2_Iage=siages 92 94 # 93 95 ######################### Export Variables ############################### … … 113 115 export VAR1_T VAR1_S VAR1_QNET VAR1_ERP VAR1_EMP VAR1_MLD 114 116 export VAR1_U VAR1_V 115 export VAR1_Ithick VAR1_Ifrac VAR1_Isnow VAR1_Isal 117 export VAR1_Ithick VAR1_Ifrac VAR1_Isnow VAR1_Isal VAR1_Iage 116 118 #===================== EXP2 ===================== 117 119 export FILE2_T FILE2_U FILE2_V FILE2_I 118 120 export VAR2_T VAR2_S VAR2_QNET VAR2_ERP VAR2_EMP VAR2_MLD 119 121 export VAR2_U VAR2_V 120 export VAR2_Ithick VAR2_Ifrac VAR2_Isnow VAR2_Isal 122 export VAR2_Ithick VAR2_Ifrac VAR2_Isnow VAR2_Isal VAR2_Iage 121 123 # -
branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_vardef.sh.last_year_5D
r4507 r4514 26 26 SAXO_DIR=/Users/sflod/SAXO_DIR 27 27 # 28 DIR_DATA= /Users/sflod/idl_PLOTS/DATA_STORE/RUN_CLIMATO/lim3_ada # path of data in NetCDF format29 DIR_CLIMATO= /Users/sflod/idl_PLOTS/CLIMATOLOGIES # path of climatological data30 DIR_MASK= /Users/sflod/idl_PLOTS/MASK # path of mask files (ex: subbasins)28 DIR_DATA=~/idl_PLOTS/DATA_STORE/RUN_CLIMATO/lim3_ada # path of data in NetCDF format 29 DIR_CLIMATO=~/idl_PLOTS/CLIMATOLOGIES # path of climatological data 30 DIR_MASK=~//idl_PLOTS/MASK # path of mask files (ex: subbasins) 31 31 # 32 32 #===================== Model GRID ===================== … … 55 55 #===================== EXP1 ===================== 56 56 # 57 #FILE1_T=l3fwb178_24910101_25001231_1Y_grid_T.nc # exp1 grid T input file58 #FILE1_U=l3fwb178_24910101_25001231_1Y_grid_U.nc # exp1 grid U input file59 #FILE1_V=l3fwb178_24910101_25001231_1Y_grid_V.nc # exp1 grid V input file60 #FILE1_I=l3fwb178_24910101_25001231_1M_icemod.nc # exp1 ice input file61 57 FILE1_T=l3fwb178_25010101_25011231_5D_grid_T.nc # exp1 grid T input file 62 58 FILE1_U=l3fwb178_25010101_25011231_5D_grid_U.nc # exp1 grid U input file 63 59 FILE1_V=l3fwb178_25010101_25011231_5D_grid_V.nc # exp1 grid V input file 64 60 FILE1_I=l3fwb178_25010101_25011231_5D_icemod.nc # exp1 ice input file 65 VAR1_T=votemper 66 VAR1_S=vosaline 61 VAR1_T=votemper 62 VAR1_S=vosaline 67 63 VAR1_QNET=sohflxdo 68 VAR1_ERP=sowafldp 69 VAR1_EMP=sowaflup 70 #VAR1_MLD=somxl010 64 VAR1_ERP=sowafldp 65 VAR1_EMP=sowaflup 71 66 VAR1_MLD=max5dmld 72 #VAR1_U=uocetr_eff73 #VAR1_U=uoce74 #VAR1_U=vozocrtx75 67 VAR1_U=vozoeftr 76 #VAR1_V=vocetr_eff77 #VAR1_V=voce78 68 VAR1_V=vomeeftr 79 #VAR1_V=vomecrty 80 VAR1_Ithick=iicethic 69 VAR1_Ithick=iicethic 81 70 VAR1_Ifrac=iiceconc 82 VAR1_Isnow=isnow thi71 VAR1_Isnow=isnowvol 83 72 VAR1_Isal=iicesali 73 VAR1_Iage=iiceages 84 74 # 85 75 #===================== EXP2 ===================== 86 76 # 87 #FILE2_T=l3fwb178_24910101_25001231_1Y_grid_T.nc # exp1 grid T input file88 #FILE2_U=l3fwb178_24910101_25001231_1Y_grid_U.nc # exp1 grid U input file89 #FILE2_V=l3fwb178_24910101_25001231_1Y_grid_V.nc # exp1 grid V input file90 #FILE2_I=l3fwb178_24910101_25001231_1M_icemod.nc # exp1 ice input file91 77 FILE2_T=l3fwb178_25010101_25011231_5D_grid_T.nc # exp1 grid T input file 92 78 FILE2_U=l3fwb178_25010101_25011231_5D_grid_U.nc # exp1 grid U input file 93 79 FILE2_V=l3fwb178_25010101_25011231_5D_grid_V.nc # exp1 grid V input file 94 80 FILE2_I=l3fwb178_25010101_25011231_5D_icemod.nc # exp1 ice input file 95 #FILE2_T=l378sal6_25010101_25011231_5D_grid_T.nc # exp1 grid T input file96 #FILE2_U=l378sal6_25010101_25011231_5D_grid_U.nc # exp1 grid U input file97 #FILE2_V=l378sal6_25010101_25011231_5D_grid_V.nc # exp1 grid V input file98 #FILE2_I=l378sal6_25010101_25011231_5D_icemod.nc # exp1 ice input file99 81 VAR2_T=votemper 100 82 VAR2_S=vosaline … … 102 84 VAR2_ERP=sowafldp 103 85 VAR2_EMP=sowaflup 104 #VAR2_MLD=somxl010105 86 VAR2_MLD=max5dmld 106 #VAR2_U=uoce_eff107 #VAR2_U=uoce108 87 VAR2_U=vozoeftr 109 #VAR2_V=vocetr_eff110 #VAR2_V=voce111 88 VAR2_V=vomeeftr 112 89 VAR2_Ithick=iicethic … … 114 91 VAR2_Isnow=isnowthi 115 92 VAR2_Isal=iicesali 93 VAR2_Iage=iiceages 116 94 # 117 95 ######################### Export Variables ############################### … … 137 115 export VAR1_T VAR1_S VAR1_QNET VAR1_ERP VAR1_EMP VAR1_MLD 138 116 export VAR1_U VAR1_V 139 export VAR1_Ithick VAR1_Ifrac VAR1_Isnow VAR1_Isal 117 export VAR1_Ithick VAR1_Ifrac VAR1_Isnow VAR1_Isal VAR1_Iage 140 118 #===================== EXP2 ===================== 141 119 export FILE2_T FILE2_U FILE2_V FILE2_I 142 120 export VAR2_T VAR2_S VAR2_QNET VAR2_ERP VAR2_EMP VAR2_MLD 143 121 export VAR2_U VAR2_V 144 export VAR2_Ithick VAR2_Ifrac VAR2_Isnow VAR2_Isal 122 export VAR2_Ithick VAR2_Ifrac VAR2_Isnow VAR2_Isal VAR2_Iage 145 123 # -
branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_vardef.sh.new_names
r4507 r4514 69 69 VAR1_Ithick=sithic 70 70 VAR1_Ifrac=siconc 71 VAR1_Isnow=sn dept71 VAR1_Isnow=snvolu 72 72 VAR1_Isal=sisali 73 VAR1_Iage=siages 73 74 # 74 75 #===================== EXP2 ===================== … … 88 89 VAR2_Ithick=sithic 89 90 VAR2_Ifrac=siconc 90 VAR2_Isnow=sn dept91 VAR2_Isnow=snvolu 91 92 VAR2_Isal=sisali 93 VAR2_Iage=siages 92 94 # 93 95 ######################### Export Variables ############################### … … 113 115 export VAR1_T VAR1_S VAR1_QNET VAR1_ERP VAR1_EMP VAR1_MLD 114 116 export VAR1_U VAR1_V 115 export VAR1_Ithick VAR1_Ifrac VAR1_Isnow VAR1_Isal 117 export VAR1_Ithick VAR1_Ifrac VAR1_Isnow VAR1_Isal VAR1_Iage 116 118 #===================== EXP2 ===================== 117 119 export FILE2_T FILE2_U FILE2_V FILE2_I 118 120 export VAR2_T VAR2_S VAR2_QNET VAR2_ERP VAR2_EMP VAR2_MLD 119 121 export VAR2_U VAR2_V 120 export VAR2_Ithick VAR2_Ifrac VAR2_Isnow VAR2_Isal 122 export VAR2_Ithick VAR2_Ifrac VAR2_Isnow VAR2_Isal VAR2_Iage 121 123 # -
branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_vardef.sh_example1
r4507 r4514 26 26 SAXO_DIR=/Users/sflod/SAXO_DIR 27 27 # 28 DIR_DATA= /Users/sflod/idl_PLOTS/DATA_STORE/RUN_CLIMATO/lim3_ada# path of data in NetCDF format29 DIR_CLIMATO= /Users/sflod/idl_PLOTS/CLIMATOLOGIES # path of climatological data30 DIR_MASK= /Users/sflod/idl_PLOTS/MASK # path of mask files (ex: subbasins)28 DIR_DATA=~/idl_PLOTS/DATA_STORE/RUN_CLIMATO/lim3_ada # path of data in NetCDF format 29 DIR_CLIMATO=~/idl_PLOTS/CLIMATOLOGIES # path of climatological data 30 DIR_MASK=~//idl_PLOTS/MASK # path of mask files (ex: subbasins) 31 31 # 32 32 #===================== Model GRID ===================== … … 55 55 #===================== EXP1 ===================== 56 56 # 57 FILE1_T= tr33beta_00910101_01001231_1Y_grid_T.nc # exp1 grid T input file58 FILE1_U= tr33beta_00910101_01001231_1Y_grid_U.nc # exp1 grid U input file59 FILE1_V= tr33beta_00910101_01001231_1Y_grid_V.nc # exp1 grid V input file60 FILE1_I= tr33beta_00910101_01001231_1M_icemod.nc # exp1 ice input file61 VAR1_T= thetao62 VAR1_S= so63 VAR1_QNET= qt64 VAR1_ERP= wfcorr65 VAR1_EMP= wfo66 VAR1_MLD= mldr10_167 VAR1_U= uocetr_eff68 VAR1_V=vo cetr_eff69 VAR1_Ithick= sit70 VAR1_Ifrac= sic71 VAR1_Isnow=isnow thi57 FILE1_T=500yfwb0_24910101_25001231_1Y_grid_T.nc # exp1 grid T input file 58 FILE1_U=500yfwb0_24910101_25001231_1Y_grid_U.nc # exp1 grid U input file 59 FILE1_V=500yfwb0_24910101_25001231_1Y_grid_V.nc # exp1 grid V input file 60 FILE1_I=500yfwb0_24910101_25001231_1M_icemod.nc # exp1 ice input file 61 VAR1_T=votemper 62 VAR1_S=vosaline 63 VAR1_QNET=sohflxdo 64 VAR1_ERP=sowafldp 65 VAR1_EMP=sowaflup 66 VAR1_MLD=somxl010 67 VAR1_U=vozoeftr 68 VAR1_V=vomeeftr 69 VAR1_Ithick=iicethic 70 VAR1_Ifrac=iiceconc 71 VAR1_Isnow=isnowvol 72 72 VAR1_Isal=iicesali 73 VAR1_Iage=iiceages 73 74 # 74 75 #===================== EXP2 ===================== … … 90 91 VAR2_Isnow=isnowthi 91 92 VAR2_Isal=iicesali 93 VAR2_Iage=iiceages 92 94 # 93 95 ######################### Export Variables ############################### … … 113 115 export VAR1_T VAR1_S VAR1_QNET VAR1_ERP VAR1_EMP VAR1_MLD 114 116 export VAR1_U VAR1_V 115 export VAR1_Ithick VAR1_Ifrac VAR1_Isnow VAR1_Isal 117 export VAR1_Ithick VAR1_Ifrac VAR1_Isnow VAR1_Isal VAR1_Iage 116 118 #===================== EXP2 ===================== 117 119 export FILE2_T FILE2_U FILE2_V FILE2_I 118 120 export VAR2_T VAR2_S VAR2_QNET VAR2_ERP VAR2_EMP VAR2_MLD 119 121 export VAR2_U VAR2_V 120 export VAR2_Ithick VAR2_Ifrac VAR2_Isnow VAR2_Isal 122 export VAR2_Ithick VAR2_Ifrac VAR2_Isnow VAR2_Isal VAR2_Iage 121 123 # -
branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_vardef.sh_example2
r4507 r4514 21 21 #idl_command=idl71 22 22 # 23 PS_DIR=$( pwd )/ 500yfwb0_ps_plot_500y24 PDF_DIR=$( pwd )/ 500yfwb0_pdf_plot_500y23 PS_DIR=$( pwd )/ps_plot 24 PDF_DIR=$( pwd )/pdf_plot 25 25 HTML_DIR=$( pwd )/html_plot 26 26 SAXO_DIR=/Users/sflod/SAXO_DIR 27 27 # 28 DIR_DATA= /Users/sflod/idl_PLOTS/DATA_STORE/RUN_CLIMATO/lim3_ada # path of data in NetCDF format29 DIR_CLIMATO= /Users/sflod/idl_PLOTS/CLIMATOLOGIES # path of climatological data30 DIR_MASK= /Users/sflod/idl_PLOTS/MASK # path of mask files (ex: subbasins)28 DIR_DATA=~/idl_PLOTS/DATA_STORE/RUN_CLIMATO/lim3_ada # path of data in NetCDF format 29 DIR_CLIMATO=~/idl_PLOTS/CLIMATOLOGIES # path of climatological data 30 DIR_MASK=~//idl_PLOTS/MASK # path of mask files (ex: subbasins) 31 31 # 32 32 #===================== Model GRID ===================== … … 55 55 #===================== EXP1 ===================== 56 56 # 57 FILE1_T= tag33_20910101_21001231_1Y_grid_T.nc # exp1 grid T input file58 FILE1_U= tag33_20910101_21001231_1Y_grid_U.nc # exp1 grid U input file59 FILE1_V= tag33_20910101_21001231_1Y_grid_V.nc # exp1 grid V input file60 FILE1_I= tag33_20910101_21001231_1M_icemod.nc # exp1 ice input file61 VAR1_T= thetao62 VAR1_S= so63 VAR1_QNET= qt64 VAR1_ERP= wfcorr65 VAR1_EMP= wfo66 VAR1_MLD= mldr10_167 VAR1_U= uocetr_eff68 VAR1_V=vo cetr_eff69 VAR1_Ithick= sit70 VAR1_Ifrac= sic71 VAR1_Isnow=isnow thi57 FILE1_T=500yfwb0_24910101_25001231_1Y_grid_T.nc # exp1 grid T input file 58 FILE1_U=500yfwb0_24910101_25001231_1Y_grid_U.nc # exp1 grid U input file 59 FILE1_V=500yfwb0_24910101_25001231_1Y_grid_V.nc # exp1 grid V input file 60 FILE1_I=500yfwb0_24910101_25001231_1M_icemod.nc # exp1 ice input file 61 VAR1_T=votemper 62 VAR1_S=vosaline 63 VAR1_QNET=sohflxdo 64 VAR1_ERP=sowafldp 65 VAR1_EMP=sowaflup 66 VAR1_MLD=somxl010 67 VAR1_U=vozoeftr 68 VAR1_V=vomeeftr 69 VAR1_Ithick=iicethic 70 VAR1_Ifrac=iiceconc 71 VAR1_Isnow=isnowvol 72 72 VAR1_Isal=iicesali 73 VAR1_Iage=iiceages 73 74 # 74 75 #===================== EXP2 ===================== 75 76 # 76 FILE2_T=tmx_2000_00910101_01001231_1Y_grid_T.nc # exp1 grid T input file 77 FILE2_U=tmx_2000_00910101_01001231_1Y_grid_U.nc # exp1 grid U input file 78 FILE2_V=tmx_2000_00910101_01001231_1Y_grid_V.nc # exp1 grid V input file 79 FILE2_I=tmx_2000_00910101_01001231_1M_icemod.nc # exp1 ice input file 80 VAR2_T=thetao 81 VAR2_S=so 82 VAR2_QNET=qns+qsr 83 VAR2_ERP=wfcorr 84 VAR2_EMP=wfo 85 VAR2_MLD=mldr10_1 86 VAR2_U=uoce_eff 87 VAR2_V=voce_eff 88 VAR2_Ithick=sit 89 VAR2_Ifrac=sic 90 VAR2_Isnow=isnowthi 91 VAR2_Isal=iicesali 77 FILE1_T=500yfwb0_24910101_25001231_1Y_grid_T.nc # exp1 grid T input file 78 FILE1_U=500yfwb0_24910101_25001231_1Y_grid_U.nc # exp1 grid U input file 79 FILE1_V=500yfwb0_24910101_25001231_1Y_grid_V.nc # exp1 grid V input file 80 FILE1_I=500yfwb0_24910101_25001231_1M_icemod.nc # exp1 ice input file 81 VAR1_T=votemper 82 VAR1_S=vosaline 83 VAR1_QNET=sohflxdo 84 VAR1_ERP=sowafldp 85 VAR1_EMP=sowaflup 86 VAR1_MLD=somxl010 87 VAR1_U=vozoeftr 88 VAR1_V=vomeeftr 89 VAR1_Ithick=iicethic 90 VAR1_Ifrac=iiceconc 91 VAR1_Isnow=isnowvol 92 VAR1_Isal=iicesali 93 VAR1_Iage=iiceages 92 94 # 93 95 ######################### Export Variables ############################### … … 113 115 export VAR1_T VAR1_S VAR1_QNET VAR1_ERP VAR1_EMP VAR1_MLD 114 116 export VAR1_U VAR1_V 115 export VAR1_Ithick VAR1_Ifrac VAR1_Isnow VAR1_Isal 117 export VAR1_Ithick VAR1_Ifrac VAR1_Isnow VAR1_Isal VAR1_Iage 116 118 #===================== EXP2 ===================== 117 119 export FILE2_T FILE2_U FILE2_V FILE2_I 118 120 export VAR2_T VAR2_S VAR2_QNET VAR2_ERP VAR2_EMP VAR2_MLD 119 121 export VAR2_U VAR2_V 120 export VAR2_Ithick VAR2_Ifrac VAR2_Isnow VAR2_Isal 122 export VAR2_Ithick VAR2_Ifrac VAR2_Isnow VAR2_Isal VAR2_Iage 121 123 # -
branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_ICE_FRAM.pro
r4504 r4514 49 49 ; 50 50 ;OBSERVATIONS : mean seasonal cycle/month 51 vol_obs = [0.261625, 0.230750, 0.325375, 0.252000, 0.172500, 0.0805000, 0.0805000, 0.0805000, 0.0805000, 0.176500, 0.148500, 0.235000] 51 ; vol_obs = [0.261625, 0.230750, 0.325375, 0.252000, 0.172500, 0.0805000, 0.0805000, 0.0805000, 0.0805000, 0.176500, 0.148500, 0.235000] 52 vol_obs = [261.625, 230.750, 325.375, 252.000, 172.500, 80.5000, 80.5000, 80.5000, 80.5000, 176.500, 148.500, 235.000] 52 53 area_obs = [0.103292, 0.0997500, 0.107625, 0.0944167, 0.0612083, 0.0262500, 0.0262500, 0.0262500, 0.0262500, 0.0843750, 0.0914583, 0.104083] 53 54 ; … … 80 81 ICE_vol_export = reform(ICE_vol_export, 12, nyr) 81 82 ICE_vol_export = total(ICE_vol_export,2)/nyr 82 ; ICE_vol_export = {arr:ICE_vol_export * 1.e-12 * 86400 * 365 , unit : '10^3 Km^3/year'} ; annual mean 83 ICE_vol_export = {arr:ICE_vol_export * 1.e-12 * 86400 * 30 , unit : '10^3 Km^3/month'} ; monthly mean 84 85 ; 86 ;title = 'Northern Hemisphere'+'!C'+prefix+' '+d1_d2+'!C'+'Global Annual Mean Ice Volume (Black SOLID simulation)' 83 ; ICE_vol_export = {arr:ICE_vol_export * 1.e-9 * 86400 * 365 , unit : '10^3 Km^3/year'} ; annual mean 84 ICE_vol_export = {arr:ICE_vol_export * 1.e-9 * 86400 * 30 , unit : '10^3 Km^3/month'} ; monthly mean 85 86 ; 87 87 title = 'Fram Strait Areal Export'+'!C'+prefix+' '+d1_d2 88 88 jpt=12 … … 95 95 tot_area_expo = total(ICE_area_export.arr) 96 96 ; 97 xyouts, julday( 6,15,1900), 0.48, 'Tot. Annual Export OBS = 0.851 million Km2', ALIGN = 0, CHARTHICK = 2, CHARSIZE=1, COLOR=298 xyouts, julday( 6,15,1900), 0.46, 'Tot. Annual Export Model = '+strtrim(tot_area_expo, 1)+' million Km2', ALIGN = 0, CHARTHICK = 2, CHARSIZE=1, COLOR=299 xyouts, julday(7,15,1900), 0.4 2, 'Data from Kwok et al.(2004), 1979-2002 ', ALIGN = 0, CHARTHICK = 2, CHARSIZE=0.8, COLOR=2100 ; title ='Southern Hemisphere' +'!C'+prefix+' '+d1_d2+' - '+'!C'+'Global Annual Mean Ice Volume (Black SOLID simulation)'97 xyouts, julday(5,15,1900), 0.49, 'Tot. Annual Export OBS = 0.851 million Km2', ALIGN = 0, CHARTHICK = 2, CHARSIZE=1, COLOR=2 98 xyouts, julday(5,15,1900), 0.46, 'Tot. Annual Export Model = '+strtrim(tot_area_expo, 1)+' million Km2', ALIGN = 0, CHARTHICK = 2, CHARSIZE=1, COLOR=2 99 xyouts, julday(7,15,1900), 0.41, 'Data from Kwok et al.(2004), 1979-2002 ', ALIGN = 0, CHARTHICK = 2, CHARSIZE=0.8, COLOR=2 100 ; 101 101 title = 'Fram Strait Volume Export'+'!C'+prefix+' '+d1_d2 102 pltt, ICE_vol_export, 't', /REMPLI, MIN = 0., MAX = .6, /NOERASE, XGRIDSTYLE = 1 $102 pltt, ICE_vol_export, 't', /REMPLI, MIN = 60., MAX = 500. , /NOERASE, XGRIDSTYLE = 1 $ 103 103 , small = [1, 2, 2], YTITLE = '10^3 Km^3/month', TITLE = title, DATE_FORMAT = '%M', _extra = ex 104 104 pltt, vol_obs, 't', /REMPLI, /NOERASE, psym = 4, THICK = 4 $ ; light blue … … 107 107 tot_vol_expo = total(ICE_vol_export.arr) 108 108 tot_vol_expo_Sv = tot_vol_expo * 1.e06 * 1/86400 * 1/365 ; annual mean in Sverdrup 109 xyouts, julday( 6,15,1900), 0.48, 'Tot. Annual Export OBS = 2124Km3/year', ALIGN = 0, CHARTHICK = 2, CHARSIZE=1, COLOR=2110 xyouts, julday( 6,15,1900), 0.46, 'Tot. Annual Export Model = '+strtrim(tot_vol_expo, 1)+'Km3/year', ALIGN = 0, CHARTHICK = 2, CHARSIZE=1, COLOR=2111 xyouts, julday(9,15,1900), 0.44, 'in Sv = '+strtrim(tot_vol_expo_Sv, 1)+' Sv', ALIGN = 0, CHARTHICK = 2, CHARSIZE=1, COLOR=2112 xyouts, julday(7,15,1900), 0.38, 'Data from Kwok et al.(2004), 1992-1998 ', ALIGN = 0, CHARTHICK = 2, CHARSIZE=0.8, COLOR=2109 xyouts, julday(5,15,1900), 490, 'Tot. Annual Export OBS = 2124 10^3 Km3/year', ALIGN = 0, CHARTHICK = 2, CHARSIZE=1, COLOR=2 110 xyouts, julday(5,15,1900), 460, 'Tot. Annual Export Model = '+strtrim(tot_vol_expo, 1)+' 10^3 Km3/year', ALIGN = 0, CHARTHICK = 2, CHARSIZE=1, COLOR=2 111 xyouts, julday(9,15,1900), 440, 'in Sv = '+strtrim(tot_vol_expo_Sv, 1)+' Sv', ALIGN = 0, CHARTHICK = 2, CHARSIZE=1, COLOR=2 112 xyouts, julday(7,15,1900), 360, 'Data from Kwok et al.(2004), 1992-1998 ', ALIGN = 0, CHARTHICK = 2, CHARSIZE=0.8, COLOR=2 113 113 114 114 htmltxt = [ htmltxt, '<hr>'+blabla, '<br><img width="80%" src='+filename+'.png /> ' ] … … 147 147 ICE_vol_export_2 = total(ICE_vol_export_2,2)/nyr 148 148 ; ICE_vol_export_2 = {arr:ICE_vol_export_2 * 1.e-12 * 86400 * 365 , unit : '10^3 Km^3/year'} ; annual mean 149 ICE_vol_export_2 = {arr:ICE_vol_export_2 * 1.e- 12* 86400 * 30 , unit : '10^3 Km^3/month'} ; monthly mean149 ICE_vol_export_2 = {arr:ICE_vol_export_2 * 1.e-9 * 86400 * 30 , unit : '10^3 Km^3/month'} ; monthly mean 150 150 ; 151 151 ; … … 165 165 tot_area_expo_2 = total(ICE_area_export_2.arr) 166 166 ; 167 xyouts, julday(6,15,1900), 0.48, 'Tot. Annual Export OBS = 0.851 million Km2', ALIGN = 0, CHARTHICK = 2, CHARSIZE=1, COLOR=2 168 xyouts, julday(6,15,1900), 0.46, 'Tot. Annual Export Model 1= '+strtrim(tot_area_expo, 1)+' million Km2', ALIGN = 0, CHARTHICK = 2, CHARSIZE=1, COLOR=2 169 xyouts, julday(6,15,1900), 0.44, 'Tot. Annual Export Model 2= '+strtrim(tot_area_expo_2, 1)+' million Km2', ALIGN = 0, CHARTHICK = 2, CHARSIZE=1, COLOR=2 170 xyouts, julday(7,15,1900), 0.42, 'Data from Kwok et al.(2004), 1979-2002 ', ALIGN = 0, CHARTHICK = 2, CHARSIZE=0.8, COLOR=2 167 xyouts, julday(5,15,1900), 0.49, 'Tot. Annual Export OBS = 0.851 million Km2', ALIGN = 0, CHARTHICK = 2, CHARSIZE=1, COLOR=2 168 xyouts, julday(5,15,1900), 0.46, 'Tot. Annual Export Model 1= '+strtrim(tot_area_expo, 1)+' million Km2', ALIGN = 0, CHARTHICK = 2, CHARSIZE=1, COLOR=2 169 xyouts, julday(5,15,1900), 0.43, 'Tot. Annual Export Model 2= '+strtrim(tot_area_expo_2, 1)+' million Km2', ALIGN = 0, CHARTHICK = 2, CHARSIZE=1, COLOR=2 170 xyouts, julday(7,15,1900), 0.41, 'Data from Kwok et al.(2004), 1979-2002 ', ALIGN = 0, CHARTHICK = 2, CHARSIZE=0.8, COLOR=2 171 ; 171 172 title = 'Fram Strait Volume Export'+'!C'+prefix+' (BLACK) - '+prefix2+' (RED) '+d1_d2_2 172 pltt, ICE_vol_export, 't', /REMPLI, MIN = 0., MAX = .6, /NOERASE, XGRIDSTYLE = 1 $173 pltt, ICE_vol_export, 't', /REMPLI, MIN = 60., MAX = 500., /NOERASE, XGRIDSTYLE = 1 $ 173 174 , small = [1, 2, 2], YTITLE = '10^3 Km^3/month', TITLE = title, DATE_FORMAT = '%M', _extra = ex 174 175 pltt, ICE_vol_export_2, 't', /REMPLI, /NOERASE $ … … 181 182 tot_vol_expo_Sv = tot_vol_expo * 1.e06 * 1/86400 * 1/365 ; annual mean in Sverdrup 182 183 tot_vol_expo_2_Sv = tot_vol_expo_2 * 1.e06 * 1/86400 * 1/365 ; annual mean in Sverdrup 183 xyouts, julday( 6,15,1900), 0.48, 'Tot. Annual Export OBS = 2124Km3/year', ALIGN = 0, CHARTHICK = 2, CHARSIZE=1, COLOR=2184 xyouts, julday( 6,15,1900), 0.46, 'Tot. Annual Export Model 1 = '+strtrim(tot_vol_expo, 1)+'Km3/year', ALIGN = 0, CHARTHICK = 2, CHARSIZE=1, COLOR=2185 xyouts, julday(9,15,1900), 0.44, 'in Sv = '+strtrim(tot_vol_expo_Sv, 1)+' Sv', ALIGN = 0, CHARTHICK = 2, CHARSIZE=1, COLOR=2186 xyouts, julday( 6,15,1900), 0.42, 'Tot. Annual Export Model 2 = '+strtrim(tot_vol_expo_2, 1)+'Km3/year', ALIGN = 0, CHARTHICK = 2, CHARSIZE=1, COLOR=2187 xyouts, julday(9,15,1900), 0.40, 'in Sv = '+strtrim(tot_vol_expo_2_Sv, 1)+' Sv', ALIGN = 0, CHARTHICK = 2, CHARSIZE=1, COLOR=2188 xyouts, julday(7,15,1900), 0.38, 'Data from Kwok et al.(2004), 1992-1998 ', ALIGN = 0, CHARTHICK = 2, CHARSIZE=0.8, COLOR=2184 xyouts, julday(5,15,1900), 490, 'Tot. Annual Export OBS = 2124 10^3 Km3/year', ALIGN = 0, CHARTHICK = 2, CHARSIZE=1, COLOR=2 185 xyouts, julday(5,15,1900), 460, 'Tot. Annual Export Model 1 = '+strtrim(tot_vol_expo, 1)+' 10^3 Km3/year', ALIGN = 0, CHARTHICK = 2, CHARSIZE=1, COLOR=2 186 xyouts, julday(9,15,1900), 440, 'in Sv = '+strtrim(tot_vol_expo_Sv, 1)+' Sv', ALIGN = 0, CHARTHICK = 2, CHARSIZE=1, COLOR=2 187 xyouts, julday(5,15,1900), 410, 'Tot. Annual Export Model 2 = '+strtrim(tot_vol_expo_2, 1)+' 10^3 Km3/year', ALIGN = 0, CHARTHICK = 2, CHARSIZE=1, COLOR=2 188 xyouts, julday(9,15,1900), 390, 'in Sv = '+strtrim(tot_vol_expo_2_Sv, 1)+' Sv', ALIGN = 0, CHARTHICK = 2, CHARSIZE=1, COLOR=2 189 xyouts, julday(7,15,1900), 360, 'Data from Kwok et al.(2004), 1992-1998 ', ALIGN = 0, CHARTHICK = 2, CHARSIZE=0.8, COLOR=2 189 190 190 191 htmltxt = [ htmltxt, '<hr>'+blabla, '<br><img width="80%" src='+filename+'_2.png /> ' ] -
branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_all.pro
r4504 r4514 78 78 IF doplot EQ cnt OR doplot EQ 0 THEN std_ts_ICE_Vol, masknp, POSTSCRIPT = postscript, _extra = ex 79 79 ; 80 ;SF cnt = cnt+1 & blabla = 'ICE Velocity '81 ;SF IF doplot EQ cnt OR doplot EQ 0 THEN std_ts_ICE_Vel, masknp, POSTSCRIPT = postscript, _extra = ex80 ;SF cnt = cnt+1 & blabla = 'ICE Velocity ' 81 ;SF IF doplot EQ cnt OR doplot EQ 0 THEN std_ts_ICE_Vel, masknp, POSTSCRIPT = postscript, _extra = ex 82 82 ; 83 83 cnt = cnt+1 & blabla = 'Sea-Ice Volume Export at Fram Strait'
Note: See TracChangeset
for help on using the changeset viewer.