Changes between oasis3_prism_2_5 and oasis3_prism_2_4

The changes between version oasis3_prism_2_5 and version oasis3_prism_2_4 delivered in December 2004 are listed here after. Please note that those modifications should not bring any difference in the interpolation results, except for SCRIPR/DISTWGT (see below).

• Bug corrections:
  • In prism/src/lib/scrip/src/scriprmp.F: initialisation of dst_array(:); bug fix announced to the mailing list diff-oasis@cerfacs.fr on 02/02/2006.

  • In prism/src/lib/psmile/src/prism_enddef_proto.F and prism/src/lib/
    clim/src/CLIM_Start_MPI.F: the call to MPI_barrier (that created a deadlock when not all processes of a component model were exchanging data with the coupler) was changed for a call to MPI_wait on the previous MPI_Isend; bug fix announced to the mailing list diff-oasis@cerfacs.fr on 02/23/2006.

  • For SCRIPR/DISTWGT, in prism/src/lib/scrip/src/remap_distwgt.f: line 190 was repeated without epsilon modification; bug fix announced to the mailing list diff-oasis@cerfacs.fr on 03/21/2006.

  • In prism/src/lib/psmile/src/mod_prism_put_proto.F90, for prism_put
    _proto_r28 and prism_put_proto_r24, the reshape of the 2d field was moved after the test checking if the field is defined in the namcouple (thanks to Arnaud Caubel from LSCE).

• Modification in SCRIP interpolations
  • For SCRIPR interpolations (see section 6.4), the value 1.0E+20 is assigned to target grid points for which no value has been calculated if prism/src/lib/scrip/src/scriprmp.f or vector.F90 (for vector interpolation) are compiled with ll_weightot = .true..

  • For SCRIPR/GAUSWGT: if routine prism/src/lib/scrip/src/remap_gauswgt.f is compiled with ll_nnei=.true., the non-masked nearest neighbour is used for target point if all original neighbours are masked (see section 6.4).

  • For SCRIPR/BICUBIC (routine prism/src/lib/scrip/src/remap_bicubic.f), the convergence criteria was modified so to ensure convergence even in single precision.

  • For SCRIPR/CONSERV (routine prism/src/lib/scrip/src/remap_conserv.f), a test was added for non-convex cell so that integration does not stall.

  • The routine prism/src/lib/scrip/src/corners.F was modified so to abort if it is called for the target grid, as the automatic calculation of corners works only for Logically-Rectangular (LR) grids and as the target grid type is unknown. If needed, the reverse remapping, in which the current target grid become the source grid, can be done .

• Other important modifications
  • A new PSMILe routine prism/src/lib/psmile/src/prism_get_freq.F was added; this routine can be used to retrieve the coupling period of field (see section 4.6.3).

  • The routines of the mpp_io library in prism/src/lib/mpp_io changed name and were merged with the OASIS4 mpp_io library.

  • Routine prism/src/mod/oasis3/src/extrap.F was modified to ensure that the extrapolation works even if the MASK value is very big (thanks to J.M. Epitalon).

  • In the namcouple, there is no need anymore to define a lag greater than 0 (e.g. LAG=+1) for fields in mode NONE.

  • Diverse modifications were included for successful compilation with NAGW compiler: non portable use of ``kind'', etc. (thanks to Luis Kornblueh from MPI).

  • In prism/src/lib/psmile/mod_prism_get_proto.F90, a potential deadlock was removed (the master process was sending a message to itself)(thanks to Luis Kornblueh from MPI).

  • Routine prism/src/lib/scrip/src/scriprmp_vector.F90 was completely rewritten for more clarity.

  • Obsolete transformations INVERT and REVERSE were removed from the toy coupled model TOYCLIM (in file prism/util/running/toyclim/input/namcouple. This change does not affect the statistics printed in the cplout but changes the orientation of some fields in the NetCDF ouput files (see the results in prism/data/toyclim/outdata).

• Other minor modifications:
  • In prism/src/lib/psmile/src/prism_enddef_proto.F, allocation is done only for rg_field_trans or dg_field_trans depending on precision for REAL (but not for both, to save memory).

  • In few routines in prism/src/lib/clim and in prism/src/mod/oasis3, parentheses were added to make sure that && has priority over $\vert$$\vert$ in CPP instructions (thanks to A. Caubel from LSCE).

  • Routines scrip/src/corners.f, netcdf.f, and scriprmp.f were renamed
    corners.F, netcdf.F, scriprmp.F and the line ``INCLUDE 'netcdf.inc' '' was changed for ``#include $<$netcdf.inc$>$ ''