Skip to content
Snippets Groups Projects
Commit f19be1f6 authored by Victor Stepanenko's avatar Victor Stepanenko
Browse files

Minor improvements

parent 77525ac3
No related branches found
No related tags found
No related merge requests found
...@@ -213,7 +213,8 @@ use PHYS_FUNC, only: & ...@@ -213,7 +213,8 @@ use PHYS_FUNC, only: &
& UNFRWAT, & & UNFRWAT, &
& WI_MAX, & & WI_MAX, &
& WL_MAX, & & WL_MAX, &
& MIXED_LAYER_CALC & MIXED_LAYER_CALC, &
& HC_CORR_CARBEQUIL
use EVOLUTION_VARIABLES, only : & use EVOLUTION_VARIABLES, only : &
& UPDATE_CURRENT_TIMESTEP, & & UPDATE_CURRENT_TIMESTEP, &
...@@ -228,7 +229,9 @@ use METH_OXYG_CONSTANTS, only : & ...@@ -228,7 +229,9 @@ use METH_OXYG_CONSTANTS, only : &
& molm3tonM, & & molm3tonM, &
& molm3toppm_co2,molm3toppm_o2,molm3toppm_ch4, & & molm3toppm_co2,molm3toppm_o2,molm3toppm_ch4, &
& molm3tomgl_co2,molm3tomkgl_ch4,molm3tomgl_o2, & & molm3tomgl_co2,molm3tomkgl_ch4,molm3tomgl_o2, &
& photic_threshold & molm3tomcM, &
& photic_threshold, &
& pH
use OXYGEN_MOD, only : & use OXYGEN_MOD, only : &
& OXYGEN, OXYGEN_PRODCONS, & & OXYGEN, OXYGEN_PRODCONS, &
...@@ -2407,7 +2410,7 @@ if (flag_print) then ! The output in ASCII files ...@@ -2407,7 +2410,7 @@ if (flag_print) then ! The output in ASCII files
& outpath, 'oxygen_water', i, ndec, .false.) & outpath, 'oxygen_water', i, ndec, .false.)
ndec = -1 ndec = -1
i = i + 2 i = i + 2
work1(1:M+1) = qwater(1:M+1,1)*molm3tomkgl_ch4 ! molm3tonM work1(1:M+1) = qwater(1:M+1,1)*molm3tomcM !molm3tomkgl_ch4 ! molm3tonM
call PROFILE_OUTPUT & call PROFILE_OUTPUT &
& (ix, iy, dnx, dny, & & (ix, iy, dnx, dny, &
& year, month, day, hour, & & year, month, day, hour, &
...@@ -2430,7 +2433,9 @@ if (flag_print) then ! The output in ASCII files ...@@ -2430,7 +2433,9 @@ if (flag_print) then ! The output in ASCII files
endif endif
ndec = -1 ndec = -1
i = i + 2 i = i + 2
work1(1:M+1) = DIC(1:M+1,1)*molm3tomgl_co2 !molm3tonM do j = 1, M+1
work1(j) = DIC(j,1) / HC_CORR_CARBEQUIL(Tw1(j)+Kelvin0,pH) * molm3tomcM !molm3tomgl_co2 !molm3tonM
enddo
call PROFILE_OUTPUT & call PROFILE_OUTPUT &
& (ix, iy, dnx, dny, & & (ix, iy, dnx, dny, &
& year, month, day, hour, & & year, month, day, hour, &
......
...@@ -2,6 +2,7 @@ MODULE OUT_MOD ...@@ -2,6 +2,7 @@ MODULE OUT_MOD
use LAKE_DATATYPES, only : ireals, iintegers use LAKE_DATATYPES, only : ireals, iintegers
use INOUT, only : CHECK_UNIT use INOUT, only : CHECK_UNIT
use DRIVING_PARAMS, only : missing_value
character, save :: outpath*60 character, save :: outpath*60
...@@ -1475,6 +1476,7 @@ endif ...@@ -1475,6 +1476,7 @@ endif
if (int(time/(dt_out*hour_sec)) > tcount) then if (int(time/(dt_out*hour_sec)) > tcount) then
allocate (Temp(1:rtemp%numarr)) allocate (Temp(1:rtemp%numarr))
Temp(:) = missing_value
do k = 1, rtemp%numarr do k = 1, rtemp%numarr
sumh = 0. sumh = 0.
...@@ -1489,6 +1491,7 @@ if (int(time/(dt_out*hour_sec)) > tcount) then ...@@ -1489,6 +1491,7 @@ if (int(time/(dt_out*hour_sec)) > tcount) then
larthick = h*sum(gsp%ddz05(itherm(i):itherm(i+1)-1)) ! Equilibrium thickness of constant-density layer larthick = h*sum(gsp%ddz05(itherm(i):itherm(i+1)-1)) ! Equilibrium thickness of constant-density layer
dlarthick = amplx*sin(pi*rtemp%arr1(1,k)/Lx) + amply*sin(pi*rtemp%arr1(2,k)/Ly) dlarthick = amplx*sin(pi*rtemp%arr1(1,k)/Lx) + amply*sin(pi*rtemp%arr1(2,k)/Ly)
sumh = sumh + larthick + dlarthick sumh = sumh + larthick + dlarthick
!print*, 'ay', hy1ml
endif if1 endif if1
if (sumh >= h - rtemp%arr1(3,k)) then if (sumh >= h - rtemp%arr1(3,k)) then
!print*, itherm(i)+1, amplx !print*, itherm(i)+1, amplx
...@@ -1499,12 +1502,14 @@ if (int(time/(dt_out*hour_sec)) > tcount) then ...@@ -1499,12 +1502,14 @@ if (int(time/(dt_out*hour_sec)) > tcount) then
if ( h_ <= sumh1 .and. h_ > sumh1 - 0.5*dsumh ) then if ( h_ <= sumh1 .and. h_ > sumh1 - 0.5*dsumh ) then
alpha = (h_ - sumh1 + 0.5*dsumh)/dsumh !Assuming regular grid and the same thickness deviation alpha = (h_ - sumh1 + 0.5*dsumh)/dsumh !Assuming regular grid and the same thickness deviation
!in adjacent layers !in adjacent layers
Temp(k) = Tw(j-1)*alpha + Tw(j)*(1. - alpha) Temp(k) = Tw(max(j-1,1))*alpha + Tw(j)*(1. - alpha)
!print*, '1', j, nvar
exit cyclrs exit cyclrs
elseif ( h_ <= sumh1 - 0.5*dsumh .and. h_ > sumh1 - dsumh ) then elseif ( h_ <= sumh1 - 0.5*dsumh .and. h_ > sumh1 - dsumh ) then
alpha = (sumh1 - 0.5*dsumh - h_)/dsumh !Assuming regular grid and the same thickness deviation alpha = (sumh1 - 0.5*dsumh - h_)/dsumh !Assuming regular grid and the same thickness deviation
!in adjacent layers !in adjacent layers
Temp(k) = Tw(j+1)*alpha + Tw(j)*(1. - alpha) Temp(k) = Tw(min(j+1,M+1))*alpha + Tw(j)*(1. - alpha)
!print*, '2', j, nvar
exit cyclrs exit cyclrs
endif endif
sumh1 = sumh1 - dsumh sumh1 = sumh1 - dsumh
......
...@@ -902,6 +902,9 @@ END FUNCTION FP_THETA ...@@ -902,6 +902,9 @@ END FUNCTION FP_THETA
! Thermocline thickness ! Thermocline thickness
x1 = 0.; x2 = ls%h1 x1 = 0.; x2 = ls%h1
u1 = 0.; u2 = 0.
v1 = 0.; v2 = 0.
Sal1 = 0.; Sal2 = 0.
do i = 1, M do i = 1, M
if (wst%Tw2(i) >= Ttherm0 .and. wst%Tw2(i+1) < Ttherm0) then if (wst%Tw2(i) >= Ttherm0 .and. wst%Tw2(i+1) < Ttherm0) then
y = gsp%ddz(i)*ls%h1 y = gsp%ddz(i)*ls%h1
...@@ -1476,6 +1479,8 @@ END FUNCTION FP_THETA ...@@ -1476,6 +1479,8 @@ END FUNCTION FP_THETA
!! \f[ !! \f[
!! HC\_CORR\_CARBEQUIL = (1+k_1 10^{pH}+k_1*k_2 10^{2pH}) !! HC\_CORR\_CARBEQUIL = (1+k_1 10^{pH}+k_1*k_2 10^{2pH})
!! \f] !! \f]
!! Concomitantly, it is a ratio of DIC to CO_2 molar concenatrations.
!! DIC molar concentration is a molar concentration of C atoms in DIC.
FUNCTION HC_CORR_CARBEQUIL(temp,pH) FUNCTION HC_CORR_CARBEQUIL(temp,pH)
implicit none implicit none
......
...@@ -30,7 +30,7 @@ nplot = input("Choose a variable to plot \n \ ...@@ -30,7 +30,7 @@ nplot = input("Choose a variable to plot \n \
nmod_max = 12 #The last item in the list above, from the modeled data nmod_max = 12 #The last item in the list above, from the modeled data
path = '../results/Mojai2016/time_series/' # path to files path = '../results/Kuiva2014septnew/time_series/' # path to files
#pathobs = os.path.expanduser('/media/victor/main/Files/main/observdata/Iseo/') #pathobs = os.path.expanduser('/media/victor/main/Files/main/observdata/Iseo/')
pathobs = os.path.expanduser('/media/victor/main/Files/main/observdata/Kuivajarvi/Kuiva2014sept_campaign/') pathobs = os.path.expanduser('/media/victor/main/Files/main/observdata/Kuivajarvi/Kuiva2014sept_campaign/')
exp = 'base' exp = 'base'
...@@ -77,8 +77,8 @@ if int(nplot) == 1: # Temperature ...@@ -77,8 +77,8 @@ if int(nplot) == 1: # Temperature
#t1disp = 7.95 # Final time for display #t1disp = 7.95 # Final time for display
#nvardisp = nlevs_water #nvardisp = nlevs_water
title = "Temperature, $^\circ$C " title = "Temperature, $^\circ$C "
levmin = 6. levmin = 8.
levmax = 30. levmax = 20.
nlevplot = 20 # Number of countour levels nlevplot = 20 # Number of countour levels
pathsave = path pathsave = path
if int(nplot) == 2: # Salinity if int(nplot) == 2: # Salinity
...@@ -116,12 +116,12 @@ elif int(nplot) == 3: # Methane ...@@ -116,12 +116,12 @@ elif int(nplot) == 3: # Methane
#t0disp = 6 # Initial time for display #t0disp = 6 # Initial time for display
#t1disp = 11 # Final time for display #t1disp = 11 # Final time for display
#nvardisp = nlevs_water #nvardisp = nlevs_water
title = "Methane, $\mu g/l$" title = "Methane, $mmol/m^3$"
text = 'maxval' text = 'maxval'
#levmin = 0. #levmin = 0.
#levmax = 600. #levmax = 5000.
nlevplot = 20 # Number of countour levels nlevplot = 20 # Number of countour levels
stepticks = 60 #stepticks = 60
pathsave = path pathsave = path
elif int(nplot) == 4: # Oxygen elif int(nplot) == 4: # Oxygen
basename = 'oxygen_water 1 1' basename = 'oxygen_water 1 1'
...@@ -156,11 +156,11 @@ elif int(nplot) == 5: # Carbon dioxide ...@@ -156,11 +156,11 @@ elif int(nplot) == 5: # Carbon dioxide
#t0disp = 6 # Initial time for display #t0disp = 6 # Initial time for display
#t1disp = 11 # Final time for display #t1disp = 11 # Final time for display
#nvardisp = nlevs_water #nvardisp = nlevs_water
title = "Carbon dioxide, mg/l" title = "Carbon dioxide, $mmol/m^3$"
#levmin = 0. #levmin = 0.
#levmax = 20. #levmax = 20.
nlevplot = 20 # Number of countour levels nlevplot = 20 # Number of countour levels
stepticks = 2 #stepticks = 2
pathsave = path pathsave = path
elif int(nplot) == 6: # Methane bubble flux elif int(nplot) == 6: # Methane bubble flux
basename = 'fbblflx_ch4 1 1' basename = 'fbblflx_ch4 1 1'
...@@ -266,7 +266,8 @@ elif int(nplot) == 11: # y-speed component ...@@ -266,7 +266,8 @@ elif int(nplot) == 11: # y-speed component
time0 = month_doy(7,ny)/mdays + 17./31. time0 = month_doy(7,ny)/mdays + 17./31.
#t0disp = 4 # Initial time for display #t0disp = 4 # Initial time for display
#t1disp = 10 # Final time for display #t1disp = 10 # Final time for display
#levmax = 1. #levmin = -0.01
#levmax = 0.01
#nvardisp = nlevs_water #nvardisp = nlevs_water
title = "y-speed component, $m/s$" title = "y-speed component, $m/s$"
pathsave = path pathsave = path
......
...@@ -159,10 +159,10 @@ def plot_timeser(plotid,nvar,legendcommon) : ...@@ -159,10 +159,10 @@ def plot_timeser(plotid,nvar,legendcommon) :
#path = '../results/Seida1_2007-2008/time_series/' # path to files #path = '../results/Seida1_2007-2008/time_series/' # path to files
# Number of experiments is limited to 10 # Number of experiments is limited to 10
pathobs = [os.path.expanduser('/home/victor/Files/main/observdata/Mojai2016/')] #pathobs = [os.path.expanduser('/home/victor/Files/main/observdata/Mojai2016/')]
#pathobs = [os.path.expanduser('~/Files/main/observdata/Valkea-Kotinen2006/')] pathobs = [os.path.expanduser('~/Files/main/observdata/Kuivajarvi/campaign2014/')]
#path_ = ['../results/Kuivajarvi2013_c2/time_series/']#, \ #path_ = ['../results/Kuivajarvi2013_c2/time_series/']#, \
path_ = ['../results/Mojai2016/time_series/'] #, \ path_ = ['../results/Kuiva2014septnew/time_series/'] #, \
#'../results/Valkea-Kotinen2006_nomomflxpart/time_series/', \ #'../results/Valkea-Kotinen2006_nomomflxpart/time_series/', \
#'../results/Valkea-Kotinen2006_nokor/time_series/', \ #'../results/Valkea-Kotinen2006_nokor/time_series/', \
#'../results/Valkea-Kotinen2006_dynpgrad4nokor/time_series/', \ #'../results/Valkea-Kotinen2006_dynpgrad4nokor/time_series/', \
...@@ -191,11 +191,11 @@ def plot_timeser(plotid,nvar,legendcommon) : ...@@ -191,11 +191,11 @@ def plot_timeser(plotid,nvar,legendcommon) :
nlevs_water = 16 nlevs_water = 16
#The time of the first data entry in time series #The time of the first data entry in time series
ny = 2016 #year ny = 2014 #year
mdays = 365/12 mdays = 365/12
nmonth = 6 nmonth = 9
nday = 25 nday = 2
nhour = 1 nhour = 0
nmin = 0 nmin = 0
timezone = None timezone = None
...@@ -253,12 +253,12 @@ def plot_timeser(plotid,nvar,legendcommon) : ...@@ -253,12 +253,12 @@ def plot_timeser(plotid,nvar,legendcommon) :
nameexp = nameexp_ nameexp = nameexp_
npaths = len(path_) npaths = len(path_)
if int(nplot) == tempmod3_: # Temperature at selected points if int(nplot) == tempmod3_: # Temperature at selected points
basename = 'temploc 1 1' basename = 'temploc 1 1' #'temploc 1 1'
fileout = 'temploc' fileout = 'temploc'
nheader = 0 nheader = 0
ncols = 3 #Total number of columns ncols = 3 #Total number of columns
ncol = [3] #[2,3] ncol = [2,3] #[2,3]
legend = ['Southern point'] #['point 1','point 2'] legend = ['T$_{7m}$, model','T$_{11m}$, model'] #['point 1','point 2']
labelx = 'Time, days' labelx = 'Time, days'
labely = 'Temperature, C$^{\circ}$' labely = 'Temperature, C$^{\circ}$'
timescale = 24. # 365/12*24 - number of hours in a month timescale = 24. # 365/12*24 - number of hours in a month
...@@ -268,6 +268,7 @@ def plot_timeser(plotid,nvar,legendcommon) : ...@@ -268,6 +268,7 @@ def plot_timeser(plotid,nvar,legendcommon) :
path = path_ path = path_
nameexp = nameexp_ nameexp = nameexp_
npaths = len(path_) npaths = len(path_)
missingval = -999.
if int(nplot) == senslat_: # Sensible and latent heat fluxes if int(nplot) == senslat_: # Sensible and latent heat fluxes
basename = 'T_fluxes 1 1' basename = 'T_fluxes 1 1'
fileout = 'Senslat_ser' fileout = 'Senslat_ser'
...@@ -712,14 +713,36 @@ def plot_timeser(plotid,nvar,legendcommon) : ...@@ -712,14 +713,36 @@ def plot_timeser(plotid,nvar,legendcommon) :
# nameexp = nameexpobs # nameexp = nameexpobs
# npaths = 1 # npaths = 1
# delim = ',' # delim = ','
#if int(nplot) == tempobs_: # Temperature measured
# #Mojai2016
# basename = 'Tsurf_obs_Gella_2016' #'LDS'
# fileout = 'Temploc_obs'
# nheader = 0 #Number of lines in the header
# ncols = 3 # LDS: 7 + 21*2 #Total number of columns
# ncol = [3] #LDS_14.8m:[37]
# legend = ['T\_surf'] #['T\_14.80m']
# #ncol = nm.arange(9, ncols, 2) #,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17] # ,20,21,22,23] #The number of columns in the file to display
# #legend = \
# #['T_0.30m', 'T_0.55m', 'T_0.80m', 'T_1.05m', 'T_1.30m', 'T_2.30m', 'T_3.30m', \
# # 'T_4.30m', 'T_5.80m', 'T_7.30m', 'T_8.80m', 'T_10.30m', 'T_11.80m', 'T_13.30m', \
# # 'T_14.80m', 'T_16.30m', 'T_17.80m', 'T_19.80m', 'T_21.80m', 'T_35.80m', 'T_49.80m'] #LDS
# labelx = 'Time, days'
# labely = 'Temperature, Celsius'
# ctime = 2
# timescale = 1. # 365/12*24 - number of hours in a month
# time0 = 0. #month_doy(7,ny)/mdays + 17./31. #Initial time
# path = pathobs
# nameexp = nameexpobs
# npaths = 1
# delim = ','
if int(nplot) == tempobs_: # Temperature measured if int(nplot) == tempobs_: # Temperature measured
#Mojai2016 #Kuivajarvi2014 campaign September
basename = 'Tsurf_obs_Gella_2016' #'LDS' basename = 'KuivaMETEO_201409' #'LDS'
fileout = 'Temploc_obs' fileout = 'Temploc_obs'
nheader = 0 #Number of lines in the header nheader = 0 #Number of lines in the header
ncols = 3 # LDS: 7 + 21*2 #Total number of columns ncols = 34 # LDS: 7 + 21*2 #Total number of columns
ncol = [3] #LDS_14.8m:[37] ncol = [15,17] #LDS_14.8m:[37]
legend = ['T\_surf'] #['T\_14.80m'] legend = ['T$_{7m}$','T$_{10m}$'] #['T\_14.80m']
#ncol = nm.arange(9, ncols, 2) #,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17] # ,20,21,22,23] #The number of columns in the file to display #ncol = nm.arange(9, ncols, 2) #,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17] # ,20,21,22,23] #The number of columns in the file to display
#legend = \ #legend = \
#['T_0.30m', 'T_0.55m', 'T_0.80m', 'T_1.05m', 'T_1.30m', 'T_2.30m', 'T_3.30m', \ #['T_0.30m', 'T_0.55m', 'T_0.80m', 'T_1.05m', 'T_1.30m', 'T_2.30m', 'T_3.30m', \
...@@ -733,7 +756,8 @@ def plot_timeser(plotid,nvar,legendcommon) : ...@@ -733,7 +756,8 @@ def plot_timeser(plotid,nvar,legendcommon) :
path = pathobs path = pathobs
nameexp = nameexpobs nameexp = nameexpobs
npaths = 1 npaths = 1
delim = ',' #delim = ','
missingval = -999.
if int(nplot) == mldobs_: # Mixed-layer depth observed if int(nplot) == mldobs_: # Mixed-layer depth observed
basename = 'Valkea-Kotinen_validation_May-Dec_2006_3_mod' #'KJdata2013' #'obs' basename = 'Valkea-Kotinen_validation_May-Dec_2006_3_mod' #'KJdata2013' #'obs'
fileout = 'mld_obs' fileout = 'mld_obs'
...@@ -783,21 +807,44 @@ def plot_timeser(plotid,nvar,legendcommon) : ...@@ -783,21 +807,44 @@ def plot_timeser(plotid,nvar,legendcommon) :
path = pathobs path = pathobs
nameexp = nameexpobs nameexp = nameexpobs
npaths = 1 npaths = 1
#if int(nplot) == co2obs1_: # Carbon dioxide measured
# basename = 'KJdata2013'
# fileout = 'co2_obs'
# nheader = 1 #Number of lines in the header
# ncols = 21 #Total number of columns
# ncol = [18]#,19,20,21] #The number of columns in the file to display
# legend = ['CO2_0.2m']#,'CO2_1.5m','CO2_2.5m','CO2_7m']
# labelx = 'Time, months'
# labely = 'Carbon dioxide, ppm'
# ctime = 1
# timescale = mdays # 365/12*24 - number of hours in a month
# time0 = month_doy(1,ny)/mdays #Initial time
# path = pathobs
# nameexp = nameexpobs
# npaths = 1
if int(nplot) == co2obs1_: # Carbon dioxide measured if int(nplot) == co2obs1_: # Carbon dioxide measured
basename = 'KJdata2013' # Kuivajarvi 2014 campaign september
basename = 'KuivaMETEO_201409' #'LDS'
fileout = 'co2_obs' fileout = 'co2_obs'
nheader = 1 #Number of lines in the header nheader = 0 #Number of lines in the header
ncols = 21 #Total number of columns ncols = 34 # LDS: 7 + 21*2 #Total number of columns
ncol = [18]#,19,20,21] #The number of columns in the file to display ncol = [30] #LDS_14.8m:[37]
legend = ['CO2_0.2m']#,'CO2_1.5m','CO2_2.5m','CO2_7m'] legend = ['CO$_2$, 7 m'] #['T\_14.80m']
labelx = 'Time, months' #ncol = nm.arange(9, ncols, 2) #,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17] # ,20,21,22,23] #The number of columns in the file to display
labely = 'Carbon dioxide, ppm' #legend = \
ctime = 1 #['T_0.30m', 'T_0.55m', 'T_0.80m', 'T_1.05m', 'T_1.30m', 'T_2.30m', 'T_3.30m', \
timescale = mdays # 365/12*24 - number of hours in a month # 'T_4.30m', 'T_5.80m', 'T_7.30m', 'T_8.80m', 'T_10.30m', 'T_11.80m', 'T_13.30m', \
time0 = month_doy(1,ny)/mdays #Initial time # 'T_14.80m', 'T_16.30m', 'T_17.80m', 'T_19.80m', 'T_21.80m', 'T_35.80m', 'T_49.80m'] #LDS
labelx = 'Time, days'
labely = 'Concentration, ppm'
ctime = 2
timescale = 1. # 365/12*24 - number of hours in a month
time0 = 0. #month_doy(7,ny)/mdays + 17./31. #Initial time
path = pathobs path = pathobs
nameexp = nameexpobs nameexp = nameexpobs
npaths = 1 npaths = 1
#delim = ','
missingval = -999.
if int(nplot) == co2obs2_: # CO_2 measured manually if int(nplot) == co2obs2_: # CO_2 measured manually
basename = 'ch4&co2man2013' basename = 'ch4&co2man2013'
fileout = 'co2_obs_man' fileout = 'co2_obs_man'
...@@ -1038,12 +1085,13 @@ def plot_timeser(plotid,nvar,legendcommon) : ...@@ -1038,12 +1085,13 @@ def plot_timeser(plotid,nvar,legendcommon) :
line = line.split(delim) line = line.split(delim)
else: else:
line = line.split() line = line.split()
print(line)
time[j] = float(line[ctime-1]) time[j] = float(line[ctime-1])
for i in range(len(ncol)): for i in range(len(ncol)):
vals[j,i] = float(line[ncol[i]]) vals[j,i] = float(line[ncol[i]])
if time_label_date: time2 = time[:]*time_unit + time0_
time[:] = time[:]/timescale + time0 time[:] = time[:]/timescale + time0
if time_label_date: time2 = time[:]*time_unit + time0_
for j in range(nvardisp) : for j in range(nvardisp) :
vard = [0]*nrows vard = [0]*nrows
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment