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 239 for trunk/NEMO/OPA_SRC/DIA/diaptr.F90 – NEMO

Ignore:
Timestamp:
2005-03-22T18:55:20+01:00 (19 years ago)
Author:
opalod
Message:

CT : UPDATE172 : remove all direct acces modules and the related cpp key key_fdir

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMO/OPA_SRC/DIA/diaptr.F90

    r190 r239  
    345345   END SUBROUTINE dia_ptr_init 
    346346 
    347 #if defined key_fdir 
    348    !!--------------------------------------------------------------------- 
    349    !!   'key_fdir'                                      direct access file 
    350    !!--------------------------------------------------------------------- 
    351  
    352    SUBROUTINE dia_ptr_wri( kt ) 
    353       !!--------------------------------------------------------------------- 
    354       !!                ***  ROUTINE dia_ptr_wri  *** 
    355       !! 
    356       !! ** Purpose :   output of poleward fluxes  
    357       !!  
    358       !! ** Method  :   NetCDF file 
    359       !! 
    360       !! History : 
    361       !!   9.0  !  03-09  (G. Madec)  Original code   
    362       !!----------------------------------------------------------------------  
    363       !! * Arguments 
    364       INTEGER, INTENT(in) ::   kt   ! ocean time-step index 
    365       REAL(wp), DIMENSION(jpj) ::   zphi, zfoo 
    366       !!----------------------------------------------------------------------  
    367  
    368       IF( kt == nit000 ) THEN 
    369  
    370          ! Reference latitude 
    371          ! ------------------ 
    372          !                                           ! ======================= 
    373          IF( cp_cfg == "orca" ) THEN                 !   ORCA configurations 
    374             !                                        ! ======================= 
    375  
    376             IF( jp_cfg == 05  )   iline = 192   ! i-line that passes near the North Pole 
    377             IF( jp_cfg == 025 )   iline = 384   ! i-line that passes near the North Pole 
    378             IF( jp_cfg == 2   )   iline =  48   ! i-line that passes near the North Pole 
    379             IF( jp_cfg == 4   )   iline =  24   ! i-line that passes near the North Pole 
    380   
    381             zphi(:) = 0.e0    
    382             DO ji = mi0(iline), mi1(iline) 
    383                zphi(:) = gphiv(ji,:)         ! if iline is in the local domain 
    384             END DO 
    385             ! provide the correct zphi to all local domains 
    386             IF( lk_mpp )   CALL mpp_sum( zphi, jpj )         
    387             ! introduce arbitray northernmost grid point to avoid netcdf error 
    388             DO jj=mj0(jpjglo), mj1(jpjglo) 
    389                zphi(jj) = 2*zphi(jj-1)-zphi(jj-2) 
    390             ENDDO 
    391  
    392             !                                        ! ======================= 
    393          ELSE                                        !   OTHER configurations 
    394             !                                        ! ======================= 
    395             zphi(:) = gphiv(1,:)             ! assume lat/lon coordinate, select the first i-line 
    396             ! 
    397          ENDIF 
    398   
    399          ! open the output file          
    400          CALL ctlopn( numptr, 'opaptr.output', 'UNKNOWN', 'UNFORMATTED', 'SEQUENTIAL', 1, numout, lwp, 1 ) 
    401  
    402          ! header of output 
    403          WRITE( numptr ) cexper, no, zdt, nf_ptr, jpj, jpk, zphi 
    404  
    405       ENDIF 
    406  
    407       IF( MOD( kt, nf_ptr ) == 0 ) THEN 
    408             IF(lwp) WRITE( numptr ) kt, tn_jk, sn_jk, v_msf,              & 
    409 # if defined key_diaeiv 
    410                &                    pht_eiv, pst_eiv, v_msf_eiv,          & 
    411 # endif 
    412                &                    pht_adv, pht_ldf, pst_adv, pst_ldf 
    413       ENDIF 
    414  
    415       IF( kt == nitend )   CLOSE( numptr ) 
    416  
    417    END SUBROUTINE dia_ptr_wri 
    418  
    419 #else 
    420347   !!--------------------------------------------------------------------- 
    421348   !!   Default option :                                       NetCDF file 
     
    598525   END SUBROUTINE dia_ptr_wri 
    599526 
    600 #endif 
    601  
    602527   !!====================================================================== 
    603528END MODULE diaptr 
Note: See TracChangeset for help on using the changeset viewer.