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

Changeset 2912


Ignore:
Timestamp:
2011-10-13T15:40:45+02:00 (12 years ago)
Author:
cbricaud
Message:

corrections after review

File:
1 edited

Legend:

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

    r2909 r2912  
    277277        READ(numdct_in,iostat=iost)isec 
    278278        IF (iost .NE. 0 )EXIT !end of file  
    279         WRITE(cltmp,'(a,i4.4,a,i4.4)')'diadct: read sections : Problem of section number: isec= ',isec,' et jsec= ',jsec 
     279        WRITE(cltmp,'(a,i4.4,a,i4.4)')'diadct: read sections : Problem of section number: isec= ',isec,' and jsec= ',jsec 
    280280        IF( jsec .NE. isec )  CALL ctl_stop( cltmp ) 
    281281 
     
    306306            WRITE(numout,*)                  "      Number of points in the section:  ",iptglo 
    307307            WRITE(numout,*)                  "      Number of classes                 ",secs(jsec)%nb_class 
    308             WRITE(numout,'(A40,10(f8.3,1X))')"      Insitu density classes :          ",secs(jsec)%zsigi 
     308!            WRITE(numout,'(A40,10(f8.3,1X))')"      Insitu density classes :          ",secs(jsec)%zsigi 
     309            WRITE(numout,'(A40,nb_class_max(f8.3,1X))')"      Insitu density classes :          ",secs(jsec)%zsigi 
    309310            WRITE(numout,'(A40,10(f8.3,1X))')"      Potential density classes :       ",secs(jsec)%zsigp 
    310311            WRITE(numout,'(A40,10(f8.3,1X))')"      Salinity classes :                ",secs(jsec)%zsal 
     
    353354                 iptloc = iptloc + 1                                                 ! count local points 
    354355                 secs(jsec)%listPoint(iptloc) = POINT_SECTION(mi0(iiglo),mj0(ijglo)) ! store local coordinates 
    355                  secs(jsec)%direction(iptloc) = directemp(jpt)               ! store local direction 
     356                 secs(jsec)%direction(iptloc) = directemp(jpt)                       ! store local direction 
    356357              ENDIF 
    357358 
     
    362363           !debug 
    363364           !----- 
    364            !IF( lwp .AND. ( jsec==nn_secdebug .OR. nn_secdebug==-1 ) )THEN 
    365            IF(  ( jsec==nn_secdebug .OR. nn_secdebug==-1 ) )THEN 
     365           IF(   lwp .AND. ( jsec==nn_secdebug .OR. nn_secdebug==-1 ) )THEN 
    366366              WRITE(numout,*)"      List of points selected by the proc:" 
    367367              DO jpt = 1,iptloc 
     
    369369                 ijglo = secs(jsec)%listPoint(jpt)%J + jpjzoom - 1 + njmpp - 1 
    370370                 WRITE(numout,*)'         # I J : ',iiglo,ijglo 
    371                  !WRITE(narea+200,*)'         # I J : ',iiglo,ijglo 
    372371              ENDDO 
    373372           ENDIF 
     
    396395        ELSE  ! iptglo = 0 
    397396           IF( lwp .AND. ( jsec==nn_secdebug .OR. nn_secdebug==-1 ) )& 
    398               WRITE(numout,*)'   Not points for this section.' 
     397              WRITE(numout,*)'   No points for this section.' 
    399398        ENDIF 
    400399 
     
    425424                itest         ,& !indice value of the side of the domain  
    426425                                 !where points could be redundant 
    427                 isgn          ,& !way of course in listpoint 
    428                 ipoint          ,& !way of course in listpoint 
     426                isgn          ,& ! isgn= 1 : scan listpoint from start to end 
     427                                 ! isgn=-1 : scan listpoint from end to start  
     428                ipoint        ,& !way of course in listpoint 
    429429                istart,iend      !first and last points selected in listpoint 
    430430     INTEGER :: jpoint   =0      !loop on list points 
    431      INTEGER,DIMENSION(nb_point_max)   :: idirec !contains temporare sec%direction 
    432      INTEGER,DIMENSION(2,nb_point_max) :: icoord !contains temporare sec%listpoint 
     431     INTEGER,DIMENSION(nb_point_max)   :: idirec !contains temporary sec%direction 
     432     INTEGER,DIMENSION(2,nb_point_max) :: icoord !contains temporary sec%listpoint 
    433433     !---------------------------------------------------------------------------- 
    434434     IF( ld_debug )WRITE(numout,*)'      -------------------------' 
     
    482482 
    483483  SUBROUTINE transport(sec,ld_debug) 
    484      !!--------------------------------------------------------------------- 
     484     !!------------------------------------------------------------------------------------------- 
    485485     !!                     ***  ROUTINE transport  *** 
    486486     !! 
     
    497497     !!              test on the density/temperature/salinity/level 
    498498     !! 
    499      !! ** Output: sec%transport: transport in the 2 direction and temperature,  
    500      !!                       salinity, density meaned by the transport 
    501      !! 
    502      !! 
    503      !!------------------------------------------------------------------- 
     499     !! ** Output: sec%transport: volume/mass/ice/heat/salt transport in the 2 directions 
     500     !! 
     501     !! 
     502     !!------------------------------------------------------------------------------------------- 
    504503     !! * Arguments 
    505504     TYPE(SECTION),INTENT(INOUT) :: sec 
     
    513512                zumid_ice    , zvmid_ice    ,&!U/V ice velocity 
    514513                zTnorm                      ,&!transport of velocity through one cell's sides 
    515                 ztransp1     , ztransp2     ,&!total        transport in sens 1 and 2 
     514                ztransp1     , ztransp2     ,&!total        transport in directions 1 and 2 
    516515                ztemp1       , ztemp2       ,&!temperature  transport     " 
    517516                zrhoi1       , zrhoi2       ,&!mass         transport     " 
    518517                zrhop1       , zrhop2       ,&!mass         transport     " 
    519518                zsal1        , zsal2        ,&!salinity     transport     " 
    520                 zice_vol_pos , zice_vol_neg ,&!volumic  ice transport     " 
    521                 zice_surf_pos, zice_surf_neg  !surfacic ice transport     " 
     519                zice_vol_pos , zice_vol_neg ,&!volume  ice transport     " 
     520                zice_surf_pos, zice_surf_neg  !surface ice transport     " 
    522521     REAL(wp):: ztn, zsn, zrhoi, zrhop, zsshn, zfsdep ! temperature/salinity/ssh/potential density /depth at u/v point 
    523522     REAL(wp) :: aj,i0,j0,i1,j1,i,j 
     
    545544        ! 
    546545        !convention: 
    547         !   non horizontal section: sens + is toward left hand of section 
    548         !       horizontal section: sens + is toward north of section 
     546        !   non horizontal section: direction + is toward left hand of section 
     547        !       horizontal section: direction + is toward north of section 
    549548        ! 
    550549        ! 
     
    552551        !       ----------------      -----------------     ---------------             -------------- 
    553552        ! 
    554         !   isgnv=1         sens +       
    555         !  ______         _____             ______(i1,j1)      (i1,j1)                                       
    556         !(i0,j0) |           //|            |                  |                   sens +    
    557         !        | isgnu=1  // |            |isgnu=1           |isgnu=1              /|\ 
    558         !        |_______  //         ______|    \\            | ---\                 | 
    559         !               |             | isgnv=-1  \\ |         | ---/ sens +       ____________(i1,j1)  
    560         !               |             |          __\\|         |               (i0,j0) 
    561         !               |(i1,j1)      |(i0,j0)   sens +        |                      isgnv=1                                  
    562         !                                                      (i0,j0) 
     553        !   isgnv=1         direction +       
     554        !  ______         _____             ______                                                    
     555        !        |           //|            |                  |                         direction +    
     556        !        | isgnu=1  // |            |isgnu=1           |isgnu=1                     /|\ 
     557        !        |_______  //         ______|    \\            | ---\                        | 
     558        !               |             | isgnv=-1  \\ |         | ---/ direction +       ____________ 
     559        !               |             |          __\\|         |                     
     560        !               |             |     direction +        |                      isgnv=1                                  
     561        !                                                       
    563562        !---------------------------------------------------------------------------------------------------- 
    564563        isgnu = 1 
     
    576575           !------------------------------------------ 
    577576           ! Select good coordinate to have velocity of the segment 
     577           ! 
     578           ! CASE( sec%direction(jseg) : 
    578579           ! 
    579580           !              |U CASE(3) 
     
    647648 
    648649 
    649                     !compute velocity with the good sens 
     650                    !compute velocity with the correct direction 
    650651                    SELECT CASE( sec%direction(jseg) ) 
    651652                    CASE(0,1)   
     
    657658                    END SELECT 
    658659 
    659                     !velocity* cell's lenght * cell's thickness 
     660                    !velocity* cell's length * cell's thickness 
    660661                    zTnorm=zumid*e2u(k%I,k%J)*  fse3u(k%I,k%J,jk)+     & 
    661662                           zvmid*e1v(k%I,k%J)*  fse3v(k%I,k%J,jk) 
     
    769770 
    770771 
    771      ELSE  !if nb_inmesh=0 
     772     ELSE  !if sec%nb_point =0 
    772773        zsum(1:2,:)=0. 
    773774        IF (sec%llstrpond) zsum(3:10,:)=0. 
    774775        zsum( 11:14,:)=0. 
    775      ENDIF   !end of nb_inmesh=0 case 
     776     ENDIF   !end of sec%nb_point =0 case 
    776777 
    777778     !-------------------------------| 
Note: See TracChangeset for help on using the changeset viewer.