ROMS_r215: /home/cholod/wd/peps/ROMS/romsagrif_r215/Roms_tools/Roms_Agrif/compute_extended

00001 !
00002 ! $Id: compute_extended_bounds.h,v 1.2 2003/12/17 13:56:00 pmarches Exp $
00003 !
00004 /* Auxiliary module "compute_extended_bounds.h":
00005 ------------------------------------------------------
00006  Bounds designed to cover interior points of an array and
00007  ghost points associated with PHYSICAL side boundaries AND
00008  also ghost points of internal computational boundaries of
00009  MPI-subdomains.
00010 
00011  In the case of shared memory code (non-MPI version)
00012  IstrR, IendR,JstrR,JendR are equivalent to that computed
00013  by "compute_auxiliary_bounds.h".
00014 */
00015       integer IstrR,IendR,JstrR,JendR
00016       if (Istr.eq.1) then
00017 # ifdef EW_PERIODIC
00018         IstrR=Istr-2
00019 # else
00020 #  ifdef MPI
00021         if (WEST_INTER) then
00022           IstrR=Istr-2
00023         else
00024           IstrR=Istr-1
00025         endif
00026 #  else
00027         IstrR=Istr-1
00028 #  endif
00029 # endif
00030       else
00031         IstrR=Istr
00032       endif
00033 
00034 #ifdef MPI
00035       if (Iend.eq.Lmmpi) then
00036 #else
00037       if (Iend.eq.Lm) then
00038 #endif      
00039 # ifdef EW_PERIODIC
00040         IendR=Iend+2
00041 # else
00042 #  ifdef MPI
00043         if (EAST_INTER) then
00044           IendR=Iend+2
00045         else
00046           IendR=Iend+1
00047         endif
00048 #  else
00049         IendR=Iend+1
00050 #  endif
00051 # endif
00052       else
00053         IendR=Iend
00054       endif
00055 
00056       if (Jstr.eq.1) then
00057 # ifdef NS_PERIODIC
00058         JstrR=Jstr-2
00059 # else
00060 #  ifdef MPI
00061         if (SOUTH_INTER) then
00062           JstrR=Jstr-2
00063         else
00064           JstrR=Jstr-1
00065         endif
00066 #  else
00067         JstrR=Jstr-1
00068 #  endif
00069 # endif
00070       else
00071         JstrR=Jstr
00072       endif
00073 
00074 #ifdef MPI
00075       if (Jend.eq.Mmmpi) then
00076 #else
00077       if (Jend.eq.Mm) then
00078 #endif 
00079 # ifdef NS_PERIODIC
00080         JendR=Jend+2
00081 # else
00082 #  ifdef MPI
00083         if (NORTH_INTER) then
00084           JendR=Jend+2
00085         else
00086           JendR=Jend+1
00087         endif
00088 #  else
00089         JendR=Jend+1
00090 #  endif
00091 # endif
00092       else
00093         JendR=Jend
00094       endif
00095
/home/cholod/wd/peps/ROMS/romsagrif_r215/Roms_tools/Roms_Agrif/compute_extended_bounds.h
