New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 1264 – NEMO

Changeset 1264


Ignore:
Timestamp:
2009-01-13T19:54:08+01:00 (15 years ago)
Author:
cetlod
Message:

clean TOP model routines to avoid warning when compiling, see ticket:303

Location:
trunk/NEMO/TOP_SRC
Files:
24 edited

Legend:

Unmodified
Added
Removed

  • trunk/NEMO/TOP_SRC/LOBSTER/trcbio.F90

    r1255 r1264  
    6161      INTEGER, INTENT( in ) ::   kt      ! ocean time-step index       
    6262      !! 
    63       INTEGER  ::   ji, jj, jk, jn, jl 
     63      INTEGER  ::   ji, jj, jk, jl 
    6464      REAL(wp) ::   zdet, zzoo, zphy, zno3, znh4, zdom      ! now concentrations 
    6565      REAL(wp) ::   zlno3, zlnh4, zle, zlt                  ! limitation terms for phyto 

  • trunk/NEMO/TOP_SRC/MY_TRC/trcctl_my_trc.F90

    r1255 r1264  
    6161   !!====================================================================== 
    6262END MODULE trcctl_my_trc 
    63        

  • trunk/NEMO/TOP_SRC/PISCES/p4zbio.F90

    r1152 r1264  
    8282           
    8383      CALL p4z_sink ( kt, jnt )     ! vertical flux of particulate organic matter 
    84  
    8584      CALL p4z_opt  ( kt, jnt )     ! Optic: PAR in the water column 
    86  
    8785      CALL p4z_lim  ( kt, jnt )     ! co-limitations by the various nutrients 
    88  
    8986      CALL p4z_prod ( kt, jnt )     ! phytoplankton growth rate over the global ocean.  
    9087      !                             ! (for each element : C, Si, Fe, Chl ) 
    91  
    9288      CALL p4z_rem  ( kt, jnt )     ! remineralization terms of organic matter+scavenging of Fe 
    93  
    9489      CALL p4z_mort ( kt, jnt )     ! phytoplankton mortality 
    95  
    9690      !                             ! zooplankton sources/sinks routines  
    9791      CALL p4z_micro( kt, jnt )           ! microzooplankton 

  • trunk/NEMO/TOP_SRC/PISCES/p4zsink.F90

    r1180 r1264  
    3232     sinkfer                   !: Small BFe sinking flux 
    3333 
     34   REAL(wp) ::   & 
     35     xstep , xstep2            !: Time step duration for biology 
     36 
    3437#if  defined key_kriest 
    3538   REAL(wp)          ::       &    
    … …  
    9295      REAL(wp) :: zdiv , zdiv1, zdiv2, zdiv3, zdiv4, zdiv5 
    9396      REAL(wp) :: zval1, zval2, zval3, zval4 
    94       REAL(wp) :: zstep 
    9597#if defined key_trc_dia3d 
    9698      REAL(wp) ::   zrfact2 
    … …  
    101103      !!--------------------------------------------------------------------- 
    102104 
    103       IF( ( kt * jnt ) == nittrc000  )   CALL p4z_sink_init   ! Initialization (first time-step only) 
    104  
    105        zstep = rfact2 / rjjss      ! Time step duration for biology 
    106  
     105      IF( ( kt * jnt ) == nittrc000  )  THEN  
     106          CALL p4z_sink_init   ! Initialization (first time-step only) 
     107          xstep  = rfact2 / rjjss      ! Time step duration for biology 
     108          xstep2 = rfact2 /  2. 
     109      ENDIF 
    107110 
    108111!     Initialisation of variables used to compute Sinking Speed 
    … …  
    252255                     &    ) 
    253256 
    254                   zaggsi = ( zagg4 + zagg5 ) * zstep / 10. 
     257                  zaggsi = ( zagg4 + zagg5 ) * xstep / 10. 
    255258 
    256259                  zagg = 0.5 * xkr_stick * ( zaggsh + zaggsi ) 
    … …  
    260263 
    261264                  zaggdoc = ( 0.4 * trn(ji,jj,jk,jpdoc)               & 
    262                      &        + 1018.  * trn(ji,jj,jk,jppoc)  ) * zstep    & 
     265                     &        + 1018.  * trn(ji,jj,jk,jppoc)  ) * xstep    & 
    263266# if defined key_off_degrad 
    264267                     &        * facvol(ji,jj,jk)                              & 
    … …  
    450453      REAL(wp) ::   zagg1, zagg2, zagg3, zagg4 
    451454      REAL(wp) ::   zagg , zaggfe, zaggdoc, zaggdoc2 
    452       REAL(wp) ::   zfact, zstep, zwsmax 
     455      REAL(wp) ::   zfact, zwsmax 
    453456#if defined key_trc_dia3d 
    454457      REAL(wp) ::   zrfact2 
    … …  
    457460      !!--------------------------------------------------------------------- 
    458461 
    459        zstep = rfact2 / rjjss      ! Timestep duration for biology 
    460  
     462      IF( ( kt * jnt ) == nittrc000  )  THEN 
     463          xstep  = rfact2 / rjjss      ! Timestep duration for biology 
     464          xstep2 = rfact2 /  2. 
     465      ENDIF 
    461466 
    462467!    Sinking speeds of detritus is increased with depth as shown 
    … …  
    490495         DO jj = 1, jpj 
    491496            DO ji = 1, jpi 
    492                zwsmax = 0.8 * fse3t(ji,jj,jk) / zstep 
     497               zwsmax = 0.8 * fse3t(ji,jj,jk) / xstep 
    493498               wsbio4(ji,jj,jk) = MIN( wsbio4(ji,jj,jk), zwsmax ) 
    494499               wsbio3(ji,jj,jk) = MIN( wsbio3(ji,jj,jk), zwsmax ) 
    … …  
    527532         DO jj = 1, jpj 
    528533            DO ji = 1, jpi 
    529                zfact = zstep * xdiss(ji,jj,jk) 
     534               zfact = xstep * xdiss(ji,jj,jk) 
    530535 
    531536!    Part I : Coagulation dependent on turbulence 
    … …  
    549554 
    550555# if defined key_off_degrad 
    551                zagg3 = 0.66 * zstep * trn(ji,jj,jk,jppoc) * trn(ji,jj,jk,jpgoc) * facvol(ji,jj,jk) 
     556               zagg3 = 0.66 * xstep * trn(ji,jj,jk,jppoc) * trn(ji,jj,jk,jpgoc) * facvol(ji,jj,jk) 
    552557# else 
    553                zagg3 = 0.66 * zstep * trn(ji,jj,jk,jppoc) * trn(ji,jj,jk,jpgoc) 
    554 # endif 
    555  
    556 # if defined key_off_degrad 
    557                zagg4 = 0.e0 * zstep * trn(ji,jj,jk,jppoc) * trn(ji,jj,jk,jppoc) * facvol(ji,jj,jk) 
     558               zagg3 = 0.66 * xstep * trn(ji,jj,jk,jppoc) * trn(ji,jj,jk,jpgoc) 
     559# endif 
     560 
     561# if defined key_off_degrad 
     562               zagg4 = 0.e0 * xstep * trn(ji,jj,jk,jppoc) * trn(ji,jj,jk,jppoc) * facvol(ji,jj,jk) 
    558563# else 
    559                zagg4 = 0.e0 * zstep * trn(ji,jj,jk,jppoc) * trn(ji,jj,jk,jppoc) 
     564               zagg4 = 0.e0 * xstep * trn(ji,jj,jk,jppoc) * trn(ji,jj,jk,jppoc) 
    560565# endif 
    561566 
    … …  
    627632      !! 
    628633      INTEGER  ::   ji, jj, jk, jn 
    629       REAL(wp) ::   zigma,zew,zstep,zign, zflx 
     634      REAL(wp) ::   zigma,zew,zign, zflx 
    630635      REAL(wp), DIMENSION(jpi,jpj,jpk) ::  ztraz, zakz 
    631636      REAL(wp), DIMENSION(jpi,jpj,jpk) ::  zwsink2 
    632637      !!--------------------------------------------------------------------- 
    633638 
    634       zstep  = rfact2 / 2. 
    635639 
    636640      ztraz(:,:,:) = 0.e0 
    … …  
    680684            DO jj = 1, jpj       
    681685               DO ji = 1, jpi     
    682                   zigma = zwsink2(ji,jj,jk+1) * zstep / fse3w(ji,jj,jk+1) 
     686                  zigma = zwsink2(ji,jj,jk+1) * xstep2 / fse3w(ji,jj,jk+1) 
    683687                  zew   = zwsink2(ji,jj,jk+1) 
    684                   psinkflx(ji,jj,jk+1) = -zew * ( trn(ji,jj,jk,jp_tra) - 0.5 * ( 1 + zigma ) * zakz(ji,jj,jk) ) * zstep 
     688                  psinkflx(ji,jj,jk+1) = -zew * ( trn(ji,jj,jk,jp_tra) - 0.5 * ( 1 + zigma ) * zakz(ji,jj,jk) ) * xstep2 
    685689               END DO 
    686690            END DO 

  • trunk/NEMO/TOP_SRC/PISCES/sms_pisces.F90

    r1180 r1264  
    8080   !!======================================================================    
    8181END MODULE sms_pisces     
    82  

  • trunk/NEMO/TOP_SRC/PISCES/trcctl_pisces.F90

    r1146 r1264  
    6868   !!====================================================================== 
    6969END MODULE trcctl_pisces 
    70        

  • trunk/NEMO/TOP_SRC/PISCES/trcini_pisces.F90

    r1146 r1264  
    5454      !! ** Purpose :   Initialisation of the PISCES biochemical model 
    5555      !!---------------------------------------------------------------------- 
    56       INTEGER ::   ji, jj, jk, jn 
     56      INTEGER ::   ji, jj, jk 
    5757      REAL(wp) ::  caralk, bicarb, co3 
    5858 

  • trunk/NEMO/TOP_SRC/SED/sedadv.F90

    r1179 r1264  
    457457#endif 
    458458END MODULE sedadv 
    459                    

  • trunk/NEMO/TOP_SRC/SED/sedco3.F90

    r1179 r1264  
    211211 
    212212END MODULE sedco3 
    213  

  • trunk/NEMO/TOP_SRC/SED/seddta.F90

    r1250 r1264  
    275275   SUBROUTINE sed_dta ( kt ) 
    276276     INTEGER, INTENT(in) :: kt 
     277     WRITE(*,*) 'sed_stp: You should not have seen this print! error?', kt  
    277278  END SUBROUTINE sed_dta 
    278279#endif 

  • trunk/NEMO/TOP_SRC/SED/sedini.F90

    r1250 r1264  
    861861 
    862862END MODULE sedini 
    863  
    864  
    865  

  • trunk/NEMO/TOP_SRC/SED/sedmbc.F90

    r1250 r1264  
    334334#endif 
    335335END MODULE sedmbc 
    336  
    337  
    338  
    339  
    340  

  • trunk/NEMO/TOP_SRC/SED/sedmodel.F90

    r1250 r1264  
    5050#else 
    5151   !!====================================================================== 
    52    !! MODULE sedbtb  :   Dummy module  
     52   !! MODULE sedmodel  :   Dummy module  
    5353   !!====================================================================== 
    5454CONTAINS 
    … …  
    6060 
    6161END MODULE sedmodel 
    62  
    63  
    64  
    65  
    66  
    67  
    68  

  • trunk/NEMO/TOP_SRC/SED/sedstp.F90

    r1250 r1264  
    5858      CALL sed_sfc( kt )         ! Give back new bottom wat chem to tracer model 
    5959#endif 
    60       CALL sed_rst_write( kt )   ! restart file output 
     60      CALL sed_rst_wri( kt )   ! restart file output 
    6161      CALL sed_wri( kt )         ! outputs 
    6262 
    … …  
    7575   END SUBROUTINE sed_stp 
    7676#endif 
    77  
    7877END MODULE sedstp 
    79  
    80  
    81  
    82  
    83  
    84  
    85  

  • trunk/NEMO/TOP_SRC/TRP/trcadv_muscl.F90

    r1193 r1264  
    8080      REAL(wp) ::   z0u, z0v, z0w 
    8181      REAL(wp) ::   zzt1, zzt2, zalpha, z2dtt 
    82       REAL(wp) ::   ztai, ztaj, zfui, zfvj 
     82      REAL(wp) ::   zfui, zfvj 
    8383      REAL(wp), DIMENSION(:,:,:), ALLOCATABLE ::   ztrtrd 
     84#if defined key_trc_diatrd 
     85      REAL(wp) ::   ztai, ztaj 
     86#endif 
    8487      CHARACTER (len=22) :: charout 
    8588      !!---------------------------------------------------------------------- 

  • trunk/NEMO/TOP_SRC/TRP/trcadv_muscl2.F90

    r1193 r1264  
    8282      REAL(wp) ::   z0u, z0v, z0w 
    8383      REAL(wp) ::   zzt1, zzt2, zalpha 
    84  
    85       REAL(wp) ::   ztai, ztaj 
    8684      REAL(wp) ::   zfui, zfvj 
    8785      REAL(wp), DIMENSION(:,:,:), ALLOCATABLE ::   ztrtrd 
     86#if defined key_trc_diatrd 
     87      REAL(wp) ::   ztai, ztaj 
     88#endif 
    8889      CHARACTER (len=22) :: charout 
    8990      !!---------------------------------------------------------------------- 

  • trunk/NEMO/TOP_SRC/TRP/trcadv_tvd.F90

    r1193 r1264  
    7777      REAL(wp) ::   zfp_ui, zfp_vj, zfp_wk 
    7878      REAL(wp) ::   zfm_ui, zfm_vj, zfm_wk 
     79#if defined key_trc_diatrd 
    7980      REAL(wp) ::   zgm, zgz 
     81#endif 
    8082      CHARACTER (len=22) :: charout 
    8183      !!---------------------------------------------------------------------- 

  • trunk/NEMO/TOP_SRC/TRP/trcbbl.F90

    r1198 r1264  
    107107      INTEGER ::   ii0, ii1, ij0, ij1                       ! temporary integers 
    108108      INTEGER ::   iku1, iku2, ikv1, ikv2                   ! temporary intergers 
    109       REAL(wp) ::   ze3u, ze3v                              ! temporary scalars 
    110       REAL(wp) ::   zsign, zt, zs, zh, zalbet 
    111       REAL(wp) ::   zgdrho, zbtr, ztra 
    112       REAL(wp), DIMENSION(jpi,jpj) ::   zki, zkj, zkx, zky  ! temporary workspace arrays 
     109      REAL(wp) ::  ze3u, ze3v                              ! temporary scalars 
     110      REAL(wp) ::  zbtr, ztra 
     111#if ! defined key_off_tra 
     112      REAL(wp) ::   zgdrho, zalbet, zsign, zt, zs, zh 
     113      REAL(wp), DIMENSION(jpi,jpj) ::   zki, zkj 
     114#endif 
     115      REAL(wp), DIMENSION(jpi,jpj) ::   zkx, zky  ! temporary workspace arrays 
    113116      REAL(wp), DIMENSION(jpi,jpj) ::   ztnb, zsnb, zdep 
    114117      REAL(wp), DIMENSION(jpi,jpj) ::   ztrb, zahu, zahv 

  • trunk/NEMO/TOP_SRC/TRP/trcldf_bilapg.F90

    r1175 r1264  
    6565      INTEGER, INTENT( in ) ::   kt                         ! ocean time-step index 
    6666      INTEGER ::   ji, jj, jk, jn                           ! dummy loop indices 
    67       REAL(wp) ::   ztra                                    ! workspace    
    6867      REAL(wp), DIMENSION(jpi,jpj,jpk,jptra) ::   wk1, wk2  ! work array used for rotated biharmonic 
    6968      !                                                     ! operator on tracers and/or momentum 

  • trunk/NEMO/TOP_SRC/TRP/trcldf_iso_zps.F90

    r1193 r1264  
    320320 
    321321         ENDIF 
    322 !CDIR END 
    323322         !                                                       ! =========== 
    324323      END DO                                                     ! tracer loop 

  • trunk/NEMO/TOP_SRC/TRP/trcldf_lap.F90

    r1175 r1264  
    7878 
    7979      REAL(wp) ::   & 
    80          ztra, ztrax, ztray               ! workspace 
     80         ztrax, ztray               ! workspace 
    8181      CHARACTER (len=22) :: charout 
    8282      REAL(wp), DIMENSION(:,:,:), ALLOCATABLE ::   ztrtrd 

  • trunk/NEMO/TOP_SRC/TRP/trczdf_imp.F90

    r1175 r1264  
    8585         zwd, zws, zwi,          &  ! ??? 
    8686         zwx, zwy, zwt              ! ??? 
     87#if defined key_trc_diatrd 
    8788      REAL(wp) ::  ztra      ! temporary scalars 
     89#endif 
    8890      REAL(wp), DIMENSION(jpi,jpj,jpk,jptra) ::   & 
    8991         ztrd 

  • trunk/NEMO/TOP_SRC/TRP/trczdf_iso_vopt.F90

    r1258 r1264  
    8585      ! ----------------------------------------------------- 
    8686 
    87       CALL trc_zdf_iso( kt ) 
     87      CALL trc_zdf_iso 
    8888 
    8989      IF( ln_ctl ) THEN    ! print mean trends (used for debugging) 
    … …  
    369369 
    370370 
    371    SUBROUTINE trc_zdf_iso( kt ) 
     371   SUBROUTINE trc_zdf_iso 
    372372      !!---------------------------------------------------------------------- 
    373373      !!                  ***  ROUTINE trc_zdf_iso  *** 
    … …  
    411411                            zwy => va      ! workspace arrays 
    412412 
    413       INTEGER, INTENT( in ) ::   kt              ! ocean time-step index 
    414413      INTEGER ::   ji, jj, jk, jn                ! dummy loop indices 
    415414      INTEGER ::   iku, ikv                       

  • trunk/NEMO/TOP_SRC/TRP/trdmld_trc.F90

    r1258 r1264  
    15981598      WRITE(*,*) '  "      "      : You should not have seen this print! error?', ctype 
    15991599      WRITE(*,*) '  "      "      : You should not have seen this print! error?', ktrd 
     1600      WRITE(*,*) '  "      "      : You should not have seen this print! error?', kjn 
    16001601   END SUBROUTINE trd_mld_trc_zint 
    16011602 
Note: See TracChangeset for help on using the changeset viewer.