Log of /trunk/phylmd/Interface_surf/interfsurf_hq.f
Parent Directory
| Revision Log
Revision
301 -
(
view)
(
annotate)
-
[select for diffs]
Modified
Thu Aug 2 17:23:07 2018 UTC
(5 years, 11 months ago)
by
guez
File length: 8241 byte(s)
Diff to
previous 300
,
to
selected 217
Move the call to conf_interface up to physiq, so there is no need to
test first call inside pbl_surface for this.
run_off_lic in fonte_neige was computed but not used. Pass it up to
pbl_surface so we can output it (following LMDZ).
Simplify the logic in interfsur_lim so we do not need debut.
Remove the tests on the order of surface types in interfsurf_hq. Just
add comments in indicesol.
Revision
300 -
(
view)
(
annotate)
-
[select for diffs]
Modified
Thu Aug 2 15:55:01 2018 UTC
(5 years, 11 months ago)
by
guez
File length: 9166 byte(s)
Diff to
previous 299
,
to
selected 217
The calls to calcul_fluxs were always done with an array argument
dif_grnd set to a scalar (and that is also the case in LMDZ). So just
use a scalar argument.
In procedure fonte_neige, the value of run_off_lic from previous call
was actually not used. So we can remove the save attribute and make it
an automatic array.
Revision
299 -
(
view)
(
annotate)
-
[select for diffs]
Modified
Thu Aug 2 14:27:11 2018 UTC
(5 years, 11 months ago)
by
guez
File length: 9602 byte(s)
Diff to
previous 297
,
to
selected 217
Use directly dtphys from module comconst when possible instead of
having it trickle down through procedure arguments.
Revision
297 -
(
view)
(
annotate)
-
[select for diffs]
Modified
Thu Jul 26 16:02:11 2018 UTC
(5 years, 11 months ago)
by
guez
File length: 9835 byte(s)
Diff to
previous 293
,
to
selected 217
Rename module interface_surf to conf_interface_m.
Move the computation of pkf out of procedure climb_hq_down into clqh.
Revision
280 -
(
view)
(
annotate)
-
[select for diffs]
Modified
Fri Jul 20 15:47:57 2018 UTC
(5 years, 11 months ago)
by
guez
File length: 10100 byte(s)
Diff to
previous 279
,
to
selected 217
Remove gamq (counter-gradient for water vapor) in procedure clqh, was
always 0 (does not appear either any longer in LMDZ).
Revision
279 -
(
view)
(
annotate)
-
[select for diffs]
Modified
Fri Jul 20 14:30:23 2018 UTC
(5 years, 11 months ago)
by
guez
File length: 10135 byte(s)
Diff to
previous 269
,
to
selected 217
fqcalving was saved in physiq and had intent inout in pbl_surface. So
we could set fqcalving to 0 only once per run. The point is fqcalving
must be defined everywhere for the computation of the average over all
surfaces, even values that get multiplied by pctsrf = 0. I find it
clearer to set fqcalving to 0 at every call of pbl_surface. This is
more expensive but allows to give intent out to fqcalving in
pbl_surface and remove the save attribute in physiq.
Add zxfqcalving output netCDF variable (following LMDZ).
Revision
269 -
(
view)
(
annotate)
-
[select for diffs]
Modified
Thu May 3 17:14:14 2018 UTC
(6 years, 2 months ago)
by
guez
File length: 10168 byte(s)
Diff to
previous 268
,
to
selected 217
Bug fix: dif_grnd for case is_sic in procedure
interfsurf_hq. (Following LMDZ. In LMDZ, the bug was fixed in revision
686.)
Revision
267 -
(
view)
(
annotate)
-
[select for diffs]
Modified
Thu May 3 16:14:08 2018 UTC
(6 years, 2 months ago)
by
guez
File length: 10268 byte(s)
Diff to
previous 254
,
to
selected 217
Rename procedure clmain to pbl_surface (following LMDZ).
Remove choice soil_model = f. This choice made the algorithm unclear
in interfsurf_hq. Also soil_model = f is never used in LMDZ. radsol
was intent inout in clqh because of its possible modification in
interfsurf_hq, but the corresponding actual argument yrads is not used
in pbl_surface. The modification of radsol in interfsurf_hq was a bad
idea. Now we can more clearly make radsol an intent in argument of
clqh and interfsurf_hq.
Revision
225 -
(
view)
(
annotate)
-
[select for diffs]
Modified
Mon Oct 16 12:35:41 2017 UTC
(6 years, 8 months ago)
by
guez
Original Path:
trunk/Sources/phylmd/Interface_surf/interfsurf_hq.f
File length: 11001 byte(s)
Diff to
previous 223
,
to
selected 217
LMDZE is now in Fortran 2003 (use of allocatable arguments).
gradsdef was not used.
Change names: [uv]10m to [uv]10m_srf in clmain, y[uv]1 to
[uv]1lay. Remove useless complication: zx_alf[12]. Do not modify
[uv]1lay after initial definition from [uv].
Add [uv]10m_srf to output.
Change names in physiq: [uv]10m to [uv]10m_srf, z[uv]10m to [uv]10m,
corresponding to NetCDF output names.
Remove unused complication couchelimite and useless variable inirnpb
in phytrac.
Revision
221 -
(
view)
(
annotate)
-
[select for diffs]
Modified
Thu Apr 20 14:44:47 2017 UTC
(7 years, 2 months ago)
by
guez
Original Path:
trunk/Sources/phylmd/Interface_surf/interfsurf_hq.f
File length: 11497 byte(s)
Diff to
previous 217
clcdrag is no longer used in LMDZ. Replaced by cdrag in LMDZ. In cdrag
in LMDZ, zxli is a symbolic constant, false. So removed case zxli true
in LMDZE.
read_sst is called zero (if no ocean point on the whole planet) time or
once per call of physiq. If mod(itap - 1, lmt_pas) == 0 then we have
advanced in time of lmt_pas and deja_lu is necessarily false.
qsat[sl] and dqsat[sl] were never called.
Added output of qsurf in histins, following LMDZ.
Last dummy argument dtime of phystokenc is always the same as first
dummy argument pdtphys, removed dtime.
Removed make rules for nag_xref95, since it does not exist any longer.
Revision
217 -
(
view)
(
annotate)
-
[selected]
Modified
Thu Mar 30 14:25:18 2017 UTC
(7 years, 3 months ago)
by
guez
Original Path:
trunk/Sources/phylmd/Interface_surf/interfsurf_hq.f
File length: 11530 byte(s)
Diff to
previous 215
run_off_lic downgraded from variable of module interface_surf to local
variable of fonte_neige.
Code could not work with ok_aie set to true, so removed this
possibility. tauae, piz_ae, cg_ae, topswai, solswai were then
0. cldtaupi was the same as cldtaupd.
In sw and procedures called by sw, flag_aer did not need to be double
precision, changed it to logical.
Downgraded re and fl from arguments of newmicro to local
variables. Added output of re and fl (following LMDZ).
Revision
215 -
(
view)
(
annotate)
-
[select for diffs]
Modified
Tue Mar 28 12:46:28 2017 UTC
(7 years, 3 months ago)
by
guez
Original Path:
trunk/Sources/phylmd/Interface_surf/interfsurf_hq.f
File length: 11668 byte(s)
Diff to
previous 214
,
to
selected 217
size(snow) is now knon in interfsurf_hq.
Renamed snow to fsnow in clmain, same name as corresponding actual
argument. We can then rename ysnow to simply snow in clmain, same name
as corresponding dummy argument of clqh. No need to initialize local
snow to 0 since it is only used with indices 1:knon and already
initialized from fsnow for each type of surface. If there is no point
for a given type of surface, fsnow should be reset to 0 for this
type. We need to give a valid value to fsnow in this case even if it
will be multiplied by pctsrf = 0 in physiq.
In physiq, no need for intermediate zxsnow for output.
Removed unused arguments tsurf, p1lay, beta, coef1lay, ps, t1lay,
q1lay, u1lay, v1lay, petAcoef, peqAcoef, petBcoef, peqBcoef of
fonte_neige, with unused computations of zx_qs and zcor. (Same was
done in LMDZ.)
Revision
209 -
(
view)
(
annotate)
-
[select for diffs]
Modified
Wed Dec 7 17:37:21 2016 UTC
(7 years, 7 months ago)
by
guez
Original Path:
trunk/Sources/phylmd/Interface_surf/interfsurf_hq.f
File length: 12681 byte(s)
Diff to
previous 208
,
to
selected 217
The program did not work with cycle_diurne set to false. mu0 in
physiq, which is supposed to be a cosine, was set to -999.999. So prmu
in swu had a value of the order of 1e3. So zrmum1 in sw2s had a value
of the order of 1e3. So zrayl in sw2s had a value of the order of
1e15. So ztray and ptauaz in swclr also had a large value. So zcorae
at line 138 in swclr had a large negative value, which resulted in
overflow at line 138.
This assignment of -999.999 to mu0 dates from somewhere between
revisions 348 and 524 of LMDZ. It was corrected in revision 1068 of
LMDZ with a call to angle which was present in revision 348. However,
procedure angle was removed from LMDZE in revision 22 because it was
not used. Hesitated to bring back angle but, finally, just removed the
option of having no diurnal cycle.
Revision
206 -
(
view)
(
annotate)
-
[select for diffs]
Modified
Tue Aug 30 12:52:46 2016 UTC
(7 years, 10 months ago)
by
guez
Original Path:
trunk/Sources/phylmd/Interface_surf/interfsurf_hq.f
File length: 12900 byte(s)
Diff to
previous 202
,
to
selected 217
Removed dimension klev of flux_[tquv] and y_flux_[tquv] in
clmain. Removed dimension klev of flux_[tquv] in physiq. Removed
dimension klev of flux_[tq] in hbtm. Removed dimension klev of
flux_[tq] in clqh and computations for layers other than the surface
layer. Removed dimension klev of flux_v in clvent and computations for
layers other than the surface layer. Values for layers other than the
surface layer were not used nor output (not even in LMDZ).
Removed argument dnwd0 of concvl. Simply write - mp in physiq
(following LMDZ).
Removed useless intermediary variables zxflux[tquv] in physiq.
Revision
202 -
(
view)
(
annotate)
-
[select for diffs]
Modified
Wed Jun 8 12:23:41 2016 UTC
(8 years, 1 month ago)
by
guez
Original Path:
trunk/Sources/phylmd/Interface_surf/interfsurf_hq.f
File length: 12979 byte(s)
Diff to
previous 191
,
to
selected 217
Promoted lmt_pas from local variable of physiq to variable of module
conf_gcm_m.
Removed variable run_off of module interface_surf. Was not
used. Called run_off_ter in LMDZ, but not used nor printed there
either.
Simplified logic in interfoce_lim. The way it was convoluted with
interfsurf_hq and clmain was quite a mess. Extracted reading of SST
into a separate procedure: read_sst. We do not need SST and pctsrf_new
at the same time: SST is not needed for sea-ice surface. I did not
like this programming: going through the procedure repeatedly for
different purposes and testing inside whether there was something to
do or it was already done. Reading is now only controlled by itap and
lmt_pas, instead of debut, jour, jour_lu and deja_lu. Now we do not
copy from pct_tmp to pctsrf_new every time step.
Simplified processing of pctsrf in clmain and below. It was quite
troubling: pctsrf_new was intent out in interfoce_lim but only defined
for ocean and sea-ice. Also the idea of having arrays for all
surfaces, pcsrf and pctsrf_new, in interfsurf_hq, which is called for
a particular surface, was troubling. pctsrf_new for all surfaces was
intent out in intefsurf_hq, but not defined for all surfaces at each
call. Removed argument pctsrf_new of clmain: was a duplicate of pctsrf
on output, and not used in physiq. Replaced pctsrf_new in clmain by
pctsrf_new_oce and pctsrf_new_sic, which were the only ones modified.
Revision
191 -
(
view)
(
annotate)
-
[select for diffs]
Modified
Mon May 9 19:56:28 2016 UTC
(8 years, 2 months ago)
by
guez
Original Path:
trunk/Sources/phylmd/Interface_surf/interfsurf_hq.f
File length: 13796 byte(s)
Diff to
previous 178
,
to
selected 217
Extracted the call to read_comdissnew out of conf_gcm.
Made ok_instan a variable of module clesphys, itau_phy a variable of
module phyetat0_m, nid_ins a variable of module ini_histins_m, itap a
variable of new module time_phylmdz, so that histwrite_phy can be
called from any procedure without the need to cascade those variables
into that procedure. Made itau_w a variable of module time_phylmdz so
that it is computed only once per time step of physics.
Extracted variables of module clesphys which were in namelist
conf_phys_nml into their own namelist, clesphys_nml, and created
procedure read_clesphys reading clesphys_nml, to avoid side effect.
No need for double precision in procedure getso4fromfile. Assume there
is a single variable for the whole year in the NetCDF file instead of
one variable per month.
Created generic procedure histwrite_phy and removed procedure
write_histins, following LMDZ. histwrite_phy has only two arguments,
can be called from anywhere, and should manage the logic of writing or
not writing into various history files with various operations. So the
test on ok_instan goes inside histwrite_phy.
Test for raz_date in phyetat0 instead of physiq to avoid side effect.
Created procedure increment_itap to avoid side effect.
Removed unnecessary differences between procedures readsulfate and
readsulfate_pi.
Revision
178 -
(
view)
(
annotate)
-
[select for diffs]
Modified
Fri Mar 11 18:47:26 2016 UTC
(8 years, 4 months ago)
by
guez
Original Path:
trunk/Sources/phylmd/Interface_surf/interfsurf_hq.f
File length: 14023 byte(s)
Diff to
previous 175
,
to
selected 217
Moved variables date0, deltat, datasz_max, ncvar_ids, point, buff_pos,
buffer, regular from module histcom_var to modules where they are
defined.
Removed procedure ioipslmpp, useless for a sequential program.
Added argument datasz_max to histwrite_real (to avoid circular
dependency with histwrite).
Removed useless variables and computations everywhere.
Changed real litteral constants from default kind to double precision
in lwb, lwu, lwvn, sw1s, swtt, swtt1, swu.
Removed unused arguments: paer of sw, sw1s, sw2s, swclr; pcldsw of
sw1s, sw2s; pdsig, prayl of swr; co2_ppm of clmain, clqh; tsol of
transp_lay; nsrf of screenp; kcrit and kknu of gwstress; pstd of
orosetup.
Added output of relative humidity.
Revision
175 -
(
view)
(
annotate)
-
[select for diffs]
Modified
Fri Feb 5 16:02:34 2016 UTC
(8 years, 5 months ago)
by
guez
Original Path:
trunk/Sources/phylmd/Interface_surf/interfsurf_hq.f
File length: 14037 byte(s)
Diff to
previous 174
,
to
selected 217
Added argument itau_phy to ini_histins, phyetat0, phytrac and
phyredem0. Removed variable itau_phy of module temps. Avoiding side
effect in etat0 and phyetat0. The procedures ini_histins, phyetat0,
phytrac and phyredem0 are all called by physiq so there is no
cascading variable penalty.
In procedure inifilr, made the condition on colat0 weaker to allow for
rounding error.
Removed arguments flux_o, flux_g and t_slab of clmain, flux_o and
flux_g of clqh and interfsurf_hq, tslab and seaice of phyetat0 and
phyredem. NetCDF variables TSLAB and SEAICE no longer in
restartphy.nc. All these variables were related to the not-implemented
slab ocean. seaice and tslab were just set to 0 in phyetat0 and never
used nor changed. flux_o and flux_g were computed in clmain but never
used in physiq.
Removed argument swnet of clqh. Was used only to compute a local
variable, swdown, which was not used.
Revision
174 -
(
view)
(
annotate)
-
[select for diffs]
Modified
Wed Nov 25 20:14:19 2015 UTC
(8 years, 7 months ago)
by
guez
Original Path:
trunk/Sources/phylmd/Interface_surf/interfsurf_hq.f
File length: 15088 byte(s)
Diff to
previous 171
,
to
selected 217
Simplifications in procedure albsno. Since veget(:, 2:) was 0,
iterations 2 to nvm of the loop computing alb_neig_grid were useless.
Useless initializations of alb_neige in procedure interfsurf_hq:
alb_neig is always computed by albsno just before being used. Useless
computation of local variable zfra in the land ice case.
Revision
171 -
(
view)
(
annotate)
-
[select for diffs]
Modified
Tue Sep 29 19:48:59 2015 UTC
(8 years, 9 months ago)
by
guez
Original Path:
trunk/Sources/phylmd/Interface_surf/interfsurf_hq.f
File length: 15272 byte(s)
Diff to
previous 156
,
to
selected 217
Removed argument ierr of abort_gcm. It was always 1 and not used.
Just encapsulated pres2lev into a module.
Removed test on run_off in procedure calcul_fluxs. Useless. The test
is always done just before in interfsurf_hq.
Removed named constants rea and repsm in module suphec: never used.
Revision
154 -
(
view)
(
annotate)
-
[select for diffs]
Modified
Tue Jul 7 17:49:23 2015 UTC
(9 years ago)
by
guez
Original Path:
trunk/Sources/phylmd/Interface_surf/interfsurf_hq.f
File length: 15374 byte(s)
Diff to
previous 150
,
to
selected 217
Removed argument dtphys of physiq. Use it directly from comconst in
physiq instead.
Donwgraded variables eignfnu, eignfnv of module inifgn_m to dummy
arguments of SUBROUTINE inifgn. They were not used elsewhere than in
the calling procedure inifilr. Renamed argument dv of inifgn to eignval_v.
Made alboc and alboc_cd independent of the size of arguments. Now we
can call them only at indices knindex in interfsurf_hq, where we need
them. Fixed a bug in alboc_cd: rmu0 was modified, and the
corresponding actual argument in interfsurf_hq is an intent(in)
argument of interfsurf_hq.
Variables of size knon instead of klon in interfsur_lim and interfsurf_hq.
Removed argument alb_new of interfsurf_hq because it was the same than
alblw. Simplified test on cycle_diurne, following LMDZ.
Moved tests on nbapp_rad from physiq to read_clesphys2. No need for
separate counter itaprad, we can use itap. Define lmt_pas and radpas
from integer input parameters instead of real-type computed values.
Revision
130 -
(
view)
(
annotate)
-
[select for diffs]
Modified
Tue Feb 24 15:43:51 2015 UTC
(9 years, 4 months ago)
by
guez
File length: 16361 byte(s)
Diff to
previous 125
,
to
selected 217
The information in argument rdayvrai of calfis was redundant with the
information in argument time. Furthermore, in the physics part of gcm,
we need separately the day number (an integer) and the time of
day. So, replaced real argument rdayvrai of calfis containing elapsed
time by integer argument dayvrai containing day number. Corresponding
change in leapfrog. In procedure physiq, replaced real argument
rdayvrai by integer argument dayvrai. In procedures readsulfate and
readsulfate_preind, replaced real argument r_day by arguments dayvrai
and time.
In procedure alboc, replaced real argument rjour by integer argument
jour. alboc was always called by interfsurf_hq with actual argument
real(jour), and the meaning of the dummy argument in alboc seems to be
that it should be an integer.
In procedure leapfrog, local variable time could not be > 1. Removed
test.
In physiq, replaced nint(rdayvrai) by dayvrai. This changes the
results since julien now changes at 0 h instead of 12 h. This follows
LMDZ, where the argument of ozonecm is days_elapsed+1.
Revision
116 -
(
view)
(
annotate)
-
[select for diffs]
Modified
Thu Dec 4 16:35:03 2014 UTC
(9 years, 7 months ago)
by
guez
File length: 16339 byte(s)
Diff to
previous 106
,
to
selected 217
In test_disvert, write output files before testing order of pressure
values, so we have more information if there is a problem.
Revision
104 -
(
view)
(
annotate)
-
[select for diffs]
Modified
Thu Sep 4 10:05:52 2014 UTC
(9 years, 10 months ago)
by
guez
File length: 16257 byte(s)
Diff to
previous 101
,
to
selected 217
Removed procedure sortvarc0. Called sortvarc with an additional
argument resetvarc instead. (Following LMDZ.) Moved current time
computations and some printing statements from sortvarc to
caldyn. Could then remove arguments itau and time_0 of sortvarc, and
could remove "use dynetat0". Better to keep "dynetat0.f" as a gcm-only
file.
Moved some variables from module ener to module sortvarc.
Split file "mathelp.f" into single-procedure files.
Removed unused argument nadv of adaptdt. Removed dimension arguments
of bernoui.
Removed unused argument nisurf of interfoce_lim. Changed the size of
argument lmt_sst of interfoce_lim from klon to knon. Removed case when
newlmt is false.
dynredem1 is called only once in each run, either ce0l or gcm. So
variable nb in call to nf95_put_var was always 1. Removed variable nb.
Removed dimension arguments of calcul_fluxs. Removed unused arguments
precip_rain, precip_snow, snow of calcul_fluxs. Changed the size of
all the arrays in calcul_fluxs from klon to knon.
Removed dimension arguments of fonte_neige. Changed the size of all
the arrays in fonte_neige from klon to knon.
Changed the size of arguments tsurf and tsurf_new of interfsurf_hq
from klon to knon. Changed the size of argument ptsrf of soil from
klon to knon.
Revision
101 -
(
view)
(
annotate)
-
[select for diffs]
Modified
Mon Jul 7 17:45:21 2014 UTC
(10 years ago)
by
guez
File length: 15716 byte(s)
Diff to
previous 100
,
to
selected 217
Removed unused files "interfoce_slab.f" and "gath2cpl.f". Removed
unused variables coastalflow and riverflow of module
interface_surf. Removed unused arguments cal, radsol, dif_grnd,
fluxlat, fluxsens, dflux_s, dflux_l of procedure fonte_neige. Removed
unused arguments tslab, seaice of procedure interfsurf_hq and
clqh. Removed unused arguments seaice of procedure clmain.
In interfsurf_hq, used variable soil_model of module clesphys2 instead
of cascading it as an argument from physiq.
In phyetat0, stop if masque not found.
Variable TS instead of "TS[0-9][0-9]" in "(re)startphy.nc", with
additional dimension nbsrf.
Revision
100 -
(
view)
(
annotate)
-
[select for diffs]
Modified
Wed Jul 2 19:07:58 2014 UTC
(10 years ago)
by
guez
File length: 21294 byte(s)
Diff to
previous 99
,
to
selected 217
Removed unused variable tmp_rlic of module interface_surf.
Removed useless call to gath2cpl in procedure interfsurf_hq. Removed
then uncalled procedure gath2cpl.
Revision
99 -
(
view)
(
annotate)
-
[select for diffs]
Modified
Wed Jul 2 18:39:15 2014 UTC
(10 years ago)
by
guez
File length: 21506 byte(s)
Diff to
previous 98
,
to
selected 217
Created procedure test_disvert (following LMDZ). Added procedures
hybrid and funcd in module disvert_m. Upgraded compute_ab from
internal procedure of disvert to module procedure. Added variables y,
ya in module disvert_m. Upgraded s from local variable of procedure
disvert to module variable.
Renamed allowed value of variable vert_sampling in procedure disvert
from "read" to "read_hybrid". Added possibility to read pressure
values, value "read_pressure". Replaced vertical distribution for
value "param" by the distribution "strato_correct" from LMDZ (but kept
the value "param"). In case "tropo", replaced 1 by dsigmin (following
LMDZ). In case "strato", replaced 0.3 by dsigmin (following LMDZ).
Changed computation of bp in procedure compute_ab.
Removed debugindex case in clmain. Removed useless argument rlon of
procedure clmain. Removed useless variables ytaux, ytauy of procedure
clmain.
Removed intermediary variables tsol, qsol, tsolsrf, tslab in procedure
etat0.
Removed variable ok_veget:. coupling with the model Orchid is not
possible. Removed variable ocean: modeling an ocean slab is not
possible.
Removed useless variables tmp_rriv and tmp_rcoa from module
interface_surf.
Moved initialization of variables da, mp, phi in procedure physiq to
to inside the test iflag_con >= 3.
Revision
98 -
(
view)
(
annotate)
-
[select for diffs]
Modified
Tue May 13 17:23:16 2014 UTC
(10 years, 2 months ago)
by
guez
File length: 24544 byte(s)
Diff to
previous 82
,
to
selected 217
Split inter_barxy.f : one procedure per module, one module per
file. Grouped the files into a directory.
Split orbite.f.
Value of raz_date read from the namelist is taken into account
(resetting the step counter) even if annee_ref == anneeref and day_ref
== dayref. raz_date is no longer modified by gcm main unit. (Following
LMDZ.)
Removed argument klon of interfsur_lim. Renamed arguments lmt_alb,
lmt_rug to alb_new, z0_new (same name as corresponding actual
arguments in interfsurf_hq).
Removed argument klon of interfsurf_hq.
Removed arguments qs and d_qs of diagetpq. Were always
zero. Downgraded arguments d_qw, d_ql of diagetpq to local variables,
they were not used in physiq. Removed all computations for solid water
in diagetpq, was just zero.
Downgraded arguments fs_bound, fq_bound of diagphy to local variables,
they were not used in physiq. Encapsulated in a test on iprt all
computations in diagphy.
Removed parameter nbtr of module dimphy. Replaced it everywhere in the
program by nqmx - 2.
Removed parameter rnpb of procedure physiq. Kept the true case in
physiq and phytrac. Could not work with false case anyway.
Removed arguments klon, llm, airephy of qcheck. Removed argument ftsol
of initrrnpb, was not used.
Revision
72 -
(
view)
(
annotate)
-
[select for diffs]
Modified
Tue Jul 23 13:00:07 2013 UTC
(10 years, 11 months ago)
by
guez
Original Path:
trunk/libf/phylmd/Interface_surf/interfsurf_hq.f90
File length: 24798 byte(s)
Diff to
previous 62
,
to
selected 217
NaN to signalling NaN in gfortran_debug.mk.
Removed unused procedures in getincom and getincom2. In procedure
conf_interface, replaced call to getincom by new namelist. Moved
procedure conf_interface into module interface_surf.
Added variables sig1 and w01 to startphy.nc and restartphy.nc, for
procedure cv_driver. Renamed (ema_)?work1 and (ema_)?work2 to sig1 and
w01 in concvl and physiq.
Deleted unused arguments of clmain, clqh and intersurf_hq, among which
(y)?sollwdown. Following LMDZ, in physiq, read sollw instead of
sollwdown from startphy.nc, write sollw instead of sollwdown to
restartphy.nc.
In procedure sw, initialized zfs[ud][pn]a[di], for runs where ok_ade
and ok_aie are false. (Following LMDZ.)
Added dimension klev to startphy.nc and restartphy.nc, and deleted
dimension horizon_vertical. Made t_ancien and q_ancien two-dimensional
NetCDF variables. Bug fix: in phyetat0, define ratqs, clwcon and
rnebcon for vertical levels >=2.
Bug fix: set mfg, p[de]n_[ud] to 0. when iflag_con >= 3. (Following LMDZ.)
Revision
62 -
(
view)
(
annotate)
-
[select for diffs]
Modified
Thu Jul 26 14:37:37 2012 UTC
(11 years, 11 months ago)
by
guez
Original Path:
trunk/libf/phylmd/Interface_surf/interfsurf_hq.f90
File length: 25989 byte(s)
Diff to
previous 54
,
to
selected 217
Changed handling of compiler in compilation system.
Removed the prefix letters "y", "p", "t" or "z" in some names of variables.
Replaced calls to NetCDF by calls to NetCDF95.
Extracted "ioget_calendar" procedures from "calendar.f90" into a
separate file.
Extracted to a separate file, "mathop2.f90", procedures that were not
part of the generic interface "mathop" in "mathop.f90".
Removed computation of "dq" in "bilan_dyn", which was not used.
In "iniadvtrac", removed schemes 20 Slopes and 30 Prather. Was not
compatible with declarations of array sizes.
In "clcdrag", "ustarhb", "vdif_kcay", "yamada4" and "coefkz", changed
the size of some arrays from "klon" to "knon".
Removed possible call to "conema3" in "physiq".
Removed unused argument "cd" in "yamada".
Revision
54 -
(
view)
(
annotate)
-
[select for diffs]
Added
Tue Dec 6 15:07:04 2011 UTC
(12 years, 7 months ago)
by
guez
Original Path:
trunk/libf/phylmd/Interface_surf/interfsurf_hq.f90
File length: 26001 byte(s)
Diff
to
selected 217
Removed Numerical Recipes procedure "ran1". Replaced calls to "ran1"
in "inidissip" by calls to intrinsic procedures.
Split file "interface_surf.f90" into a file with a module containing
only variables, "interface_surf", and single-procedure files. Gathered
files into directory "Interface_surf".
Added argument "cdivu" to "gradiv" and "gradiv2", "cdivh" to
"divgrad2" and "divgrad", and "crot" to "nxgraro2" and
"nxgrarot". "dissip" now uses variables "cdivu", "cdivh" and "crot"
from module "inidissip_m", so it can pass them to "gradiv2",
etc. Thanks to this modification, we avoid a circular dependency
betwwen "inidissip.f90" and "gradiv2.f90", etc. The value -1. used by
"gradiv2", for instance, during computation of eigenvalues is not the
value "cdivu" computed by "inidissip".
Extracted procedure "start_inter_3d" from module "startdyn", to its
own module.
In "inidissip", unrolled loop on "ii". I find it clearer now.
Moved variables "matriceun", "matriceus", "matricevn", "matricevs",
"matrinvn" and "matrinvs" from module "parafilt" to module
"inifilr_m". Moved variables "jfiltnu", "jfiltnv", "jfiltsu",
"jfiltsv" from module "coefils" to module "inifilr_m".
This form allows you to request diffs between any two revisions of this file.
For each of the two "sides" of the diff,
enter a numeric revision.