!-------------------------------------- LICENCE BEGIN ------------------------------------ !Environment Canada - Atmospheric Science and Technology License/Disclaimer, ! version 3; Last Modified: May 7, 2008. !This is free but copyrighted software; you can use/redistribute/modify it under the terms !of the Environment Canada - Atmospheric Science and Technology License/Disclaimer !version 3 or (at your option) any later version that should be found at: !http://collaboration.cmc.ec.gc.ca/science/rpn.comm/license.html ! !This software is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; !without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. !See the above mentioned License/Disclaimer for more details. !You should have received a copy of the License/Disclaimer along with this software; !if not, you can write to: EC-RPN COMM Group, 2121 TransCanada, suite 500, Dorval (Quebec), !CANADA, H9P 1J3; or send e-mail to service.rpn@ec.gc.ca !-------------------------------------- LICENCE END -------------------------------------- ***s/r set_cori - allocate and compute coriolis factor * #include "model_macros_f.h"*
subroutine set_cori 1,1 #include "impnone.cdk"
* *author * michel roch - rpn - aug 95 * *revision * v2_00 - Desgagne/Lee - initial MPI version (from setcori v1_03) * v3_00 - Desgagne & Lee - Lam configuration * *object * See above id. * *arguments * none * *implicits #include "glb_ld.cdk"
#include "lun.cdk"
#include "grd.cdk"
#include "geomg.cdk"
#include "cori.cdk"
* ** * integer dim,dyn_alloc_error real*8 ONE parameter ( ONE = 1.0 ) * * --------------------------------------------------------------- * if (Lun_out.gt.0) write(lun_out,1000) * * *C 1. allocate memory for coriolis factor on U and V grids * _____________________________________________________ * dim = LARRAY2D*2 DYNAMIC_ALLOC(Cori_fcoru_8,dim) DYNAMIC_ALLOC(Cori_fcorv_8,dim) * * * 2. compute coriolis factors on U and V grids * _________________________________________ * call coriol8
( cori_fcoru_8, cori_fcorv_8, geomg_x_8, geomg_y_8, $ geomg_xu_8, geomg_yv_8, ONE, Grd_rot_8, $ LDIST_DIM) * 1000 format( %/,'ALLOCATE AND COMPUTE CORIOLIS FACTOR (S/R SET_CORI)', %/,'==================================================') * * --------------------------------------------------------------- * return end