Changes between oasis3_3 and oasis3_prism_2_5

The changes between version oasis3_3 and version oasis3_prism_2_5 delivered in September 2006 are the following:

• Bug corrections:
  • In oasis3/lib/scrip/src/remap_bilinear.f, remap_bicubic.f,
    remap_bilinear_reduced.f, remap_bicubic_reduced.F90: (r2084) we observed a wrong behaviour of routines remap_bilinear.f and remap_bicubic.f on the NEC SX9 when compiled with NEC SX compiler revision 400. We got round this problem by adding an explicit instruction to prevent the vectorisation of one loop.

    As remap_bilinear_reduced.f and remap_bicubic_reduced.F90 have a very similar loop, we introduced the same instruction in these routines, although nothing specific was observed with them.

  • oasis3/lib/scrip/src/remap_bicubic_reduced.F90: (r1883) Two bugfixes: 1) Memory fault when a target point was falling on the before-last latitude circle or in the before-last latitude band of the source reduced grid; 2) Wrong neighbours for target points south of the before-last latitude circle (i.e in the last latitude band or Southern).

  • oasis3/lib/psmile/src/mod_psmile_io.F90: (r2380) correction to ensure that when INVERT is used, the corner latitudes and longitudes are also inverted (and not only the center latitudes and longitudes as before).

  • oasis3/lib/scrip/src/scriprmp.F: (r1547) the calculation of the average for the line of points at the pole was bugged and did not have any effect. It is now debugged but commented.

  • In oasis3/lib/scrip/src/vector.F90:

    - correction of wrong sequences in declarations, at least for Intel & NAG compilers (thanks to L. Kornblueh from MPI); bug fix announced to the mailing list diff-oasis@cerfacs.fr on 31/10/2006.

    - (r1698 - 2008-08-20) bugfix to make sure that OASIS does not automatically calculates corners of target grid as this calculation is correct only for LR grids and target grid type is not known (thanks to S. Calmanti from Météo-France)

    - Modifications so that last 4 arguments of call to grid_init are always arrays even when corners are not defined (error detected with Intel Fortran V10.1.012 by Mike Rezny, SGI, Australia)

  • In oasis3/lib/clim.GSIP/src/CLIM_Init_Oasis.F, correction of a wrongly positioned #endif (thanks to L. Kornblueh from MPI); bug fix announced to the mailing list diff-oasis@cerfacs.fr on 31/10/2006.

  • In oasis3/lib/psmile/src/prism_enddef_proto.F, the call to MPI_Errhandler_set was moved after the test on the value of mpi_err returned by MPI_Buffer_Detach (thanks to I. Bethke from NERSC); bug fix announced to the mailing list diff-oasis@cerfacs.fr on 14/12/2006.

  • Routine oasis3/lib/psmile/src/prism_terminate_proto.F90 was modified to ensure proper deallocation of all allocated arrays (thanks to Adam Ralph from ICHEC)

  • In oasis3/lib/scrip/remap_conserv.F, small bugfix having no impact on the results, it just avoids misleading messages of type "Error sum wts map1:grid2_add ..." to be printed in the cplout log file.

  • In oasis3/src/getfld.F, givfld.F, driver.f and closerst.F, correction of a bug observed by A.Caubel from CEA for coupling fields having a sequence index greater than 1. For first iteration, closing of netcdf restart files is done in driver.f by calling new routine closerst.F. Closing is therefore removed from getfld.F. A minor correction (useless opening and closing of first netcdf restart file) was also added to givfld.F. Bug fix announced to the mailing list diff-oasis@cerfacs.fr on 09/02/2006.

  • In oasis3/src/inipar.F, bugfix for SEQMODE greater than 9 (thanks to T. Silva, Oregon State U.) and to avoid array overbound.

  • In oasis3/util/make_dir/TopMakefileOasis3: typo error: libmpp_io.a instead of libmppio.a (thanks to J.M. Epitalon from CERFACS)

  • In oasis3/lib/psmile/src/prism_init_comp_proto.F: initialisation of iprcou (bug fix thanks to J.M. Epitalon).

  • In oasis3/src/filling.f: rewind of file in the "AN" case (thanks to S. Calmanti from Météo-France)

  • In oasis3/lib/psmile/src/mod_prism_put_proto.F90: bug fix to avoid array overbounds (bug identified T. van Noije from KNMI).

  • In oasis3/lib/psmile/src/mod_psmile_io.F90: thanks to A. Caubel from CEA, modification so that the grid point longitudes and latitudes do not have to appear before the corner longitudes and latitudes in the grids.nc file (revision 13/01/2009).

  • Bugfix in oasis3/examples/testNONE/PROG/calc_errorfield.f90 to use the absolute value of the error in the non masked point mean error calculation.

  • Modifications in oasis3/lib/clim/src/CLIM_Init_Oasis.F and oasis3/lib/
    psmile/src/prism_init_comp_proto.F to allow communication log file (*.prt* files) for OASIS and/or component models to run on up to 9999 processes.

• Other major modifications
  • New directory structure. See section 3 for details.

  • Modification of many routines to allow using more than one OASIS3 executable in a coupled model resulting in pseudo-parallelisation of OASIS3 on a field-per-field basis. See section 8.2 for more detail.

  • New directory structure and update of compiling environement to work with the new directory structure. In particular, the location of directory created for compilation (see ARCHDIR in the Makefile headers make.xxx in oasis3/util/make_dir) can be arbitrarily chosen by the user.

  • Modification of routine oasis3/lib/scrip/src/remap_distwgt.F so that
    SCRIPR/DISTWGT that has by default the same behaviour than SCRIPR/BILINEAR, /BICUBIC and /CONSERV (with FRACNNEI option) i.e. the non-masked nearest neighbour is used for target grid points having their N nearest neighbour all masked. To reproduce the previous default behaviour, one has to compile with CPP key NOT_NNEIGHBOUR. See section 6.4 for details.

  • Inclusion of an additional number below the $NFIELDS keyword in the namcouple (after the total number of fields exchanged, on the same line). This number, corresponding to the maximum number of prism_def_var_proto called by ANY component model in the coupled system, is needed only if it is greater than twice the number of fields listed in the namcouple; this may be the case if OASIS3 is used in pseudo-parallel mode or if fields declared with prism_def_var_proto call (and corresponding to sending - prism_put_proto call- or receiving - prism_get_proto call - actions in the component models) do not appear in the namcouple (in this case, the sending and receiving calls simply return without any action performed).

  • added options GLBPOS, BASPOS, BASBAL for the cooking stage transformation CONSERV. For details, see section 6.5.

  • New CPP key TREAT_OVERLAY : to ensure that if two cells of the source grid overlay, at least the one with the greater numerical index is masked (they also can be both masked). For example, if the grid line with i=1 overlaps the grid line with i=imax, it is the latter that must be masked. When this is not the case with the mask defined in masks.nc, this CPP key ensures that these rules are respected. This is mandatory for SCRIPR/CONSERV remapping, see section 6.4.

  • Optimisation of SCRIPR (XXX to be detailed) interpolation weigths-and-address files, thanks to CMCC. The weights and addresses are now read once per run and stored.

  • mode stats consommation Eric M (XXX to be detailed)

  • mod_prism_get_comm (XXX to be detailed)

  • Release of a new toy coupled model TOYOASIS3. This new toy model is described in 8.3. The toy model sources are available in oasis3/src/mod/oceoa3/, atmoa3, cheoa3. It running environment is available in oasis3/src/mod/oasis3/examples/toyoasis3. The grids of the TOYOASIS3 component models and the interpolation performed have been updated compared to the previous toy coupled model TOYCLIM (which is no longer distributed with the official release).

  • In oasis3/src/mod/oasis3/src/extrap.F: optimisation (thanks to T. Schoenemeyer from NEC) and idoitold changed from 1000000 to 10000000 to support higher resolution grids (thanks to E.Maisonnave from CERFACS). XXX and CMCC

  • In oasis3/lib/psmile/src/mod_prism_grids_writing.F90:

    - routine prism_write_angle was added to to allow a component model to write the angle of its grid (see section 4.2 for more detail).

    - changed logical netcdf for l_netcdf (to run on NEC SX9, thanks to E. Maisonnave, CERFACS)

  • In may 2007, we moved from CVS to SVN for source management.

• Other modifications
  • The names of the log files for the communication information *.prt* are now always ending with 4 digits indicating the rank of the component process (e.g. model1.prt0002 or model1.prt9999) or the rank of the oasis process (e.g. Oasis.prt0001 or Oasis.prt0010). These modifications impacted routines oasis3/lib/clim/src/CLIM_Init_Oasis.F and
    oasis3/lib/psmile/src/prism_init_comp_proto.F.

  • The following routines were modified for the NAG compiler: oasis3
    /src/getfld.F, inipar.F, interp.F and oasis3/lib/scrip/src/
    remap_write.F

  • Some routines in oasis3/lib/mpp_io/src/ were modified so that all debug and log messages are now written to stdout unit and to include modifications done in the OASIS4 version for CRAY pointers and for bundles.

  • Added CPP key __SILENT CPP key to reduce log outputs to .prt files during the psmile library exchanges (thanks to S. Lorenz from MPI)

  • In oasis3/src/chkfld.f, modified misleading comment about masked points written to cplout log file (thanks to T. Craig from BOM).

  • In oasis3/lib/psmile/src/mpp_io/src/mpp_io_mod_oa.F90 modified "lowercase" function so to avoid using "transfer" function (thanks to Mike Rezny for SGI Melbourne) which causes problem with pgf90 5.2.4, 6.1.3 or 7.0 in 64 bit mode, and with gfortran 4.2.1 and 4.2.5 SUSE Linux.