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 8127 – NEMO

Changeset 8127


Ignore:
Timestamp:
2017-06-02T09:54:18+02:00 (7 years ago)
Author:
cbricaud
Message:

fix in nemo_v3_6_STABLE ticket #1906

Location:
branches/2015/nemo_v3_6_STABLE/NEMOGCM
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • branches/2015/nemo_v3_6_STABLE/NEMOGCM/NEMO/OPA_SRC/DIA/diadct.F90

    r5505 r8127  
    249249           !debug this section computing ? 
    250250           lldebug=.FALSE. 
    251            IF( (jsec==nn_secdebug .OR. nn_secdebug==-1) .AND.  kt==nit000+nn_dct-1 .AND. lwp ) lldebug=.TRUE.  
     251           IF( (jsec==nn_secdebug .OR. nn_secdebug==-1) .AND.  kt==nit000+nn_dct-1 ) lldebug=.TRUE.  
    252252 
    253253           !Compute transport through section   
     
    258258        IF( MOD(kt,nn_dctwri)==0 )THEN 
    259259 
    260            IF( lwp .AND. kt==nit000+nn_dctwri-1 )WRITE(numout,*)"      diadct: average transports and write at kt = ",kt          
     260           IF( kt==nit000+nn_dctwri-1 )WRITE(numout,*)"      diadct: average transports and write at kt = ",kt          
    261261   
    262262           !! divide arrays by nn_dctwri/nn_dct to obtain average  
     
    344344     DO jsec=1,nb_sec_max      !loop on the nb_sec sections 
    345345 
    346         IF ( lwp .AND. ( jsec==nn_secdebug .OR. nn_secdebug==-1 ) ) & 
     346        IF ( jsec==nn_secdebug .OR. nn_secdebug==-1 ) & 
    347347           & WRITE(numout,*)'debuging for section number: ',jsec  
    348348 
     
    364364        IF( jsec .NE. isec )  CALL ctl_stop( cltmp ) 
    365365 
    366         IF( lwp .AND. ( jsec==nn_secdebug .OR. nn_secdebug==-1 ) )WRITE(numout,*)"isec ",isec  
     366        IF( jsec==nn_secdebug .OR. nn_secdebug==-1 )WRITE(numout,*)"isec ",isec  
    367367 
    368368        READ(numdct_in)secs(jsec)%name 
     
    383383        !----- 
    384384 
    385         IF( lwp .AND. ( jsec==nn_secdebug .OR. nn_secdebug==-1 ) )THEN 
     385        IF( jsec==nn_secdebug .OR. nn_secdebug==-1 )THEN 
    386386           
    387387            WRITE(clformat,'(a,i2,a)') '(A40,', nb_class_max,'(f8.3,1X))'  
     
    416416           !debug 
    417417           !----- 
    418            IF( lwp .AND. ( jsec==nn_secdebug .OR. nn_secdebug==-1 ) )THEN 
     418           IF( jsec==nn_secdebug .OR. nn_secdebug==-1 )THEN 
    419419              WRITE(numout,*)"      List of points in global domain:" 
    420420              DO jpt=1,iptglo 
     
    450450           !debug 
    451451           !----- 
    452            IF(   lwp .AND. ( jsec==nn_secdebug .OR. nn_secdebug==-1 ) )THEN 
     452           IF( jsec==nn_secdebug .OR. nn_secdebug==-1 )THEN 
    453453              WRITE(numout,*)"      List of points selected by the proc:" 
    454454              DO jpt = 1,iptloc 
     
    468468           !remove redundant points between processors 
    469469           !------------------------------------------ 
    470            lldebug = .FALSE. ; IF ( (jsec==nn_secdebug .OR. nn_secdebug==-1) .AND. lwp ) lldebug = .TRUE. 
     470           lldebug = .FALSE. ; IF ( jsec==nn_secdebug .OR. nn_secdebug==-1 ) lldebug = .TRUE. 
    471471           IF( iptloc .NE. 0 )THEN 
    472472              CALL removepoints(secs(jsec),'I','top_list',lldebug) 
     
    484484           !debug 
    485485           !----- 
    486            IF( lwp .AND. ( jsec==nn_secdebug .OR. nn_secdebug==-1 ) )THEN 
     486           IF( jsec==nn_secdebug .OR. nn_secdebug==-1 )THEN 
    487487              WRITE(numout,*)"      List of points after removepoints:" 
    488488              iptloc = secs(jsec)%nb_point 
     
    496496 
    497497        ELSE  ! iptglo = 0 
    498            IF( lwp .AND. ( jsec==nn_secdebug .OR. nn_secdebug==-1 ) )& 
     498           IF( jsec==nn_secdebug .OR. nn_secdebug==-1 )& 
    499499              WRITE(numout,*)'   No points for this section.' 
    500500        ENDIF 
  • branches/2015/nemo_v3_6_STABLE/NEMOGCM/TOOLS/SECTIONS_DIADCT/src/compute_sections.f90

    r3632 r8127  
    681681     IF( sec%nb_point .ne. 0 )THEN 
    682682        IF (  sec%listPoint(sec%nb_point)%I .NE.  sec%listPoint(1)%I ) THEN 
    683            sec%slopeSection = ( sec%listPoint(sec%nb_point)%J - sec%listPoint(1)%J ) /  & 
    684                               ( sec%listPoint(sec%nb_point)%I - sec%listPoint(1)%I )       
     683           sec%slopeSection = REAL( ( sec%listPoint(sec%nb_point)%J - sec%listPoint(1)%J ) , wp )/  & 
     684                              REAL( ( sec%listPoint(sec%nb_point)%I - sec%listPoint(1)%I ) , wp ) 
     685 
    685686        ELSE 
    686687           sec%slopeSection = 10000._wp 
Note: See TracChangeset for help on using the changeset viewer.