Changeset 12377 for NEMO/trunk/src/TOP/C14/trcatm_c14.F90
- Timestamp:
- 2020-02-12T15:39:06+01:00 (4 years ago)
- Location:
- NEMO/trunk
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/trunk
- Property svn:externals
-
old new 3 3 ^/utils/build/mk@HEAD mk 4 4 ^/utils/tools@HEAD tools 5 ^/vendors/AGRIF/dev @HEAD ext/AGRIF5 ^/vendors/AGRIF/dev_r11615_ENHANCE-04_namelists_as_internalfiles_agrif@HEAD ext/AGRIF 6 6 ^/vendors/FCM@HEAD ext/FCM 7 7 ^/vendors/IOIPSL@HEAD ext/IOIPSL
-
- Property svn:externals
-
NEMO/trunk/src/TOP/C14/trcatm_c14.F90
r10069 r12377 21 21 PUBLIC trc_atm_c14_ini ! called in trcini_c14.F90 22 22 ! 23 !! * Substitutions 24 # include "do_loop_substitute.h90" 23 25 !!---------------------------------------------------------------------- 24 26 !! NEMO/TOP 4.0 , NEMO Consortium (2018) … … 118 120 IF( ierr3 /= 0 ) CALL ctl_stop( 'STOP', 'trc_atm_c14_ini: unable to allocate fareaz' ) 119 121 ! 120 DO jj = 1 , jpj ! from C14b package 121 DO ji = 1 , jpi 122 IF( gphit(ji,jj) >= yn40 ) THEN 123 fareaz(ji,jj,1) = 0. 124 fareaz(ji,jj,2) = 0. 125 fareaz(ji,jj,3) = 1. 126 ELSE IF( gphit(ji,jj ) <= ys40) THEN 127 fareaz(ji,jj,1) = 1. 128 fareaz(ji,jj,2) = 0. 129 fareaz(ji,jj,3) = 0. 130 ELSE IF( gphit(ji,jj) >= yn20 ) THEN 131 fareaz(ji,jj,1) = 0. 132 fareaz(ji,jj,2) = 2. * ( 1. - gphit(ji,jj) / yn40 ) 133 fareaz(ji,jj,3) = 2. * gphit(ji,jj) / yn40 - 1. 134 ELSE IF( gphit(ji,jj) <= ys20 ) THEN 135 fareaz(ji,jj,1) = 2. * gphit(ji,jj) / ys40 - 1. 136 fareaz(ji,jj,2) = 2. * ( 1. - gphit(ji,jj) / ys40 ) 137 fareaz(ji,jj,3) = 0. 138 ELSE 139 fareaz(ji,jj,1) = 0. 140 fareaz(ji,jj,2) = 1. 141 fareaz(ji,jj,3) = 0. 142 ENDIF 143 END DO 144 END DO 122 DO_2D_11_11 123 IF( gphit(ji,jj) >= yn40 ) THEN 124 fareaz(ji,jj,1) = 0. 125 fareaz(ji,jj,2) = 0. 126 fareaz(ji,jj,3) = 1. 127 ELSE IF( gphit(ji,jj ) <= ys40) THEN 128 fareaz(ji,jj,1) = 1. 129 fareaz(ji,jj,2) = 0. 130 fareaz(ji,jj,3) = 0. 131 ELSE IF( gphit(ji,jj) >= yn20 ) THEN 132 fareaz(ji,jj,1) = 0. 133 fareaz(ji,jj,2) = 2. * ( 1. - gphit(ji,jj) / yn40 ) 134 fareaz(ji,jj,3) = 2. * gphit(ji,jj) / yn40 - 1. 135 ELSE IF( gphit(ji,jj) <= ys20 ) THEN 136 fareaz(ji,jj,1) = 2. * gphit(ji,jj) / ys40 - 1. 137 fareaz(ji,jj,2) = 2. * ( 1. - gphit(ji,jj) / ys40 ) 138 fareaz(ji,jj,3) = 0. 139 ELSE 140 fareaz(ji,jj,1) = 0. 141 fareaz(ji,jj,2) = 1. 142 fareaz(ji,jj,3) = 0. 143 ENDIF 144 END_2D 145 145 ! 146 146 ENDIF
Note: See TracChangeset
for help on using the changeset viewer.