Changeset 4518
- Timestamp:
- 2014-02-25T16:20:43+01:00 (10 years ago)
- Location:
- branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts
- Files:
-
- 7 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_IceAge.pro
r4516 r4518 1 pro std_plot_IceAge, Iage1, Iage2, ARC = arc, ANT = ant, MARCH = march, SEPT = sept, POSTSCRIPT = postscript, _extra = ex1 pro std_plot_IceAge, Iage1, Iage2, 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(febr) THEN var = var+'Febr' 14 15 IF keyword_set(march) THEN var = var+'March' 15 16 IF keyword_set(sept) THEN var = var+'Sept' -
branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_IceFrac.pro
r4320 r4518 1 pro std_plot_IceFrac, Ifra1, Ifra2, ARC = arc, ANT = ant, MARCH = march, SEPT = sept, POSTSCRIPT = postscript, _extra = ex1 pro std_plot_IceFrac, Ifra1, Ifra2, 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(febr) THEN var = var+'Febr' 14 15 IF keyword_set(march) THEN var = var+'March' 15 16 IF keyword_set(sept) THEN var = var+'Sept' -
branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_IceThick.pro
r4320 r4518 1 pro std_plot_IceThick, Ithi1, Ithi2, ARC = arc, ANT = ant, APRIL = april, SEPT = sept, POSTSCRIPT = postscript, _extra = ex1 pro std_plot_IceThick, Ithi1, Ithi2, ARC = arc, ANT = ant, APRIL = april, JAN = jan, SEPT = sept, POSTSCRIPT = postscript, _extra = ex 2 2 3 3 compile_opt idl2, strictarrsubs … … 13 13 IF keyword_set(ant) THEN var = var+'_Ant_' 14 14 IF keyword_set(april) THEN var = var+'April' 15 IF keyword_set(jan) THEN var = var+'Jan' 15 16 IF keyword_set(sept) THEN var = var+'Sept' 16 17 -
branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_SnowThick.pro
r4514 r4518 1 pro std_plot_SnowThick, Isnow1, Isnow2, ARC = arc, ANT = ant, MARCH = march, SEPT = sept, POSTSCRIPT = postscript, _extra = ex1 pro std_plot_SnowThick, Isnow1, Isnow2, 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(febr) THEN var = var+'Febr' 14 15 IF keyword_set(march) THEN var = var+'March' 15 16 IF keyword_set(sept) THEN var = var+'Sept' -
branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_all.pro
r4517 r4518 103 103 ;climatology 104 104 QNET = read_ncdf(getenv('VAR_FLUX'), filename = std_file_oaflux ) 105 105 106 106 ; erp (evaporation damping) 107 107 ERP1 = read_ncdf(getenv('VAR1_ERP'), allrecords = allrec, direc = 't', filename = std_file1_T ) … … 170 170 april = where(mm EQ 4, cnt) 171 171 Ithi_april_1 = {arr:1./float(cnt) * total(reform(Ithi_1.arr[*, *, temporary(april)],nxt,nyt,cnt), 3), unit:Ithi_1.unit} 172 jan = where(mm EQ 1, cnt) 172 173 sept = where(mm EQ 9, cnt) 174 Ithi_jan_1 = {arr:1./float(cnt) * total(reform(Ithi_1.arr[*, *, temporary(jan)],nxt,nyt,cnt), 3), unit:Ithi_1.unit} 173 175 Ithi_sept_1 = {arr:1./float(cnt) * total(reform(Ithi_1.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Ithi_1.unit} 174 176 undefine, Ithi_1 … … 179 181 april = where(mm EQ 4, cnt) 180 182 Ithi_april_2 = {arr:1./float(cnt) * total(reform(Ithi_2.arr[*, *, temporary(april)],nxt,nyt,cnt), 3), unit:Ithi_2.unit} 183 jan = where(mm EQ 1, cnt) 181 184 sept = where(mm EQ 9, cnt) 182 185 Ithi_sept_2 = {arr:1./float(cnt) * total(reform(Ithi_2.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Ithi_2.unit} 186 Ithi_jan_2 = {arr:1./float(cnt) * total(reform(Ithi_2.arr[*, *, temporary(jan)],nxt,nyt,cnt), 3), unit:Ithi_2.unit} 183 187 undefine, Ithi_2 184 188 ENDIF ELSE BEGIN 185 189 Ithi_april_2 = {arr:-1} 186 190 Ithi_sept_2 = {arr:-1} 191 Ithi_jan_2 = {arr:-1} 187 192 ENDELSE 188 193 ; 189 194 Iage_1 = read_ncdf(getenv('VAR1_Iage'), allrecords = allrec, filename = std_file1_I ) 190 195 caldat, time, mm 196 febr = where(mm EQ 2, cnt) 197 Iage_febr_1 = {arr:1./float(cnt) * total(reform(Iage_1.arr[*, *, temporary(febr)],nxt,nyt,cnt), 3), unit:Iage_1.unit} 191 198 march = where(mm EQ 3, cnt) 192 199 Iage_march_1 = {arr:1./float(cnt) * total(reform(Iage_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Iage_1.unit} … … 198 205 Iage_2 = read_ncdf(getenv('VAR2_Iage'), allrecords = allrec, filename = std_file2_I ) 199 206 caldat, time, mm 207 febr = where(mm EQ 2, cnt) 208 Iage_febr_2 = {arr:1./float(cnt) * total(reform(Iage_2.arr[*, *, temporary(febr)],nxt,nyt,cnt), 3), unit:Iage_2.unit} 200 209 march = where(mm EQ 3, cnt) 201 210 Iage_march_2 = {arr:1./float(cnt) * total(reform(Iage_2.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Iage_2.unit} … … 204 213 undefine, Iage_2 205 214 ENDIF ELSE BEGIN 215 Iage_febr_2 = {arr:-1} 206 216 Iage_march_2 = {arr:-1} 207 217 Iage_sept_2 = {arr:-1} … … 210 220 Ifra_1 = read_ncdf(getenv('VAR1_Ifrac'), allrecords = allrec, filename = std_file1_I ) 211 221 caldat, time, mm 222 febr = where(mm EQ 2, cnt) 223 Ifra_febr_1 = {arr:1./float(cnt) * total(reform(Ifra_1.arr[*, *, temporary(febr)],nxt,nyt,cnt), 3), unit:Ifra_1.unit} 212 224 march = where(mm EQ 3, cnt) 213 225 Ifra_march_1 = {arr:1./float(cnt) * total(reform(Ifra_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Ifra_1.unit} … … 219 231 Ifra_2 = read_ncdf(getenv('VAR2_Ifrac'), allrecords = allrec, filename = std_file2_I ) 220 232 caldat, time, mm 233 febr = where(mm EQ 2, cnt) 234 Ifra_febr_2 = {arr:1./float(cnt) * total(reform(Ifra_2.arr[*, *, temporary(febr)],nxt,nyt,cnt), 3), unit:Ifra_2.unit} 221 235 march = where(mm EQ 3, cnt) 222 236 Ifra_march_2 = {arr:1./float(cnt) * total(reform(Ifra_2.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Ifra_2.unit} … … 225 239 undefine, Ifra_2 226 240 ENDIF ELSE BEGIN 241 Ifra_febr_2 = {arr:-1} 227 242 Ifra_march_2 = {arr:-1} 228 243 Ifra_sept_2 = {arr:-1} … … 231 246 Isnow_1 = read_ncdf(getenv('VAR1_Isnow'), allrecords = allrec, filename = std_file1_I ) 232 247 caldat, time, mm 248 febr = where(mm EQ 2, cnt) 249 Isnow_febr_1 = {arr:1./float(cnt) * total(reform(Isnow_1.arr[*, *, temporary(febr)],nxt,nyt,cnt), 3), unit:Isnow_1.unit} 233 250 march = where(mm EQ 3, cnt) 234 251 Isnow_march_1 = {arr:1./float(cnt) * total(reform(Isnow_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Isnow_1.unit} … … 240 257 Isnow_2 = read_ncdf(getenv('VAR2_Isnow'), allrecords = allrec, filename = std_file2_I ) 241 258 caldat, time, mm 259 febr = where(mm EQ 2, cnt) 260 Isnow_febr_2 = {arr:1./float(cnt) * total(reform(Isnow_2.arr[*, *, temporary(febr)],nxt,nyt,cnt), 3), unit:Isnow_2.unit} 242 261 march = where(mm EQ 3, cnt) 243 262 Isnow_march_2 = {arr:1./float(cnt) * total(reform(Isnow_2.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Isnow_2.unit} … … 246 265 ; undefine, Isnow_2 247 266 ENDIF ELSE BEGIN 267 Isnow_febr_2 = {arr:-1} 248 268 Isnow_march_2 = {arr:-1} 249 269 Isnow_sept_2 = {arr:-1} … … 428 448 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_Med_Sdepth, S1, S2, SLev, 38, _extra = ex 429 449 ; 430 ; cnt = cnt+1 & blabla = 'Arctic Ice Thickness: MARCH'431 ; IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_march_1, Ithi_march_2, /ARC, /MARCH, _extra = ex450 cnt = cnt+1 & blabla = 'Arctic Ice Thickness: JAN' 451 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_jan_1, Ithi_jan_2, /ARC, /JAN, _extra = ex 432 452 ; 433 453 cnt = cnt+1 & blabla = 'Arctic Ice Thickness: APRIL' 434 454 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_april_1, Ithi_april_2, /ARC, /APRIL, _extra = ex 435 ; 436 cnt = cnt+1 & blabla = 'Arctic Ice Thickness: SEPT' 437 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_sept_1, Ithi_sept_2, /ARC, /SEPT, _extra = ex 438 ; 439 ; cnt = cnt+1 & blabla = 'Antarctic Ice Thickness: MARCH' 440 ; IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_march_1, Ithi_march_2, /ANT, /MARCH, _extra = ex 441 ; 455 ; 442 456 cnt = cnt+1 & blabla = 'Antarctic Ice Thickness: APRIL' 443 457 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_april_1, Ithi_april_2, /ANT, /APRIL, _extra = ex … … 453 467 ; 454 468 cnt = cnt+1 & blabla = 'Antarctic Ice Age: MARCH' 455 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceAge, Iage_ march_1, Iage_march_2, /ANT, /MARCH, _extra = ex469 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceAge, Iage_febr_1, Iage_febr_2, /ANT, /FEBR, _extra = ex 456 470 ; 457 471 cnt = cnt+1 & blabla = 'Antarctic Ice Age: SEPT' … … 465 479 ; 466 480 cnt = cnt+1 & blabla = 'Antarctic Ice Fraction: MARCH' 467 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceFrac, Ifra_ march_1, Ifra_march_2, /ANT, /MARCH, _extra = ex481 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceFrac, Ifra_febr_1, Ifra_febr_2, /ANT, /FEBR, _extra = ex 468 482 ; 469 483 cnt = cnt+1 & blabla = 'Antarctic Ice Fraction: SEPT' … … 477 491 ; 478 492 cnt = cnt+1 & blabla = 'Antarctic SNOW Thickness: MARCH' 479 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_ march_1, Isnow_march_2, /ANT, /MARCH, _extra = ex493 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_febr_1, Isnow_febr_2, /ANT, /FEBR, _extra = ex 480 494 ; 481 495 cnt = cnt+1 & blabla = 'Antarctic SNOW Thickness: SEPT' -
branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_all_1Y.pro
r4517 r4518 1 pro std_plot_all , doplot = doplot, _extra = ex1 pro std_plot_all_1Y, doplot = doplot, _extra = ex 2 2 3 3 compile_opt idl2, strictarrsubs … … 103 103 ;climatology 104 104 QNET = read_ncdf(getenv('VAR_FLUX'), filename = std_file_oaflux ) 105 105 106 106 ; erp (evaporation damping) 107 107 ERP1 = read_ncdf(getenv('VAR1_ERP'), allrecords = allrec, direc = 't', filename = std_file1_T ) … … 170 170 april = where(mm EQ 4, cnt) 171 171 Ithi_april_1 = {arr:1./float(cnt) * total(reform(Ithi_1.arr[*, *, temporary(april)],nxt,nyt,cnt), 3), unit:Ithi_1.unit} 172 jan = where(mm EQ 1, cnt) 172 173 sept = where(mm EQ 9, cnt) 174 Ithi_jan_1 = {arr:1./float(cnt) * total(reform(Ithi_1.arr[*, *, temporary(jan)],nxt,nyt,cnt), 3), unit:Ithi_1.unit} 173 175 Ithi_sept_1 = {arr:1./float(cnt) * total(reform(Ithi_1.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Ithi_1.unit} 174 176 undefine, Ithi_1 … … 179 181 april = where(mm EQ 4, cnt) 180 182 Ithi_april_2 = {arr:1./float(cnt) * total(reform(Ithi_2.arr[*, *, temporary(april)],nxt,nyt,cnt), 3), unit:Ithi_2.unit} 183 jan = where(mm EQ 1, cnt) 181 184 sept = where(mm EQ 9, cnt) 182 185 Ithi_sept_2 = {arr:1./float(cnt) * total(reform(Ithi_2.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Ithi_2.unit} 186 Ithi_jan_2 = {arr:1./float(cnt) * total(reform(Ithi_2.arr[*, *, temporary(jan)],nxt,nyt,cnt), 3), unit:Ithi_2.unit} 183 187 undefine, Ithi_2 184 188 ENDIF ELSE BEGIN 185 189 Ithi_april_2 = {arr:-1} 186 190 Ithi_sept_2 = {arr:-1} 191 Ithi_jan_2 = {arr:-1} 187 192 ENDELSE 188 193 ; 189 194 Iage_1 = read_ncdf(getenv('VAR1_Iage'), allrecords = allrec, filename = std_file1_I ) 190 195 caldat, time, mm 196 febr = where(mm EQ 2, cnt) 197 Iage_febr_1 = {arr:1./float(cnt) * total(reform(Iage_1.arr[*, *, temporary(febr)],nxt,nyt,cnt), 3), unit:Iage_1.unit} 191 198 march = where(mm EQ 3, cnt) 192 199 Iage_march_1 = {arr:1./float(cnt) * total(reform(Iage_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Iage_1.unit} … … 198 205 Iage_2 = read_ncdf(getenv('VAR2_Iage'), allrecords = allrec, filename = std_file2_I ) 199 206 caldat, time, mm 207 febr = where(mm EQ 2, cnt) 208 Iage_febr_2 = {arr:1./float(cnt) * total(reform(Iage_2.arr[*, *, temporary(febr)],nxt,nyt,cnt), 3), unit:Iage_2.unit} 200 209 march = where(mm EQ 3, cnt) 201 210 Iage_march_2 = {arr:1./float(cnt) * total(reform(Iage_2.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Iage_2.unit} … … 204 213 undefine, Iage_2 205 214 ENDIF ELSE BEGIN 215 Iage_febr_2 = {arr:-1} 206 216 Iage_march_2 = {arr:-1} 207 217 Iage_sept_2 = {arr:-1} … … 210 220 Ifra_1 = read_ncdf(getenv('VAR1_Ifrac'), allrecords = allrec, filename = std_file1_I ) 211 221 caldat, time, mm 222 febr = where(mm EQ 2, cnt) 223 Ifra_febr_1 = {arr:1./float(cnt) * total(reform(Ifra_1.arr[*, *, temporary(febr)],nxt,nyt,cnt), 3), unit:Ifra_1.unit} 212 224 march = where(mm EQ 3, cnt) 213 225 Ifra_march_1 = {arr:1./float(cnt) * total(reform(Ifra_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Ifra_1.unit} … … 219 231 Ifra_2 = read_ncdf(getenv('VAR2_Ifrac'), allrecords = allrec, filename = std_file2_I ) 220 232 caldat, time, mm 233 febr = where(mm EQ 2, cnt) 234 Ifra_febr_2 = {arr:1./float(cnt) * total(reform(Ifra_2.arr[*, *, temporary(febr)],nxt,nyt,cnt), 3), unit:Ifra_2.unit} 221 235 march = where(mm EQ 3, cnt) 222 236 Ifra_march_2 = {arr:1./float(cnt) * total(reform(Ifra_2.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Ifra_2.unit} … … 225 239 undefine, Ifra_2 226 240 ENDIF ELSE BEGIN 241 Ifra_febr_2 = {arr:-1} 227 242 Ifra_march_2 = {arr:-1} 228 243 Ifra_sept_2 = {arr:-1} … … 231 246 Isnow_1 = read_ncdf(getenv('VAR1_Isnow'), allrecords = allrec, filename = std_file1_I ) 232 247 caldat, time, mm 248 febr = where(mm EQ 2, cnt) 249 Isnow_febr_1 = {arr:1./float(cnt) * total(reform(Isnow_1.arr[*, *, temporary(febr)],nxt,nyt,cnt), 3), unit:Isnow_1.unit} 233 250 march = where(mm EQ 3, cnt) 234 251 Isnow_march_1 = {arr:1./float(cnt) * total(reform(Isnow_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Isnow_1.unit} … … 240 257 Isnow_2 = read_ncdf(getenv('VAR2_Isnow'), allrecords = allrec, filename = std_file2_I ) 241 258 caldat, time, mm 259 febr = where(mm EQ 2, cnt) 260 Isnow_febr_2 = {arr:1./float(cnt) * total(reform(Isnow_2.arr[*, *, temporary(febr)],nxt,nyt,cnt), 3), unit:Isnow_2.unit} 242 261 march = where(mm EQ 3, cnt) 243 262 Isnow_march_2 = {arr:1./float(cnt) * total(reform(Isnow_2.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Isnow_2.unit} … … 246 265 ; undefine, Isnow_2 247 266 ENDIF ELSE BEGIN 267 Isnow_febr_2 = {arr:-1} 248 268 Isnow_march_2 = {arr:-1} 249 269 Isnow_sept_2 = {arr:-1} … … 428 448 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_Med_Sdepth, S1, S2, SLev, 38, _extra = ex 429 449 ; 430 ; cnt = cnt+1 & blabla = 'Arctic Ice Thickness: MARCH'431 ; IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_march_1, Ithi_march_2, /ARC, /MARCH, _extra = ex450 cnt = cnt+1 & blabla = 'Arctic Ice Thickness: JAN' 451 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_jan_1, Ithi_jan_2, /ARC, /JAN, _extra = ex 432 452 ; 433 453 cnt = cnt+1 & blabla = 'Arctic Ice Thickness: APRIL' 434 454 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_april_1, Ithi_april_2, /ARC, /APRIL, _extra = ex 435 ; 436 cnt = cnt+1 & blabla = 'Arctic Ice Thickness: SEPT' 437 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_sept_1, Ithi_sept_2, /ARC, /SEPT, _extra = ex 438 ; 439 ; cnt = cnt+1 & blabla = 'Antarctic Ice Thickness: MARCH' 440 ; IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_march_1, Ithi_march_2, /ANT, /MARCH, _extra = ex 441 ; 455 ; 442 456 cnt = cnt+1 & blabla = 'Antarctic Ice Thickness: APRIL' 443 457 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_april_1, Ithi_april_2, /ANT, /APRIL, _extra = ex … … 453 467 ; 454 468 cnt = cnt+1 & blabla = 'Antarctic Ice Age: MARCH' 455 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceAge, Iage_ march_1, Iage_march_2, /ANT, /MARCH, _extra = ex469 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceAge, Iage_febr_1, Iage_febr_2, /ANT, /FEBR, _extra = ex 456 470 ; 457 471 cnt = cnt+1 & blabla = 'Antarctic Ice Age: SEPT' … … 465 479 ; 466 480 cnt = cnt+1 & blabla = 'Antarctic Ice Fraction: MARCH' 467 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceFrac, Ifra_ march_1, Ifra_march_2, /ANT, /MARCH, _extra = ex481 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceFrac, Ifra_febr_1, Ifra_febr_2, /ANT, /FEBR, _extra = ex 468 482 ; 469 483 cnt = cnt+1 & blabla = 'Antarctic Ice Fraction: SEPT' … … 477 491 ; 478 492 cnt = cnt+1 & blabla = 'Antarctic SNOW Thickness: MARCH' 479 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_ march_1, Isnow_march_2, /ANT, /MARCH, _extra = ex493 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_febr_1, Isnow_febr_2, /ANT, /FEBR, _extra = ex 480 494 ; 481 495 cnt = cnt+1 & blabla = 'Antarctic SNOW Thickness: SEPT' -
branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_all_last_year_5D.pro
r4517 r4518 105 105 ;climatology 106 106 QNET = read_ncdf(getenv('VAR_FLUX'), filename = std_file_oaflux ) 107 107 108 108 ; erp (evaporation damping) 109 109 ERP1 = read_ncdf(getenv('VAR1_ERP'), allrecords = allrec, direc = 't', filename = std_file1_T ) … … 186 186 april = where(mm EQ 4, cnt) 187 187 Ithi_april_1 = {arr:1./float(cnt) * total(reform(Ithi_1.arr[*, *, temporary(april)],nxt,nyt,cnt), 3), unit:Ithi_1.unit} 188 jan = where(mm EQ 1, cnt) 188 189 sept = where(mm EQ 9, cnt) 190 Ithi_jan_1 = {arr:1./float(cnt) * total(reform(Ithi_1.arr[*, *, temporary(jan)],nxt,nyt,cnt), 3), unit:Ithi_1.unit} 189 191 Ithi_sept_1 = {arr:1./float(cnt) * total(reform(Ithi_1.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Ithi_1.unit} 190 192 undefine, Ithi_1 … … 195 197 april = where(mm EQ 4, cnt) 196 198 Ithi_april_2 = {arr:1./float(cnt) * total(reform(Ithi_2.arr[*, *, temporary(april)],nxt,nyt,cnt), 3), unit:Ithi_2.unit} 199 jan = where(mm EQ 1, cnt) 197 200 sept = where(mm EQ 9, cnt) 198 201 Ithi_sept_2 = {arr:1./float(cnt) * total(reform(Ithi_2.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Ithi_2.unit} 202 Ithi_jan_2 = {arr:1./float(cnt) * total(reform(Ithi_2.arr[*, *, temporary(jan)],nxt,nyt,cnt), 3), unit:Ithi_2.unit} 199 203 undefine, Ithi_2 200 204 ENDIF ELSE BEGIN 201 205 Ithi_april_2 = {arr:-1} 202 206 Ithi_sept_2 = {arr:-1} 207 Ithi_jan_2 = {arr:-1} 203 208 ENDELSE 204 209 ; 205 210 Iage_1 = read_ncdf(getenv('VAR1_Iage'), allrecords = allrec, filename = std_file1_I ) 206 211 caldat, time, mm 212 febr = where(mm EQ 2, cnt) 213 Iage_febr_1 = {arr:1./float(cnt) * total(reform(Iage_1.arr[*, *, temporary(febr)],nxt,nyt,cnt), 3), unit:Iage_1.unit} 207 214 march = where(mm EQ 3, cnt) 208 215 Iage_march_1 = {arr:1./float(cnt) * total(reform(Iage_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Iage_1.unit} … … 214 221 Iage_2 = read_ncdf(getenv('VAR2_Iage'), allrecords = allrec, filename = std_file2_I ) 215 222 caldat, time, mm 223 febr = where(mm EQ 2, cnt) 224 Iage_febr_2 = {arr:1./float(cnt) * total(reform(Iage_2.arr[*, *, temporary(febr)],nxt,nyt,cnt), 3), unit:Iage_2.unit} 216 225 march = where(mm EQ 3, cnt) 217 226 Iage_march_2 = {arr:1./float(cnt) * total(reform(Iage_2.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Iage_2.unit} … … 220 229 undefine, Iage_2 221 230 ENDIF ELSE BEGIN 231 Iage_febr_2 = {arr:-1} 222 232 Iage_march_2 = {arr:-1} 223 233 Iage_sept_2 = {arr:-1} … … 225 235 ; 226 236 Ifra_1 = read_ncdf(getenv('VAR1_Ifrac'), allrecords = allrec, filename = std_file1_I ) 227 ;; help, Ifra_1.arr, jpi,jpj228 237 caldat, time, mm 238 febr = where(mm EQ 2, cnt) 239 Ifra_febr_1 = {arr:1./float(cnt) * total(reform(Ifra_1.arr[*, *, temporary(febr)],nxt,nyt,cnt), 3), unit:Ifra_1.unit} 229 240 march = where(mm EQ 3, cnt) 230 241 Ifra_march_1 = {arr:1./float(cnt) * total(reform(Ifra_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Ifra_1.unit} … … 236 247 Ifra_2 = read_ncdf(getenv('VAR2_Ifrac'), allrecords = allrec, filename = std_file2_I ) 237 248 caldat, time, mm 249 febr = where(mm EQ 2, cnt) 250 Ifra_febr_2 = {arr:1./float(cnt) * total(reform(Ifra_2.arr[*, *, temporary(febr)],nxt,nyt,cnt), 3), unit:Ifra_2.unit} 238 251 march = where(mm EQ 3, cnt) 239 252 Ifra_march_2 = {arr:1./float(cnt) * total(reform(Ifra_2.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Ifra_2.unit} … … 242 255 undefine, Ifra_2 243 256 ENDIF ELSE BEGIN 257 Ifra_febr_2 = {arr:-1} 244 258 Ifra_march_2 = {arr:-1} 245 259 Ifra_sept_2 = {arr:-1} … … 248 262 Isnow_1 = read_ncdf(getenv('VAR1_Isnow'), allrecords = allrec, filename = std_file1_I ) 249 263 caldat, time, mm 264 febr = where(mm EQ 2, cnt) 265 Isnow_febr_1 = {arr:1./float(cnt) * total(reform(Isnow_1.arr[*, *, temporary(febr)],nxt,nyt,cnt), 3), unit:Isnow_1.unit} 250 266 march = where(mm EQ 3, cnt) 251 267 Isnow_march_1 = {arr:1./float(cnt) * total(reform(Isnow_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Isnow_1.unit} … … 257 273 Isnow_2 = read_ncdf(getenv('VAR2_Isnow'), allrecords = allrec, filename = std_file2_I ) 258 274 caldat, time, mm 275 febr = where(mm EQ 2, cnt) 276 Isnow_febr_2 = {arr:1./float(cnt) * total(reform(Isnow_2.arr[*, *, temporary(febr)],nxt,nyt,cnt), 3), unit:Isnow_2.unit} 259 277 march = where(mm EQ 3, cnt) 260 278 Isnow_march_2 = {arr:1./float(cnt) * total(reform(Isnow_2.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Isnow_2.unit} … … 263 281 ; undefine, Isnow_2 264 282 ENDIF ELSE BEGIN 283 Isnow_febr_2 = {arr:-1} 265 284 Isnow_march_2 = {arr:-1} 266 285 Isnow_sept_2 = {arr:-1} … … 450 469 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_Med_Sdepth, S1, S2, SLev, 38, _extra = ex 451 470 ; 452 ; cnt = cnt+1 & blabla = 'Arctic Ice Thickness: MARCH'453 ; IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_march_1, Ithi_march_2, /ARC, /MARCH, _extra = ex471 cnt = cnt+1 & blabla = 'Arctic Ice Thickness: JAN' 472 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_jan_1, Ithi_jan_2, /ARC, /JAN, _extra = ex 454 473 ; 455 474 cnt = cnt+1 & blabla = 'Arctic Ice Thickness: APRIL' 456 475 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_april_1, Ithi_april_2, /ARC, /APRIL, _extra = ex 457 ; 458 cnt = cnt+1 & blabla = 'Arctic Ice Thickness: SEPT' 459 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_sept_1, Ithi_sept_2, /ARC, /SEPT, _extra = ex 460 ; 461 ; cnt = cnt+1 & blabla = 'Antarctic Ice Thickness: MARCH' 462 ; IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_march_1, Ithi_march_2, /ANT, /MARCH, _extra = ex 463 ; 476 ; 464 477 cnt = cnt+1 & blabla = 'Antarctic Ice Thickness: APRIL' 465 478 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_april_1, Ithi_april_2, /ANT, /APRIL, _extra = ex … … 475 488 ; 476 489 cnt = cnt+1 & blabla = 'Antarctic Ice Age: MARCH' 477 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceAge, Iage_ march_1, Iage_march_2, /ANT, /MARCH, _extra = ex490 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceAge, Iage_febr_1, Iage_febr_2, /ANT, /FEBR, _extra = ex 478 491 ; 479 492 cnt = cnt+1 & blabla = 'Antarctic Ice Age: SEPT' … … 487 500 ; 488 501 cnt = cnt+1 & blabla = 'Antarctic Ice Fraction: MARCH' 489 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceFrac, Ifra_ march_1, Ifra_march_2, /ANT, /MARCH, _extra = ex502 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceFrac, Ifra_febr_1, Ifra_febr_2, /ANT, /FEBR, _extra = ex 490 503 ; 491 504 cnt = cnt+1 & blabla = 'Antarctic Ice Fraction: SEPT' … … 499 512 ; 500 513 cnt = cnt+1 & blabla = 'Antarctic SNOW Thickness: MARCH' 501 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_ march_1, Isnow_march_2, /ANT, /MARCH, _extra = ex514 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_febr_1, Isnow_febr_2, /ANT, /FEBR, _extra = ex 502 515 ; 503 516 cnt = cnt+1 & blabla = 'Antarctic SNOW Thickness: SEPT'
Note: See TracChangeset
for help on using the changeset viewer.