- Timestamp:
- 2019-02-25T15:19:53+01:00 (5 years ago)
- Location:
- branches/UKMO/dev_r5518_DMP_TOOLS/NEMOGCM/TOOLS/DMP_TOOLS/src
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/dev_r5518_DMP_TOOLS/NEMOGCM/TOOLS/DMP_TOOLS/src/coast_dist.F90
r10199 r10722 9 9 CONTAINS 10 10 11 SUBROUTINE coast_dist_weight( presto, ln_read , klev)11 SUBROUTINE coast_dist_weight( presto, zdct, ln_read , klev, ln_distcoast_calc ) 12 12 !!---------------------------------------------------------------------- 13 13 !! *** ROUTINE coast_dist_weight *** … … 20 20 IMPLICIT NONE 21 21 REAL(wp), DIMENSION(jpi,jpj), INTENT( inout ) :: presto 22 REAL(wp), DIMENSION(jpi,jpj), INTENT( inout ) :: zdct 22 23 LOGICAL, INTENT( in ) :: ln_read 23 24 INTEGER, INTENT( in ) :: klev 25 LOGICAL, INTENT( inout ) :: ln_distcoast_calc 24 26 REAL(wp), DIMENSION(jpi,jpj) :: zdct 25 27 REAL(wp) :: zinfl = 1000.e3_wp ! Distance of influence of coast line (could be … … 37 39 ELSE 38 40 39 IF ( klev == 1 ) THEN ! Since only athe surface distance to coast is used, only calculate once.41 IF (ln_distcoast_calc ) THEN ! Since only the surface distance to coast is used, only calculate once. 40 42 CALL cofdis( zdct ) 43 ln_distcoast_calc=.false. 41 44 ENDIF 42 45 -
branches/UKMO/dev_r5518_DMP_TOOLS/NEMOGCM/TOOLS/DMP_TOOLS/src/make_dmp_file.F90
r10199 r10722 32 32 CHARACTER(LEN=200) :: outfile = 'resto.nc' ! output file 33 33 REAL(wp) :: zlat, zlat2, zlat0 34 LOGICAL :: ln_distcoast_calc ! logical for calculating the distance 35 ! to coast on the first depth. 36 REAL(wp), DIMENSION(:,:), ALLOCATABLE :: zdct! distance to coast field 34 37 35 38 ! Read namelist … … 48 51 WRITE(numout, *) 'jpk = ',jpk 49 52 53 ln_distcoast_calc=.true. 54 50 55 ALLOCATE( resto(jpi, jpj) ) 56 ALLOCATE( zdct(jpi, jpj) ) 51 57 52 58 !Create output file … … 94 100 IF (ln_coast) THEN 95 101 ! Reduce damping in vicinity of coastlines 96 CALL coast_dist_weight(resto, ln_readdistcoast, jk)102 CALL coast_dist_weight(resto, zdct, ln_readdistcoast, jk, ln_distcoast_calc) 97 103 ENDIF 98 104 ENDIF
Note: See TracChangeset
for help on using the changeset viewer.