Changeset 7192
- Timestamp:
- 2016-11-04T12:23:02+01:00 (8 years ago)
- Location:
- branches/2016/dev_r7012_ROBUST5_CNRS/NEMOGCM/NEMO/TOP_SRC
- Files:
-
- 8 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2016/dev_r7012_ROBUST5_CNRS/NEMOGCM/NEMO/TOP_SRC/C14/trcatm_c14.F90
r7041 r7192 160 160 DO jn = 1, nrecc14 161 161 READ(inum2,*) iyear,incom,incom,atmc14(jn) 162 tyrc14(jn)=1950._wp- float(iyear) ! BP to AD dates162 tyrc14(jn)=1950._wp-REAL(iyear,wp) ! BP to AD dates 163 163 END DO 164 164 CLOSE(inum2) -
branches/2016/dev_r7012_ROBUST5_CNRS/NEMOGCM/NEMO/TOP_SRC/CFC/trcsms_cfc.F90
r7041 r7192 113 113 ! time interpolation at time kt 114 114 DO jm = 1, jphem 115 zpatm(jm,jl) = ( p_cfc(iyear_beg, jm, jl) * FLOAT (im1) &116 & + p_cfc(iyear_end, jm, jl) * FLOAT (im2) ) / 12.115 zpatm(jm,jl) = ( p_cfc(iyear_beg, jm, jl) * REAL(im1, wp) & 116 & + p_cfc(iyear_end, jm, jl) * REAL(im2, wp) ) / 12. 117 117 END DO 118 118 -
branches/2016/dev_r7012_ROBUST5_CNRS/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zfechem.F90
r7177 r7192 136 136 DO jj = 1, jpj 137 137 DO ji = 1, jpi 138 zlight = etot(ji,jj,jk) * zstrn(ji,jj) * float(2-jn)139 zzstrn2 = zstrn2(ji,jj) * float(2-jn) + (1. - zstrn2(ji,jj) ) * float(jn-1)138 zlight = etot(ji,jj,jk) * zstrn(ji,jj) * REAL( 2-jn, wp ) 139 zzstrn2 = zstrn2(ji,jj) * REAL( 2-jn, wp ) + (1. - zstrn2(ji,jj) ) * REAL( jn-1, wp ) 140 140 ! Calculate ligand concentrations : assume 2/3rd of excess goes to 141 141 ! strong ligands (L1) and 1/3rd to weak ligands (L2) … … 192 192 zphi = ACOS( zfff ) 193 193 DO jic = 1, 3 194 zfunc = -2 * zr * COS( zphi / 3. + 2. * FLOAT( jic - 1) * rpi / 3. ) - za2 / 3.194 zfunc = -2 * zr * COS( zphi / 3. + 2. * REAL( jic - 1, wp ) * rpi / 3. ) - za2 / 3. 195 195 IF( zfunc > 0. .AND. zfunc <= ztfe) zxs = zfunc 196 196 END DO -
branches/2016/dev_r7012_ROBUST5_CNRS/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zpoc.F90
r7162 r7192 529 529 IF (jcpoc > 1) THEN 530 530 ! 531 remindelta = log(4. * 1000. ) / float(jcpoc-1)532 reminup = 1./ 400. * exp(remindelta)531 remindelta = LOG(4. * 1000. ) / REAL(jcpoc-1, wp) 532 reminup = 1./ 400. * EXP(remindelta) 533 533 ! 534 534 ! Discretization based on incomplete gamma functions … … 540 540 reminp(1) = gamain(reminup, rshape+1.0, ifault) * xremip / alphan(1) 541 541 DO jn = 2, jcpoc-1 542 reminup = 1./ 400. * exp(float(jn) * remindelta)543 remindown = 1. / 400. * exp(float(jn-1) * remindelta)542 reminup = 1./ 400. * EXP( REAL(jn, wp) * remindelta) 543 remindown = 1. / 400. * EXP( REAL(jn-1, wp) * remindelta) 544 544 alphan(jn) = gamain(reminup, rshape, ifault) - gamain(remindown, rshape, ifault) 545 545 reminp(jn) = gamain(reminup, rshape+1.0, ifault) - gamain(remindown, rshape+1.0, ifault) 546 546 reminp(jn) = reminp(jn) * xremip / alphan(jn) 547 547 END DO 548 remindown = 1. / 400. * exp(float(jcpoc-1) * remindelta)548 remindown = 1. / 400. * EXP( REAL(jcpoc-1, wp) * remindelta) 549 549 alphan(jcpoc) = 1.0 - gamain(remindown, rshape, ifault) 550 550 reminp(jcpoc) = 1.0 - gamain(remindown, rshape+1.0, ifault) -
branches/2016/dev_r7012_ROBUST5_CNRS/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zsbc.F90
r7175 r7192 321 321 END DO 322 322 CALL iom_close( numdust ) 323 ztimes_dust = 1._wp / FLOAT( ntimes_dust)323 ztimes_dust = 1._wp / REAL(ntimes_dust, wp) 324 324 sumdepsi = 0.e0 325 325 DO jm = 1, ntimes_dust … … 386 386 END DO 387 387 CALL iom_close( numriv ) 388 ztimes_riv = 1._wp / FLOAT(ntimes_riv)388 ztimes_riv = 1._wp / REAL(ntimes_riv, wp) 389 389 DO jm = 1, ntimes_riv 390 390 rivinput(ifpr) = rivinput(ifpr) + glob_sum( zriver(:,:,jm) * tmask(:,:,1) * ztimes_riv ) … … 433 433 END DO 434 434 CALL iom_close( numdepo ) 435 ztimes_ndep = 1._wp / FLOAT( ntimes_ndep)435 ztimes_ndep = 1._wp / REAL(ntimes_ndep, wp) 436 436 nitdepinput = 0._wp 437 437 DO jm = 1, ntimes_ndep -
branches/2016/dev_r7012_ROBUST5_CNRS/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zsink.F90
r7162 r7192 136 136 IF( tmask(ji,jj,jk) == 1 ) THEN 137 137 zwsmax = 0.5 * e3t_n(ji,jj,jk) / xstep 138 wsbio3(ji,jj,jk) = MIN( wsbio3(ji,jj,jk), zwsmax * FLOAT( iiter1) )139 wsbio4(ji,jj,jk) = MIN( wsbio4(ji,jj,jk), zwsmax * FLOAT( iiter2) )138 wsbio3(ji,jj,jk) = MIN( wsbio3(ji,jj,jk), zwsmax * REAL( iiter1, wp ) ) 139 wsbio4(ji,jj,jk) = MIN( wsbio4(ji,jj,jk), zwsmax * REAL( iiter2, wp ) ) 140 140 ENDIF 141 141 END DO … … 194 194 IF( tmask(ji,jj,jk) == 1 ) THEN 195 195 zwsmax = 0.5 * e3t_n(ji,jj,jk) / xstep 196 wsfep(ji,jj,jk) = MIN( wsfep(ji,jj,jk), zwsmax * FLOAT( iiter1) )196 wsfep(ji,jj,jk) = MIN( wsfep(ji,jj,jk), zwsmax * REAL( iiter1, wp ) ) 197 197 ENDIF 198 198 END DO … … 310 310 CALL wrk_alloc( jpi, jpj, jpk, ztraz, zakz, zwsink2, ztrb ) 311 311 312 zstep = rfact2 / FLOAT( kiter) / 2.312 zstep = rfact2 / REAL( kiter, wp ) / 2. 313 313 314 314 ztraz(:,:,:) = 0.e0 -
branches/2016/dev_r7012_ROBUST5_CNRS/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zsms.F90
r7162 r7192 88 88 IF( ( ln_top_euler .AND. kt == nittrc000 ) .OR. ( .NOT.ln_top_euler .AND. kt <= nittrc000 + nn_dttrc ) ) THEN 89 89 rfactr = 1. / rfact 90 rfact2 = rfact / FLOAT( nrdttrc)90 rfact2 = rfact / REAL( nrdttrc, wp ) 91 91 rfact2r = 1. / rfact2 92 92 xstep = rfact2 / rday ! Time step duration for biology -
branches/2016/dev_r7012_ROBUST5_CNRS/NEMOGCM/NEMO/TOP_SRC/TRP/trcrad.F90
r7068 r7192 214 214 IF( l_trdtrc ) THEN 215 215 ! 216 zs2rdt = 1. / ( 2. * rdt * FLOAT( nn_dttrc) )216 zs2rdt = 1. / ( 2. * rdt * REAL( nn_dttrc, wp ) ) 217 217 ztrtrdb(:,:,:) = ( ptrb(:,:,:,jn) - ztrtrdb(:,:,:) ) * zs2rdt 218 218 ztrtrdn(:,:,:) = ( ptrn(:,:,:,jn) - ztrtrdn(:,:,:) ) * zs2rdt
Note: See TracChangeset
for help on using the changeset viewer.