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 2864 for branches/2011/dev_r2802_MERCATOR10_diadct/NEMOGCM/NEMO/OPA_SRC/DIA/diadct.F90 – NEMO

Ignore:
Timestamp:
2011-09-26T17:02:11+02:00 (13 years ago)
Author:
cbricaud
Message:

add dia_dct_init Dummy routine


File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2011/dev_r2802_MERCATOR10_diadct/NEMOGCM/NEMO/OPA_SRC/DIA/diadct.F90

    r2857 r2864  
    335335                 secs(jsec)%listPoint(iptloc) = POINT_SECTION(mi0(iiglo),mj0(ijglo)) ! store local coordinates 
    336336                 secs(jsec)%direction(iptloc) = directemp(jpt)               ! store local direction 
     337                 !WRITE(narea+200,*)'         # I J : ',iiloc,ijloc 
    337338              ENDIF 
    338339 
     
    343344           !debug 
    344345           !----- 
    345            IF( lwp .AND. ( jsec==nn_secdebug .OR. nn_secdebug==-1 ) )THEN 
     346           !IF( lwp .AND. ( jsec==nn_secdebug .OR. nn_secdebug==-1 ) )THEN 
     347           IF(  ( jsec==nn_secdebug .OR. nn_secdebug==-1 ) )THEN 
    346348              WRITE(numout,*)"      List of points selected by the proc:" 
    347349              DO jpt = 1,iptloc 
     
    349351                 ijglo = secs(jsec)%listPoint(jpt)%J + jpjzoom - 1 + njmpp - 1 
    350352                 WRITE(numout,*)'         # I J : ',iiglo,ijglo 
     353                 !WRITE(narea+200,*)'         # I J : ',iiglo,ijglo 
    351354              ENDDO 
    352355           ENDIF 
     
    582585           CASE(3)  ;   k = POINT_SECTION(sec%listPoint(jseg)%I,sec%listPoint(jseg)%J+1) 
    583586           END SELECT 
    584                
     587 
    585588           !------------------------------- 
    586589           !  LOOP ON THE DENSITY CLASSES | 
     
    620623                 !TEST ON THE DENSITY/SALINITY/TEMPERATURE/LEVEL!  
    621624                 !----------------------------------------------! 
     625           IF(ld_debug)write(narea+200,*)k ,jk 
    622626  
    623627                 IF ( (    ((( zrhop .GE. (sec%zsigp(jclass)+1000.  )) .AND.    & 
    624628                           (   zrhop .LE. (sec%zsigp(jclass+1)+1000. ))) .OR.    & 
    625629                           ( sec%zsigp(jclass) .EQ. 99.)) .AND.                 & 
    626                            ((( zrhoi .GE. (sec%zsigi(jclass) )) .AND.    & 
    627                            (   zrhoi .LE. (sec%zsigi(jclass+1)))) .OR.    & 
     630                           ((( zrhoi .GE. (sec%zsigi(jclass) + 1000. )) .AND.    & 
     631                           (   zrhoi .LE. (sec%zsigi(jclass+1)+1000. ))) .OR.    & 
    628632                           ( sec%zsigi(jclass) .EQ. 99.)) .AND.                 & 
    629633                           ((( zsn .GT. sec%zsal(jclass)) .AND.                & 
     
    659663                    ENDIF 
    660664#endif 
     665           !IF(ld_debug)write(narea+200,*)k ,jk,zTnorm              
    661666                    !COMPUTE TRANSPORT  
    662667                    !zTnorm=transport through one cell for one class 
     
    10541059   !!---------------------------------------------------------------------- 
    10551060   LOGICAL, PUBLIC, PARAMETER ::   lk_diadct = .FALSE.    !: diamht flag 
     1061   PUBLIC  
    10561062CONTAINS 
     1063 
     1064   SUBROUTINE dia_dct_init          ! Dummy routine 
     1065      WRITE(*,*) 'dia_dct_init: You should not have seen this print! error?', kt 
     1066   END SUBROUTINE dia_dct_init 
     1067 
    10571068   SUBROUTINE dia_dct( kt )           ! Dummy routine 
    10581069      INTEGER, INTENT( in ) ::   kt      ! ocean time-step index 
Note: See TracChangeset for help on using the changeset viewer.