Changeset 51 for tags


Ignore:
Timestamp:
12/12/17 12:09:14 (7 years ago)
Author:
vancop
Message:

N_ICE_RL config modifs

Location:
tags/LIM1D_v3.20
Files:
7 edited

Legend:

Unmodified
Added
Removed

  • tags/LIM1D_v3.20/GRAPHICS/ALL_contours.pro

    r6 r51  
    129129ENDIF 
    130130 
     131IF ( site EQ 'N_ICE_RL' ) THEN BEGIN 
     132   mon = numt * dt / 86400. / 30.41666666 + 1 
     133ENDIF 
     134 
    131135; 
    132136;------------------------------------------------------------------------------ 
    … …  
    199203ENDIF 
    200204 
     205IF ( site EQ 'N_ICE_RL' ) THEN BEGIN ; all wrong here, just a silly copy paste 
     206   day_obs = [ 120,  120,   120,   121,   121,   121,   124,   124,   124,   124,   125,   125,   125,   125,   125,   125,   125,   125,   125,   126,   126,   126,   127,   127,   127,   127,   130,   132,   134,   136,   138,   140,   142,   144,   146,   149,   152,   154,   127,   130,   132,   134,   136,   138,   140,   142,   144,   146,   149,   152,   154,   127,   130,   132,   134,   136,   138,   140,   142,   144,   146,   149,   152,   154,   127,   130,   132,   134,   136,   138,   140,   142,   144,   146,   149,   152,   154,   127,   130,   132,   134,   136,   138,   140,   142,   144,   146,   149,   152,   154 ] 
     207 
     208   hi_obs  = [   16,   17,   13,   17,   17,   17,   16,   16,   16,   17,   16,   17.5,   20,   14,   17,   14,   18,   30,   19,   15,   15,   15,   19,   18,   19,   17,   19.35,   21.25,   22,   22.1,   22.2,   22,   21.5,   22,   22,   23,   21,   19,   24,   22.25,   23.25,   25.25,   26,   26,   26.5,   24.5,   24,   24,   23,   23,   20.5,   21.75,   23.25,   24.75,   25,   26.5,   27,   26,   25,   25,   25,   24,   23,   21,   21.5,   22.75,   24.5,   25.4,   25.75,   25.7,   25.5,   24,   24.5,   25,   24.5,   23,   21.5,   20.25,   22.25,   23.88,   24.75,   25,   24.7,   25,   24,   24.5,   24.5,   24,   23,   20 ]   / 100. 
     209   hi_obs_std = REPLICATE( 0., N_ELEMENTS(day_obs) ) 
     210   hs_obs     = REPLICATE( 0., N_ELEMENTS(day_obs) ) 
     211   hs_obs_std = REPLICATE( 0., N_ELEMENTS(day_obs) ) 
     212   mon_obs    = day_obs / 30.467 + 1 
     213ENDIF 
     214 
    201215; 
    202216;------------------------------------------------------------------------------ 
    … …  
    231245ENDIF 
    232246 
     247IF ( site EQ 'N_ICE_RL' ) THEN BEGIN 
     248   zmin = -0.5 
     249   zmax = 0.0 
     250ENDIF 
     251 
    233252;------------------- 
    234253; Salinity contours 
    … …  
    416435; contour parameters 
    417436nlevels = 10 
     437 
    418438colors = [   20,  40,   60,   80,  100,  120, 140, 160, 180, 200, 220 ] 
    419439IF ( site EQ 'BARROW' ) THEN BEGIN 
    420440   levels = [   0.,   0.2,  0.4,   0.6,   0.8,  1.0,  1.2, 1.4, 1.6, 1.8, 2.0] 
    421441ENDIF 
     442 
    422443IF ( site EQ 'ISPOL' ) THEN BEGIN 
    423444   levels = [   0.,     2.,   4.,    6.,    8.,  10.,   12.,14., 16., 18., 20.] 
    … …  
    426447   colors = [ 20, 60, 100, 140, 180, 220, 240 ] 
    427448ENDIF 
    428 IF ( site EQ 'YROSIAE' ) THEN BEGIN 
     449 
     450IF ( site EQ 'YROSIAE' OR site EQ 'N_ICE_RL' ) THEN BEGIN 
    429451   nlevels = 7. 
    430452   levels = [ 0., 1.6, 4., 6., 8., 10., 12., 20.] 
    … …  
    463485   colors = [ 20, 60, 100, 140, 180, 220, 240 ] 
    464486ENDIF 
    465 IF ( site EQ 'YROSIAE' ) THEN BEGIN 
     487IF ( site EQ 'YROSIAE' OR site EQ 'N_ICE_RL' ) THEN BEGIN 
    466488   nlevels = 7 
    467489   levels = [ 0., 0.24, 0.4, 0.6, 0.8, 1., 2. ] 
    … …  
    503525   colors = [ 20, 60, 100, 140, 180, 220, 240 ] 
    504526ENDIF 
    505 IF ( site EQ 'YROSIAE' ) THEN BEGIN 
     527IF ( site EQ 'YROSIAE' OR site EQ 'N_ICE_RL' ) THEN BEGIN 
    506528   nlevels = 7. 
    507529   levels = [ 0., 3.9, 8., 12., 16., 20., 24., 40.] 
    … …  
    12251247colorbarn, horpos, vertpos, cb_length, 0.02, MIN(levels), MAX(levels), leveling, nlevels+1, levels, colors, format 
    12261248 
    1227 ; 
    1228 ;-------------- 
    1229 ; dFe          
    1230 ;-------------- 
    1231 ; 
    1232 !P.MULTI=[0,numplot_x, numplot_y] 
    1233 i_plot = 0 
    1234 i_plot_h = i_plot MOD numplot_x 
    1235 i_plot_v = i_plot / numplot_x + 1 
    1236  
    1237 ; contour parameters 
    1238 nlevels = 10 
    1239 nlevels = 10 
    1240 fmin = 0.0 & fmax = MAX(dFeb) ; 
    1241 cmin = 000 & cmax = 255 
    1242 do_levels, fmin, fmax, cmin, cmax, nlevels, levels, colors 
    1243  
    1244 ct = 3 ; colortable 
    1245 title = ' dFe (mumol/m3)' 
    1246  
    1247 icecontour, dFeb, h_i, h_s, nlevels, colors, levels, ct, z_ib, mon, nlay_bio, nts, $ 
    1248             tmin, tmax, zmin, zmax, title 
    1249  
    1250 ; colorbar 
    1251 format='(f5.1)' 
    1252 horpos   = 0.02 + i_plot_h * dh 
    1253 vertpos  = 0.015 + ( numplot_y - i_plot_v ) * dv ; position of colorbar. 
    1254 colorbarn, horpos, vertpos, cb_length, 0.02, MIN(levels), MAX(levels), leveling, nlevels+1, levels, colors, format 
    1255  
    1256 ; 
    1257 ;-------------- 
    1258 ; aFe          
    1259 ;-------------- 
    1260 ; 
    1261 i_plot = i_plot + 1 
    1262 i_plot_h = i_plot MOD numplot_x 
    1263 i_plot_v = i_plot / numplot_x + 1 
    1264  
    1265 ; contour parameters 
    1266 nlevels = 10 
    1267 nlevels = 10 
    1268 fmin = 0.0 & fmax = MAX(aFeb) ; 
    1269 cmin = 000 & cmax = 255 
    1270 do_levels, fmin, fmax, cmin, cmax, nlevels, levels, colors 
    1271  
    1272 ct = 3 ; colortable 
    1273 title = ' aFe (mumol/m3)' 
    1274  
    1275 icecontour, aFeb, h_i, h_s, nlevels, colors, levels, ct, z_ib, mon, nlay_bio, nts, $ 
    1276             tmin, tmax, zmin, zmax, title 
    1277  
    1278 ; colorbar 
    1279 format='(f5.1)' 
    1280 horpos   = 0.02 + i_plot_h * dh 
    1281 vertpos  = 0.015 + ( numplot_y - i_plot_v ) * dv ; position of colorbar. 
    1282 colorbarn, horpos, vertpos, cb_length, 0.02, MIN(levels), MAX(levels), leveling, nlevels+1, levels, colors, format 
    1283  
    1284 ; 
    1285 ;-------------- 
    1286 ; eFe          
    1287 ;-------------- 
    1288 ; 
    1289 i_plot = i_plot + 1 
    1290 i_plot_h = i_plot MOD numplot_x 
    1291 i_plot_v = i_plot / numplot_x + 1 
    1292  
    1293 ; contour parameters 
    1294 nlevels = 10 
    1295 nlevels = 10 
    1296 fmin = 0.0 & fmax = MAX(eFeb) ; 
    1297 cmin = 000 & cmax = 255 
    1298 do_levels, fmin, fmax, cmin, cmax, nlevels, levels, colors 
    1299  
    1300 ct = 3 ; colortable 
    1301 title = ' eFe (mumol/m3)' 
    1302  
    1303 icecontour, eFeb, h_i, h_s, nlevels, colors, levels, ct, z_ib, mon, nlay_bio, nts, $ 
    1304             tmin, tmax, zmin, zmax, title 
    1305  
    1306 ; colorbar 
    1307 format='(f5.1)' 
    1308 horpos   = 0.02 + i_plot_h * dh 
    1309 vertpos  = 0.015 + ( numplot_y - i_plot_v ) * dv ; position of colorbar. 
    1310 colorbarn, horpos, vertpos, cb_length, 0.02, MIN(levels), MAX(levels), leveling, nlevels+1, levels, colors, format 
    1311  
    1312 ; 
    1313 ;------------------------------------------------------------------------------ 
    1314 ; End of the script 
    1315 ;------------------------------------------------------------------------------ 
    1316 ; 
    13171249IF ( device EQ 'PS' ) THEN BEGIN 
    13181250   DEVICE, /CLOSE 
    … …  
    13201252   !P.MULTI=[0,2,2] 
    13211253ENDIF 
     1254 
     1255STOP 
     1256 
     1257; 
     1258;-------------- 
     1259; dFe          
     1260;-------------- 
     1261; 
     1262!P.MULTI=[0,numplot_x, numplot_y] 
     1263i_plot = 0 
     1264i_plot_h = i_plot MOD numplot_x 
     1265i_plot_v = i_plot / numplot_x + 1 
     1266 
     1267; contour parameters 
     1268nlevels = 10 
     1269nlevels = 10 
     1270fmin = 0.0 & fmax = 60. 
     1271cmin = 000 & cmax = 255 
     1272do_levels, fmin, fmax, cmin, cmax, nlevels, levels, colors 
     1273 
     1274ct = 3 ; colortable 
     1275title = ' dFe (mumol/m3)' 
     1276 
     1277icecontour, dFeb, h_i, h_s, nlevels, colors, levels, ct, z_ib, mon, nlay_bio, nts, $ 
     1278            tmin, tmax, zmin, zmax, title 
     1279 
     1280; colorbar 
     1281format='(f5.1)' 
     1282horpos   = 0.02 + i_plot_h * dh 
     1283vertpos  = 0.015 + ( numplot_y - i_plot_v ) * dv ; position of colorbar. 
     1284colorbarn, horpos, vertpos, cb_length, 0.02, MIN(levels), MAX(levels), leveling, nlevels+1, levels, colors, format 
     1285 
     1286; 
     1287;-------------- 
     1288; aFe          
     1289;-------------- 
     1290; 
     1291i_plot = i_plot + 1 
     1292i_plot_h = i_plot MOD numplot_x 
     1293i_plot_v = i_plot / numplot_x + 1 
     1294 
     1295; contour parameters 
     1296nlevels = 10 
     1297nlevels = 10 
     1298fmin = 0.0 & fmax = MAX(aFeb) ; 
     1299cmin = 000 & cmax = 255 
     1300do_levels, fmin, fmax, cmin, cmax, nlevels, levels, colors 
     1301 
     1302ct = 3 ; colortable 
     1303title = ' aFe (mumol/m3)' 
     1304 
     1305icecontour, aFeb, h_i, h_s, nlevels, colors, levels, ct, z_ib, mon, nlay_bio, nts, $ 
     1306            tmin, tmax, zmin, zmax, title 
     1307 
     1308; colorbar 
     1309format='(f5.1)' 
     1310horpos   = 0.02 + i_plot_h * dh 
     1311vertpos  = 0.015 + ( numplot_y - i_plot_v ) * dv ; position of colorbar. 
     1312colorbarn, horpos, vertpos, cb_length, 0.02, MIN(levels), MAX(levels), leveling, nlevels+1, levels, colors, format 
     1313 
     1314; 
     1315;-------------- 
     1316; eFe          
     1317;-------------- 
     1318; 
     1319i_plot = i_plot + 1 
     1320i_plot_h = i_plot MOD numplot_x 
     1321i_plot_v = i_plot / numplot_x + 1 
     1322 
     1323; contour parameters 
     1324nlevels = 10 
     1325nlevels = 10 
     1326fmin = 0.0 & fmax = MAX(eFeb) ; 
     1327cmin = 000 & cmax = 255 
     1328do_levels, fmin, fmax, cmin, cmax, nlevels, levels, colors 
     1329 
     1330ct = 3 ; colortable 
     1331title = ' eFe (mumol/m3)' 
     1332 
     1333icecontour, eFeb, h_i, h_s, nlevels, colors, levels, ct, z_ib, mon, nlay_bio, nts, $ 
     1334            tmin, tmax, zmin, zmax, title 
     1335 
     1336; colorbar 
     1337format='(f5.1)' 
     1338horpos   = 0.02 + i_plot_h * dh 
     1339vertpos  = 0.015 + ( numplot_y - i_plot_v ) * dv ; position of colorbar. 
     1340colorbarn, horpos, vertpos, cb_length, 0.02, MIN(levels), MAX(levels), leveling, nlevels+1, levels, colors, format 
     1341 
     1342; 
     1343;------------------------------------------------------------------------------ 
     1344; End of the script 
     1345;------------------------------------------------------------------------------ 
     1346; 
     1347IF ( device EQ 'PS' ) THEN BEGIN 
     1348   DEVICE, /CLOSE 
     1349   SET_PLOT, "X" 
     1350   !P.MULTI=[0,2,2] 
     1351ENDIF 
    13221352END 
    13231353 

  • tags/LIM1D_v3.20/GRAPHICS/ALL_stocks.pro

    r6 r51  
    213213; YROSIAE 
    214214; taken from yrosiae ASPECT log sheets (Fripiat et al, oct 2015) 
     215   doy_obs    = [ 263 , 299,  318 , 279 , 306 , 327,  292 , 312  ,335 ] 
     216   Ichla_obs  = [ 1.4683  , 4.6038  ,45.4413, 128.0884,  44.5830 ,219.1416 ,208.6254 ,108.0657 , 87.9038 ] 
     217   ISi_obs    = [ 22.6631 , 21.1690 , 19.2174,  18.8949 , 18.9943,  18.9517,  20.5963,  16.6855,  13.8374 ] 
     218   IN_obs     = [ 9.0613  , 5.6130  , 9.4723,  10.0763,  11.4708 , 14.6193 , 14.4107 ,  4.7585 ,  2.1876 ] 
     219   IP_obs     = [ 0.4976  , 0.4067  , 1.1311,   1.8459,   2.3452 ,  3.7857 ,  3.6948 ,  1.3270 ,  1.2602 ] 
     220ENDIF 
     221 
     222IF ( site EQ 'N_ICE_RL' ) THEN BEGIN 
     223; N_ICE_RL 
     224; ALLLLLL WRONG !!! 
    215225   doy_obs    = [ 263 , 299,  318 , 279 , 306 , 327,  292 , 312  ,335 ] 
    216226   Ichla_obs  = [ 1.4683  , 4.6038  ,45.4413, 128.0884,  44.5830 ,219.1416 ,208.6254 ,108.0657 , 87.9038 ] 
    … …  
    492502   XYOUTS, xmin+ (xmax-xmin)*0.10, ymax*0.80, ztitle 
    493503 
    494    ; Iron stocks 
    495    ymin = 0.0000   & ymax = MAX(dFet)*1.2 
    496  
    497    PLOT, [xmin, xmax], [ ymin, ymax ], /NODATA, charsize = cs, $ 
    498          XTITLE = 'mon', YTITLE = 'Iron (mumol/m2)', XSTYLE = 1, YSTYLE = 1 
    499  
    500    FOR i_run = nruns - 1, 0, -1 DO BEGIN 
    501       LOADCT, ct 
    502       OPLOT, mon(i_run,0:ntss(i_run)-1), dFet(0:ntss(i_run)-1), color = 150, thick = 2 
    503       OPLOT, mon(i_run,0:ntss(i_run)-1), aFet(0:ntss(i_run)-1), color = 0, thick = 2 
    504       OPLOT, mon(i_run,0:ntss(i_run)-1), eFet(0:ntss(i_run)-1), color = 0, thick = 1 
    505    ENDFOR 
    506    LOADCT, 0 
    507  
    508    ztitle = 'Max dFe stock: '+STRCOMPRESS(STRING(MAX(dFet), FORMAT='(F6.2)'),/REMOVE_ALL)+' mumol/m2' 
    509    XYOUTS, xmin+ (xmax-xmin)*0.10, ymax*0.90, ztitle 
     504;  ; Iron stocks 
     505;  ymin = 0.0000   & ymax = 500. ;ymax = MAX(dFet)*1.2 
     506 
     507;  PLOT, [xmin, xmax], [ ymin, ymax ], /NODATA, charsize = cs, $ 
     508;        XTITLE = 'mon', YTITLE = 'Iron (mumol/m2)', XSTYLE = 1, YSTYLE = 1 
     509 
     510;  FOR i_run = nruns - 1, 0, -1 DO BEGIN 
     511;     LOADCT, ct 
     512;     OPLOT, mon(i_run,0:ntss(i_run)-1), dFet(0:ntss(i_run)-1), color = 150, thick = 2 
     513;     OPLOT, mon(i_run,0:ntss(i_run)-1), aFet(0:ntss(i_run)-1), color = 0, thick = 2 
     514;     OPLOT, mon(i_run,0:ntss(i_run)-1), eFet(0:ntss(i_run)-1), color = 0, thick = 1 
     515;  ENDFOR 
     516;  LOADCT, 0 
     517 
     518;  ztitle = 'Max dFe stock: '+STRCOMPRESS(STRING(MAX(dFet), FORMAT='(F6.2)'),/REMOVE_ALL)+' mumol/m2' 
     519;  XYOUTS, xmin+ (xmax-xmin)*0.10, ymax*0.90, ztitle 
    510520 
    511521 

  • tags/LIM1D_v3.20/INPUT/N_ICE_RL/source_3.20/bio.param

    r49 r51  
    116116# 
    117117# nn_phs      : type of photosynthesis parameterization (1-4) 
    118 4 
     1183 
    119119# ek_bio      : Light adaptation parameter (nn_phs = 1, muE m-2 s-1) - DEFLT 2 
    1201202 

  • tags/LIM1D_v3.20/INPUT/N_ICE_RL/source_3.20/forcing.param.1h

    r49 r51  
    323299      1.0     0.0 
    3333# oceanic heat flux (W/m2) ; 5. is the ref 
    34 0       1.0     5. 
     3499      1.0     40. 
    3535# snowfall (dimensions=length/timestep) 0.012 - 1h 
    363699      1.0     0.000 

  • tags/LIM1D_v3.20/INPUT/N_ICE_RL/source_3.20/icephys.param

    r49 r51  
    5757# 
    5858# emig          ! Surface emissivity (past 0.995) 
    59 0.995 
     590.95 
    6060# fpar_fsw      ! Visible energetic ratio (Fpar/Fsw, J/J) = 0.475 (Frouin Pinker 95) 
    61610.475 

  • tags/LIM1D_v3.20/INPUT/N_ICE_RL/source_3.20/tracer.param

    r49 r51  
    2323   'dIN'   'nox'   3   10.44  26.90    0.       0. 
    2424   'dIP'   'po4'   3    0.91   2.35    0.       0. 
    25    'AoC'   'AoC'   3    0.01   0.05    0.       0. 
    26    'eoC'   'eoC'   3      0.    0.5    0.       0. 
     25   'AoC'   'AoC'   1    0.01   0.05    0.       0. 
     26   'eoC'   'eoC'   1      0.    0.5    0.       0. 
    2727   'DIC'   'DIC'   1    800.   2000.   0.       0. 
    2828   'Alk'   'Alk'   1    880.   2200.   0.       0. 

  • tags/LIM1D_v3.20/SCRIPTS/lim1d_n_ice_rl.nqs

    r50 r51  
    1717SUBDIR=`pwd` 
    1818 
    19 MAINDIR="$HOME/...yourdir..."                 ### CHANGE 
     19MAINDIR="/net/hera/usr/hera/locean/temp/mvlod/2017/LIM1D_v3.20_BEPSII"                 ### CHANGE 
    2020 
    2121SCRATCHDIR="$MAINDIR/SCRATCH/$EXP_ID"        # SCRATCH - temporary directory on which code is run 
    … …  
    2424INFILEDIR="$MAINDIR/INPUT/$CONF/$SOURCE"     # INPUT - initialization, forcing & param files 
    2525 
    26 GRAPHDIR="$HOME/Boulot/SCIENCE/PLOT_SCRIPTS/LIM1D_BIO/IDL/$CONF" # IDL plots 
     26GRAPHDIR="$MAINDIR/GRAPHICS" # IDL plots 
    2727 
    2828####################################################################################### 
    … …  
    5454# rename files 
    5555mv "forcing_"$CONF"_"$TS".nc"      forcing.nc 
     56  forcing_N_ICE_1h.nc 
    5657mv forcing.param.$TS               forcing.param 
    5758mv run.param.$TS                   run.param 
    … …  
    114115echo 
    115116cat exp_id.dat 
    116 cp exp_id.dat .. 
     117#cp exp_id.dat .. 
    117118echo 
    118119echo  
Note: See TracChangeset for help on using the changeset viewer.