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 5282 for branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/OPA_SRC/DIA – NEMO

Ignore:
Timestamp:
2015-05-18T17:19:50+02:00 (9 years ago)
Author:
diovino
Message:

Dev. branch CMCC4_simplification ticket #1456

Location:
branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/OPA_SRC/DIA
Files:
2 deleted
6 edited

Legend:

Unmodified
Added
Removed

  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/OPA_SRC/DIA/diaar5.F90

    r4990 r5282  
    2121   USE timing         ! preformance summary 
    2222   USE wrk_nemo       ! working arrays 
     23   USE fldread        ! type FLD_N 
     24   USE phycst         ! physical constant 
     25   USE in_out_manager  ! I/O manager 
    2326 
    2427   IMPLICIT NONE 
    … …  
    197200      REAL(wp) ::   zztmp   
    198201      REAL(wp), POINTER, DIMENSION(:,:,:,:) ::   zsaldta   ! Jan/Dec levitus salinity 
     202      ! reading initial file 
     203      LOGICAL  ::   ln_tsd_init      !: T & S data flag 
     204      LOGICAL  ::   ln_tsd_tradmp    !: internal damping toward input data flag 
     205      CHARACTER(len=100)            ::   cn_dir 
     206      TYPE(FLD_N)                   ::  sn_tem,sn_sal 
     207      INTEGER  ::   ios=0 
     208 
     209      NAMELIST/namtsd/ ln_tsd_init,ln_tsd_tradmp,cn_dir,sn_tem,sn_sal 
     210      ! 
     211 
     212      REWIND( numnam_ref )              ! Namelist namtsd in reference namelist : 
     213      READ  ( numnam_ref, namtsd, IOSTAT = ios, ERR = 901) 
     214901   IF( ios /= 0 ) CALL ctl_nam ( ios , ' namtsd in reference namelist for dia_ar5', lwp ) 
     215      REWIND( numnam_cfg )              ! Namelist namtsd in configuration namelist : Parameters of the run 
     216      READ  ( numnam_cfg, namtsd, IOSTAT = ios, ERR = 902 ) 
     217902   IF( ios /= 0 ) CALL ctl_nam ( ios , ' namtsd in configuration namelist for dia_ar5', lwp ) 
     218      IF(lwm) WRITE ( numond, namtsd ) 
     219      ! 
    199220      !!---------------------------------------------------------------------- 
    200221      ! 
    … …  
    216237      END DO 
    217238      IF( lk_mpp )   CALL mpp_sum( vol0 ) 
    218        
    219       CALL iom_open ( 'data_1m_salinity_nomask', inum ) 
    220       CALL iom_get  ( inum, jpdom_data, 'vosaline', zsaldta(:,:,:,1), 1  ) 
    221       CALL iom_get  ( inum, jpdom_data, 'vosaline', zsaldta(:,:,:,2), 12 ) 
     239 
     240      CALL iom_open ( TRIM( cn_dir )//TRIM(sn_sal%clname), inum ) 
     241      CALL iom_get  ( inum, jpdom_data, TRIM(sn_sal%clvar), zsaldta(:,:,:,1), 1  ) 
     242      CALL iom_get  ( inum, jpdom_data, TRIM(sn_sal%clvar), zsaldta(:,:,:,2), 12 ) 
    222243      CALL iom_close( inum ) 
    223244      sn0(:,:,:) = 0.5_wp * ( zsaldta(:,:,:,1) + zsaldta(:,:,:,2) )         

  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/OPA_SRC/DIA/diafwb.F90

    r4990 r5282  
    115115          
    116116         ! Conversion in m3 
    117          a_fwf    = a_fwf * rdttra(1) * 1.e-3  
     117         a_fwf    = a_fwf * rdt * 1.e-3 
    118118          
    119119         ! fwf correction to bring back the mean ssh to zero 
    … …  
    382382         WRITE(inum,*) 
    383383         WRITE(inum,*)    'Net freshwater budget ' 
    384          WRITE(inum,9010) '  fwf    = ',a_fwf,   ' m3 =', a_fwf   /(FLOAT(nitend-nit000+1)*rdttra(1)) * 1.e-6,' Sv' 
     384         WRITE(inum,9010) '  fwf    = ',a_fwf,   ' m3 =', a_fwf   /(FLOAT(nitend-nit000+1)*rdt) * 1.e-6,' Sv' 
    385385         WRITE(inum,*) 
    386386         WRITE(inum,9010) '  zarea =',zarea 

  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/OPA_SRC/DIA/diaharm.F90

    r4990 r5282  
    2323   USE ioipsl          ! NetCDF IPSL library 
    2424   USE lbclnk          ! ocean lateral boundary conditions (or mpp link) 
    25    USE diadimg         ! To write dimg 
    2625   USE timing          ! preformance summary 
    2726   USE wrk_nemo        ! working arrays 
    … …  
    388387      !!---------------------------------------------------------------------- 
    389388 
    390 #if defined key_dimgout 
    391       cdfile_name_T = TRIM(cexper)//'_Tidal_harmonics_gridT.dimgproc' 
    392       cdfile_name_U = TRIM(cexper)//'_Tidal_harmonics_gridU.dimgproc' 
    393       cdfile_name_V = TRIM(cexper)//'_Tidal_harmonics_gridV.dimgproc' 
    394 #endif 
     389!#if  
     390!#endif 
    395391 
    396392      IF(lwp) WRITE(numout,*) '  ' 
    397393      IF(lwp) WRITE(numout,*) 'dia_wri_harm : Write harmonic analysis results' 
    398 #if defined key_dimgout 
    399       IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~~  Output files: ', TRIM(cdfile_name_T) 
    400       IF(lwp) WRITE(numout,*) '                             ', TRIM(cdfile_name_U) 
    401       IF(lwp) WRITE(numout,*) '                             ', TRIM(cdfile_name_V) 
    402 #endif 
     394!#if  
     395!#endif 
    403396      IF(lwp) WRITE(numout,*) '  ' 
    404397 
    … …  
    406399      !///////////// 
    407400      ! 
    408 #if defined key_dimgout 
    409       cltext='Elevation amplitude and phase' 
    410       CALL dia_wri_dimg(TRIM(cdfile_name_T), TRIM(cltext), out_eta, 2*nb_ana, '2') 
    411 #else 
     401!#if  
     402!#else 
    412403      DO jh = 1, nb_ana 
    413404      CALL iom_put( TRIM(tname(jh))//'x', out_eta(:,:,jh) ) 
    414405      CALL iom_put( TRIM(tname(jh))//'y', out_eta(:,:,nb_ana+jh) ) 
    415406      END DO 
    416 #endif 
     407!#endif 
    417408 
    418409      ! B) ubar 
    419410      !///////// 
    420411      ! 
    421 #if defined key_dimgout 
    422       cltext='ubar amplitude and phase' 
    423       CALL dia_wri_dimg(TRIM(cdfile_name_U), TRIM(cltext), out_u, 2*nb_ana, '2') 
    424 #else 
     412!#if  
     413!#else 
    425414      DO jh = 1, nb_ana 
    426415      CALL iom_put( TRIM(tname(jh))//'x_u', out_u(:,:,jh) ) 
    427416      CALL iom_put( TRIM(tname(jh))//'y_u', out_u(:,:,nb_ana+jh) ) 
    428417      END DO 
    429 #endif 
     418!#endif 
    430419 
    431420      ! C) vbar 
    432421      !///////// 
    433422      ! 
    434 #if defined key_dimgout 
    435       cltext='vbar amplitude and phase' 
    436       CALL dia_wri_dimg(TRIM(cdfile_name_V), TRIM(cltext), out_v, 2*nb_ana, '2') 
    437 #else 
     423!#if  
     424!#else 
    438425      DO jh = 1, nb_ana 
    439426         CALL iom_put( TRIM(tname(jh))//'x_v', out_v(:,:,jh       ) ) 
    440427         CALL iom_put( TRIM(tname(jh))//'y_v', out_v(:,:,jh+nb_ana) ) 
    441428      END DO 
    442 #endif 
     429!#endif 
    443430      ! 
    444431   END SUBROUTINE dia_wri_harm 

  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/OPA_SRC/DIA/dianam.F90

    r2528 r5282  
    7272 
    7373      IF( llfsec .OR. kfreq < 0 ) THEN   ;   inbsec = kfreq                       ! output frequency already in seconds 
    74       ELSE                               ;   inbsec = kfreq * NINT( rdttra(1) )   ! from time-step to seconds 
     74      ELSE                               ;   inbsec = kfreq * NINT( rdt )   ! from time-step to seconds 
    7575      ENDIF 
    7676      iddss = NINT( rday          )                                         ! number of seconds in 1 day 
    … …  
    116116      ! date of the beginning and the end of the run 
    117117 
    118       zdrun = rdttra(1) / rday * REAL( nitend - nit000, wp )                ! length of the run in days 
    119       zjul  = fjulday - rdttra(1) / rday 
     118      zdrun = rdt / rday * REAL( nitend - nit000, wp )                ! length of the run in days 
     119      zjul  = fjulday - rdt / rday 
    120120      CALL ju2ymds( zjul        , iyear1, imonth1, iday1, zsec1 )           ! year/month/day of the beginning of run 
    121121      CALL ju2ymds( zjul + zdrun, iyear2, imonth2, iday2, zsec2 )           ! year/month/day of the end       of run 

  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/OPA_SRC/DIA/diaptr.F90

    r4990 r5282  
    585585         niter = ( nit000 - 1 ) / nn_fptr 
    586586         zdt = rdt 
    587          IF( nacc == 1 )   zdt = rdtmin 
    588587         ! 
    589588         IF(lwp) THEN 

  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/OPA_SRC/DIA/diawri.F90

    r4990 r5282  
    4343   USE lbclnk          ! ocean lateral boundary conditions (or mpp link) 
    4444   USE in_out_manager  ! I/O manager 
    45    USE diadimg         ! dimg direct access file format output 
    4645   USE iom 
    4746   USE ioipsl 
    … …  
    9796  END FUNCTION dia_wri_alloc 
    9897 
    99 #if defined key_dimgout 
    100    !!---------------------------------------------------------------------- 
    101    !!   'key_dimgout'                                      DIMG output file 
    102    !!---------------------------------------------------------------------- 
    103 #   include "diawri_dimg.h90" 
    104  
    105 #else 
    10698   !!---------------------------------------------------------------------- 
    10799   !!   Default option                                   NetCDF output file 
    108100   !!---------------------------------------------------------------------- 
    109 # if defined key_iomput 
     101#if defined key_iomput 
    110102   !!---------------------------------------------------------------------- 
    111103   !!   'key_iomput'                                        use IOM library 
    … …  
    418410      ! Define frequency of output and means 
    419411      zdt = rdt 
    420       IF( nacc == 1 ) zdt = rdtmin 
    421412      IF( ln_mskland )   THEN   ;   clop = "only(x)"   ! put 1.e+20 on land (very expensive!!) 
    422413      ELSE                      ;   clop = "x"         ! no use of the mask value (require less cpu time) 
    … …  
    762753      IF( lk_vvl ) THEN 
    763754         zw3d(:,:,:) = ( ( fse3t_n(:,:,:) - e3t_0(:,:,:) ) / e3t_0(:,:,:) * 100 * tmask(:,:,:) ) ** 2 
    764          CALL histwrite( nid_T, "vovvle3t", it, fse3t_n (:,:,:) , ndim_T , ndex_T  )   ! level thickness 
     755         CALL histwrite( nid_T, "vovvle3t", it, fse3t_n(:,:,:) , ndim_T , ndex_T  )   ! level thickness 
    765756         CALL histwrite( nid_T, "vovvldep", it, fsdept_n(:,:,:) , ndim_T , ndex_T  )   ! t-point depth 
    766757         CALL histwrite( nid_T, "vovvldef", it, zw3d             , ndim_T , ndex_T  )   ! level thickness deformation 
    … …  
    914905      ! 
    915906   END SUBROUTINE dia_wri 
    916 # endif 
    917  
    918 #endif 
     907#endif 
     908 
     909 
    919910 
    920911   SUBROUTINE dia_wri_state( cdfile_name, kt ) 
    … …  
    10351026      ! ----------------- 
    10361027      CALL histclo( id_i ) 
    1037 #if ! defined key_iomput && ! defined key_dimgout 
     1028#if ! defined key_iomput 
    10381029      IF( ninist /= 1  ) THEN 
    10391030         CALL histclo( nid_T ) 
Note: See TracChangeset for help on using the changeset viewer.