pkg/atm_ocn_coupler/atm_to_ocn_maprunoff.F

C $Header: /u/gcmpack/MITgcm/pkg/aim_ocn_coupler/atm_to_ocn_maprunoff.F,v 1.1 2003/12/15 02:28:00 jmc Exp $
C $Name:  $

#include "CPP_OPTIONS.h"

CStartOfInterface
      SUBROUTINE ATM_TO_OCN_MAPRUNOFF( 
     I           fld_atm, 
     U           fld_ocn )
C     *==========================================================*
C     | SUBROUTINE ATM_TO_OCN_MAPRUNOFF
C     | o maps RunOff from atmos. grid (land) to ocean grid
C     |    and Incorporate to ocean FreshWater flux
C     *==========================================================*
      IMPLICIT NONE

#include "ATMSIZE.h"
#include "OCNSIZE.h"
#include "CPL_MAP2GRIDS.h"

C     == Routine arguments ==
C     fld_atm - Atmos. grid field
C     fld_ocn - Ocean grid field (atmos field will be mapped to this field)
      _RL  fld_atm(Nx_atm*Ny_atm)
      _RL  fld_ocn(Nx_ocn*Ny_ocn)
CEndOfInterface

C     == Local variables ==
      INTEGER n, ijo, ija

C       Set to zero's everywhere
c       DO ijo=1,Nx_ocn*Ny_ocn
c         fld_ocn(ijo)=0.
c       ENDDO

C       Maps Runoff from land grid point to ocean grid point 
C       and update oceanic FreshWater flux (substract RunOff)
        DO n=1,nROmap
          ijo=ijROocn(n)
          ija=ijROatm(n)
          fld_ocn(ijo)=fld_ocn(ijo)
     &                -fld_atm(ija)*arROmap(n)
        ENDDO

      RETURN
      END
1	C $Header: /u/gcmpack/MITgcm/pkg/aim_ocn_coupler/atm_to_ocn_maprunoff.F,v 1.1 2003/12/15 02:28:00 jmc Exp $
2	C $Name: $
3
4	#include "CPP_OPTIONS.h"
5
6	CStartOfInterface
7	SUBROUTINE ATM_TO_OCN_MAPRUNOFF(
8	I fld_atm,
9	U fld_ocn )
10	C ==========================================================
11	C \| SUBROUTINE ATM_TO_OCN_MAPRUNOFF
12	C \| o maps RunOff from atmos. grid (land) to ocean grid
13	C \| and Incorporate to ocean FreshWater flux
14	C ==========================================================
15	IMPLICIT NONE
16
17	#include "ATMSIZE.h"
18	#include "OCNSIZE.h"
19	#include "CPL_MAP2GRIDS.h"
20
21	C == Routine arguments ==
22	C fld_atm - Atmos. grid field
23	C fld_ocn - Ocean grid field (atmos field will be mapped to this field)
24	_RL fld_atm(Nx_atm*Ny_atm)
25	_RL fld_ocn(Nx_ocn*Ny_ocn)
26	CEndOfInterface
27
28	C == Local variables ==
29	INTEGER n, ijo, ija
30
31	C Set to zero's everywhere
32	c DO ijo=1,Nx_ocn*Ny_ocn
33	c fld_ocn(ijo)=0.
34	c ENDDO
35
36	C Maps Runoff from land grid point to ocean grid point
37	C and update oceanic FreshWater flux (substract RunOff)
38	DO n=1,nROmap
39	ijo=ijROocn(n)
40	ija=ijROatm(n)
41	fld_ocn(ijo)=fld_ocn(ijo)
42	& -fld_atm(ija)*arROmap(n)
43	ENDDO
44
45	RETURN
46	END