Changeset 7192


Ignore:
Timestamp:
2016-11-04T12:23:02+01:00 (4 years ago)
Author:
cetlod
Message:

new top interface : minor bugfixes

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  
    160160            DO jn = 1, nrecc14 
    161161               READ(inum2,*) iyear,incom,incom,atmc14(jn) 
    162                tyrc14(jn)=1950._wp-float(iyear)                ! BP to AD dates 
     162               tyrc14(jn)=1950._wp-REAL(iyear,wp)                ! BP to AD dates 
    163163            END DO 
    164164            CLOSE(inum2) 
  • branches/2016/dev_r7012_ROBUST5_CNRS/NEMOGCM/NEMO/TOP_SRC/CFC/trcsms_cfc.F90

    r7041 r7192  
    113113         ! time interpolation at time kt 
    114114         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. 
    117117         END DO 
    118118          
  • branches/2016/dev_r7012_ROBUST5_CNRS/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zfechem.F90

    r7177 r7192  
    136136            DO jj = 1, jpj 
    137137               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 ) 
    140140                  ! Calculate ligand concentrations : assume 2/3rd of excess goes to 
    141141                  ! strong ligands (L1) and 1/3rd to weak ligands (L2) 
     
    192192                        zphi = ACOS( zfff ) 
    193193                        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. 
    195195                           IF( zfunc > 0. .AND. zfunc <= ztfe)  zxs = zfunc 
    196196                        END DO 
  • branches/2016/dev_r7012_ROBUST5_CNRS/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zpoc.F90

    r7162 r7192  
    529529      IF (jcpoc > 1) THEN 
    530530         ! 
    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) 
    533533         ! 
    534534         ! Discretization based on incomplete gamma functions 
     
    540540         reminp(1) = gamain(reminup, rshape+1.0, ifault) * xremip / alphan(1) 
    541541         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) 
    544544            alphan(jn) = gamain(reminup, rshape, ifault) - gamain(remindown, rshape, ifault) 
    545545            reminp(jn) = gamain(reminup, rshape+1.0, ifault) - gamain(remindown, rshape+1.0, ifault) 
    546546            reminp(jn) = reminp(jn) * xremip / alphan(jn) 
    547547         END DO 
    548          remindown = 1. / 400. * exp(float(jcpoc-1) * remindelta) 
     548         remindown = 1. / 400. * EXP( REAL(jcpoc-1, wp) * remindelta) 
    549549         alphan(jcpoc) = 1.0 - gamain(remindown, rshape, ifault) 
    550550         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  
    321321            END DO 
    322322            CALL iom_close( numdust ) 
    323             ztimes_dust = 1._wp / FLOAT( ntimes_dust )  
     323            ztimes_dust = 1._wp / REAL(ntimes_dust, wp)  
    324324            sumdepsi = 0.e0 
    325325            DO jm = 1, ntimes_dust 
     
    386386               END DO 
    387387               CALL iom_close( numriv ) 
    388                ztimes_riv = 1._wp / FLOAT(ntimes_riv)  
     388               ztimes_riv = 1._wp / REAL(ntimes_riv, wp)  
    389389               DO jm = 1, ntimes_riv 
    390390                  rivinput(ifpr) = rivinput(ifpr) + glob_sum( zriver(:,:,jm) * tmask(:,:,1) * ztimes_riv )  
     
    433433            END DO 
    434434            CALL iom_close( numdepo ) 
    435             ztimes_ndep = 1._wp / FLOAT( ntimes_ndep )  
     435            ztimes_ndep = 1._wp / REAL(ntimes_ndep, wp)  
    436436            nitdepinput = 0._wp 
    437437            DO jm = 1, ntimes_ndep 
  • branches/2016/dev_r7012_ROBUST5_CNRS/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zsink.F90

    r7162 r7192  
    136136               IF( tmask(ji,jj,jk) == 1 ) THEN 
    137137                 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 ) ) 
    140140               ENDIF 
    141141            END DO 
     
    194194                  IF( tmask(ji,jj,jk) == 1 ) THEN 
    195195                    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 ) ) 
    197197                  ENDIF 
    198198               END DO 
     
    310310      CALL wrk_alloc( jpi, jpj, jpk, ztraz, zakz, zwsink2, ztrb ) 
    311311 
    312       zstep = rfact2 / FLOAT( kiter ) / 2. 
     312      zstep = rfact2 / REAL( kiter, wp ) / 2. 
    313313 
    314314      ztraz(:,:,:) = 0.e0 
  • branches/2016/dev_r7012_ROBUST5_CNRS/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zsms.F90

    r7162 r7192  
    8888      IF( ( ln_top_euler .AND. kt == nittrc000 )  .OR. ( .NOT.ln_top_euler .AND. kt <= nittrc000 + nn_dttrc ) ) THEN 
    8989         rfactr  = 1. / rfact 
    90          rfact2  = rfact / FLOAT( nrdttrc ) 
     90         rfact2  = rfact / REAL( nrdttrc, wp ) 
    9191         rfact2r = 1. / rfact2 
    9292         xstep = rfact2 / rday         ! Time step duration for biology 
  • branches/2016/dev_r7012_ROBUST5_CNRS/NEMOGCM/NEMO/TOP_SRC/TRP/trcrad.F90

    r7068 r7192  
    214214            IF( l_trdtrc ) THEN 
    215215               ! 
    216                zs2rdt = 1. / ( 2. * rdt * FLOAT( nn_dttrc ) ) 
     216               zs2rdt = 1. / ( 2. * rdt * REAL( nn_dttrc, wp ) ) 
    217217               ztrtrdb(:,:,:) = ( ptrb(:,:,:,jn) - ztrtrdb(:,:,:) ) * zs2rdt 
    218218               ztrtrdn(:,:,:) = ( ptrn(:,:,:,jn) - ztrtrdn(:,:,:) ) * zs2rdt  
Note: See TracChangeset for help on using the changeset viewer.