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 6032 for branches/2014/dev_r4650_UKMO10_Tidally_Meaned_Diagnostics/NEMOGCM/NEMO/OPA_SRC/DIA/dia25h.F90 – NEMO

Ignore:
Timestamp:
2015-12-11T09:24:14+01:00 (8 years ago)
Author:
timgraham
Message:

Remove insitu temperature for now to be replaced with a versios consistent with EOS later

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2014/dev_r4650_UKMO10_Tidally_Meaned_Diagnostics/NEMOGCM/NEMO/OPA_SRC/DIA/dia25h.F90

    r5989 r6032  
    88   USE oce             ! ocean dynamics and tracers variables 
    99   USE dom_oce         ! ocean space and time domain 
    10    USE diainsitutem, ONLY: rinsitu_t, theta2t 
    1110   USE in_out_manager  ! I/O units 
    1211   USE iom             ! I/0 library 
     
    2928 
    3029  !! * variables for calculating 25-hourly means 
    31    REAL(wp),SAVE, ALLOCATABLE,   DIMENSION(:,:,:) ::   tn_25h  , sn_25h, rinsitu_t_25h   
     30   REAL(wp),SAVE, ALLOCATABLE,   DIMENSION(:,:,:) ::   tn_25h  , sn_25h 
    3231   REAL(wp),SAVE, ALLOCATABLE,   DIMENSION(:,:)   ::   sshn_25h  
    3332   REAL(wp),SAVE, ALLOCATABLE,   DIMENSION(:,:,:) ::   un_25h  , vn_25h  , wn_25h 
     
    9594         CALL ctl_stop( 'dia_25h: unable to allocate sn_25h' )   ;   RETURN 
    9695      ENDIF 
    97       ALLOCATE( rinsitu_t_25h(jpi,jpj,jpk), STAT=ierror ) 
    98       IF( ierror > 0 ) THEN 
    99          CALL ctl_stop( 'dia_25h: unable to allocate rinsitu_t_25h' )   ;   RETURN 
    100       ENDIF 
    10196      ALLOCATE( un_25h(jpi,jpj,jpk), STAT=ierror ) 
    10297      IF( ierror > 0 ) THEN 
     
    142137      sn_25h(:,:,:) = tsb(:,:,:,jp_sal) 
    143138      CALL theta2t 
    144       rinsitu_t_25h(:,:,:) = rinsitu_t(:,:,:) 
    145139      sshn_25h(:,:) = sshb(:,:) 
    146140      un_25h(:,:,:) = ub(:,:,:) 
     
    235229         tn_25h(:,:,:)        = tn_25h(:,:,:) + tsn(:,:,:,jp_tem) 
    236230         sn_25h(:,:,:)        = sn_25h(:,:,:) + tsn(:,:,:,jp_sal) 
    237          CALL theta2t 
    238          rinsitu_t_25h(:,:,:)  = rinsitu_t_25h(:,:,:) + rinsitu_t(:,:,:) 
    239231         sshn_25h(:,:)        = sshn_25h(:,:) + sshn (:,:) 
    240232         un_25h(:,:,:)        = un_25h(:,:,:) + un(:,:,:) 
     
    267259            tn_25h(:,:,:)        = tn_25h(:,:,:) / 25.0_wp 
    268260            sn_25h(:,:,:)        = sn_25h(:,:,:) / 25.0_wp 
    269             rinsitu_t_25h(:,:,:)  = rinsitu_t_25h(:,:,:) / 25.0_wp 
    270261            sshn_25h(:,:)        = sshn_25h(:,:) / 25.0_wp 
    271262            un_25h(:,:,:)        = un_25h(:,:,:) / 25.0_wp 
     
    286277            zw3d(:,:,:) = tn_25h(:,:,:)*tmask(:,:,:) + zmdi*(1.0-tmask(:,:,:)) 
    287278            CALL iom_put("temper25h", zw3d)   ! potential temperature 
    288             CALL theta2t                                                                    ! calculate insitu temp 
    289             zw3d(:,:,:) = rinsitu_t_25h(:,:,:)*tmask(:,:,:) + zmdi*(1.0-tmask(:,:,:)) 
    290             CALL iom_put("tempis25h", zw3d)   ! in-situ temperature 
    291279            zw3d(:,:,:) = sn_25h(:,:,:)*tmask(:,:,:) + zmdi*(1.0-tmask(:,:,:)) 
    292280            CALL iom_put( "salin25h", zw3d  )   ! salinity 
     
    319307            tn_25h(:,:,:) = tsn(:,:,:,jp_tem) 
    320308            sn_25h(:,:,:) = tsn(:,:,:,jp_sal) 
    321             CALL theta2t 
    322             rinsitu_t_25h(:,:,:) = rinsitu_t(:,:,:) 
    323309            sshn_25h(:,:) = sshn (:,:) 
    324310            un_25h(:,:,:) = un(:,:,:) 
Note: See TracChangeset for help on using the changeset viewer.