Changeset 2878


Ignore:
Timestamp:
2011-09-30T08:26:28+02:00 (9 years ago)
Author:
cbricaud
Message:
 
Location:
branches/2011/dev_r2802_MERCATOR10_diadct/NEMOGCM/TOOLS/SECTIONS_DIADCT/src
Files:
7 edited

Legend:

Unmodified
Added
Removed
  • branches/2011/dev_r2802_MERCATOR10_diadct/NEMOGCM/TOOLS/SECTIONS_DIADCT/src/compute_sections.f90

    r2877 r2878  
    22   !!===================================================================== 
    33   !!                       ***  MODULE  diadct  *** 
    4    !! Ocean diagnostics: Compute the transport trough a sec. 
     4   !! Ocean diagnostics: Compute the transport trough a section 
    55   !! 
     6   !! History: 2011: Clement Bricaud, Mercator-Ocean 
    67   !! 
    7    !! History: 2011: cbricaud Mercator-Ocean 
    88   !!=============================================================== 
    99   !! * Modules used 
     
    301301           DO ji=1,jpi   
    302302              coord_t=COORD_SECTION(glamf(ji,jj),gphif(ji,jj)) 
    303 !cbr              zdistFirst = distance(coord_t,coordFirst) 
    304303              zdistFirst = distance2(coord_t,coordFirst) 
    305 !cbr              zdistLast = distance(coord_t,coordLast) 
    306304              zdistLast = distance2(coord_t,coordLast) 
    307305              IF( zdistFirst .LT. zdistmesh .AND. zdistFirst .LT. zdistante )THEN 
     
    370368 
    371369     IF( sec%nb_point .ne. 0 )THEN 
    372      !IF( nb_inmesh .NE. 0 )THEN 
    373370 
    374371        !The serie of mesh's points that form the section will 'link'  
     
    393390 
    394391        !initialize distnew value (with distance between section's extremities)  
    395 !cbr        zdistnew  = distance(coordFirst,coordLast,sec%ll_date_line)  
    396392        zdistnew  = distance2(coordFirst,coordLast,sec%ll_date_line)  
    397393        prevPoint  = POINT_SECTION(0,0) 
     
    596592                 CALL write_debug(jsec,cltmp) 
    597593                    
    598                  ! on passe à l'ancienne methode: le point parmies les 4 pts (NSWE)  qui se rapproche 
    599                  ! le + du dernier pt 
    600                  !-----------------------------  
    601594                 !be carreful! we can't go backward. 
    602  
    603595                 zdistNorth = zdistNorth2    ; zdistSouth = zdistSouth2 
    604596                 zdistEst   = zdistEst2      ; zdistWest  = zdistWest2  
  • branches/2011/dev_r2802_MERCATOR10_diadct/NEMOGCM/TOOLS/SECTIONS_DIADCT/src/declarations.f90

    r2858 r2878  
    22  !!---------------------------------------------------------------------------- 
    33  !!             *** declarations global varibles  
     4  !! 
     5  !! History: 2011: Clement Bricaud, Mercator-Ocean 
    46  !! 
    57  !!---------------------------------------------------------------------------- 
  • branches/2011/dev_r2802_MERCATOR10_diadct/NEMOGCM/TOOLS/SECTIONS_DIADCT/src/diadct_sections.f90

    r2858 r2878  
    22     !!============================================================================== 
    33     !!                       ***  PROGRAM generate_sections    *** 
     4     !! 
    45     !! create a binary file containig the IJ positions of sections in global  
    56     !! coordinates for the diagnostic routine diadct.F90 of NEMO 
    6      !! 
    77     !! 
    88     !! 
  • branches/2011/dev_r2802_MERCATOR10_diadct/NEMOGCM/TOOLS/SECTIONS_DIADCT/src/readcoordmesh.f90

    r2877 r2878  
    33  !!                       ***  MODULE  readcoordmesh  *** 
    44  !!  
    5   !! History: 2011: cbricaud Mercator-Ocean 
     5  !! History: 2011: Clement Bricaud, Mercator-Ocean 
    66  !!  
    77  !!===================================================================== 
  • branches/2011/dev_r2802_MERCATOR10_diadct/NEMOGCM/TOOLS/SECTIONS_DIADCT/src/readsections.f90

    r2877 r2878  
    33   !!                       ***  MODULE  readsections  *** 
    44   !!  
    5    !! History: 2011: cbricaud Mercator-Ocean 
     5   !! History: 2011: Clement Bricaud, Mercator-Ocean 
    66   !!  
    77   !!===================================================================== 
  • branches/2011/dev_r2802_MERCATOR10_diadct/NEMOGCM/TOOLS/SECTIONS_DIADCT/src/sections_tools.f90

    r2877 r2878  
    33   !!                       ***  MODULE  sections_tools  *** 
    44   !!  
    5    !! History: 2011: cbricaud Mercator-Ocean 
     5   !! History: 2011: Clement Bricaud, Mercator-Ocean 
    66   !!  
    77   !!===================================================================== 
     
    1414   !! * Routine accessibility 
    1515   PUBLIC pointToCoordF         ! define a point with geographical coordinates 
    16    PUBLIC distance              ! compute distance between 2 points    
    1716   PUBLIC distance2             ! compute distance between 2 points 
    1817   PUBLIC distance3             ! compute distance between a point and a line 
     
    4241 
    4342  END FUNCTION pointToCoordF 
    44  
    45   REAL(wp) FUNCTION distance(coordA,coordB,ld_date_line) 
    46      !!--------------------------------------------------------------------- 
    47      !!         ***  FUNCTION distance  *** 
    48      !! 
    49      !! ** Purpose: Compute distance between coordA and coordB 
    50      !!             We add 360° to coordB%long if the line (coordA,coordB) 
    51      !!             crosses the date-line. 
    52      !! 
    53      !!--------------------------------------------------------------------- 
    54      !! * arguments 
    55      TYPE(COORD_SECTION), INTENT(IN) :: coordA,coordB 
    56      LOGICAL,INTENT(IN),OPTIONAL     :: ld_date_line 
    57  
    58      !! * Local declarations  
    59      INTEGER ::idateline             ! local integer 
    60  
    61      !----------------------------------------------- 
    62      ! 
    63      idateline=0 
    64      IF( PRESENT( ld_date_line))THEN 
    65        IF( ld_date_line )idateline=1 
    66      ENDIF 
    67      ! 
    68      !function output 
    69      distance = (sqrt((coordA%lon-coordB%lon-360*idateline )**2 + (coordA%lat-coordB%lat )**2)) 
    70      ! 
    71      RETURN 
    72      ! 
    73   END FUNCTION distance 
    7443 
    7544  REAL(wp) FUNCTION distance2(coordA,coordB,ld_date_line) 
  • branches/2011/dev_r2802_MERCATOR10_diadct/NEMOGCM/TOOLS/SECTIONS_DIADCT/src/writesections.f90

    r2858 r2878  
    33   !!                       ***  MODULE  writesections  *** 
    44   !!  
    5    !! History: 2011: cbricaud Mercator-Ocean 
     5   !! History: 2011: Clement Bricaud, Mercator-Ocean 
    66   !!  
    77   !!===================================================================== 
Note: See TracChangeset for help on using the changeset viewer.