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

Changeset 6074


Ignore:
Timestamp:
2015-12-16T17:42:24+01:00 (8 years ago)
Author:
mathiot
Message:

change variable names

Location:
branches/2015/dev_merge_2015/NEMOGCM/NEMO/OPA_SRC
Files:
4 edited

Legend:

Unmodified
Added
Removed

  • branches/2015/dev_merge_2015/NEMOGCM/NEMO/OPA_SRC/DOM/iscplhsb.F90

    r6069 r6074  
    103103 
    104104                  ! volume differences 
    105                   zde3t = fse3t_n(ji,jj,jk) * tmask(ji,jj,jk) - pe3t_b(ji,jj,jk) * ptmask_b(ji,jj,jk) 
     105                  zde3t = e3t_n(ji,jj,jk) * tmask(ji,jj,jk) - pe3t_b(ji,jj,jk) * ptmask_b(ji,jj,jk) 
    106106 
    107107                  ! heat diff 
    108                   zdtem = tsn(ji,jj,jk,jp_tem) * fse3t_n(ji,jj,jk) *  tmask  (ji,jj,jk)   & 
     108                  zdtem = tsn(ji,jj,jk,jp_tem) * e3t_n(ji,jj,jk) *  tmask  (ji,jj,jk)   & 
    109109                        - tsb(ji,jj,jk,jp_tem) * pe3t_b (ji,jj,jk) * ptmask_b(ji,jj,jk) 
    110110                  ! salt diff 
    111                   zdsal = tsn(ji,jj,jk,jp_sal) * fse3t_n(ji,jj,jk) *  tmask  (ji,jj,jk)   & 
     111                  zdsal = tsn(ji,jj,jk,jp_sal) * e3t_n(ji,jj,jk) *  tmask  (ji,jj,jk)   & 
    112112                        - tsb(ji,jj,jk,jp_sal) * pe3t_b (ji,jj,jk) * ptmask_b(ji,jj,jk) 
    113113                
    … …  
    313313         DO jj = 1, jpj 
    314314            DO ji = 1, jpi 
    315                phdivn(ji,jj,jk) = phdivn(ji,jj,jk) + hdiv_iscpl(ji,jj,jk) / fse3t_n(ji,jj,jk) 
     315               phdivn(ji,jj,jk) = phdivn(ji,jj,jk) + hdiv_iscpl(ji,jj,jk) / e3t_n(ji,jj,jk) 
    316316            END DO 
    317317         END DO 

  • branches/2015/dev_merge_2015/NEMOGCM/NEMO/OPA_SRC/DOM/iscplrst.F90

    r6069 r6074  
    6969      CALL iom_get( numror, jpdom_autoglo, 'vmask'  , zvmask_b   ) ! need to correct barotropic velocity 
    7070      CALL iom_get( numror, jpdom_autoglo, 'smask'  , zsmask_b   ) ! need to correct barotropic velocity 
    71       CALL iom_get( numror, jpdom_autoglo, 'fse3t_n' , ze3t_b(:,:,:) )  ! need to compute temperature correction 
    72       CALL iom_get( numror, jpdom_autoglo, 'fse3u_n' , ze3u_b(:,:,:) )  ! need to correct barotropic velocity 
    73       CALL iom_get( numror, jpdom_autoglo, 'fse3v_n' , ze3v_b(:,:,:) )  ! need to correct barotropic velocity 
    74       CALL iom_get( numror, jpdom_autoglo, 'fsdepw_n', zdepw_b(:,:,:) ) ! need to interpol vertical profile (vvl) 
     71      CALL iom_get( numror, jpdom_autoglo, 'e3t_n' , ze3t_b(:,:,:) )  ! need to compute temperature correction 
     72      CALL iom_get( numror, jpdom_autoglo, 'e3u_n' , ze3u_b(:,:,:) )  ! need to correct barotropic velocity 
     73      CALL iom_get( numror, jpdom_autoglo, 'e3v_n' , ze3v_b(:,:,:) )  ! need to correct barotropic velocity 
     74      CALL iom_get( numror, jpdom_autoglo, 'gdepw_n', zdepw_b(:,:,:) ) ! need to interpol vertical profile (vvl) 
    7575 
    7676      !! read namelist 
    … …  
    114114 
    115115      !! set _b and _n vertical scale factor equal 
    116       fse3t_b (:,:,:) = fse3t_n (:,:,:) 
    117       fse3u_b (:,:,:) = fse3u_n (:,:,:) 
    118       fse3v_b (:,:,:) = fse3v_n (:,:,:) 
     116      e3t_b (:,:,:) = e3t_n (:,:,:) 
     117      e3u_b (:,:,:) = e3u_n (:,:,:) 
     118      e3v_b (:,:,:) = e3v_n (:,:,:) 
    119119 
    120120      IF ( lk_vvl ) THEN 
    121          fse3uw_b(:,:,:) = fse3uw_n(:,:,:) 
    122          fse3vw_b(:,:,:) = fse3vw_n(:,:,:) 
    123          fsdept_b(:,:,:) = fsdept_n(:,:,:) 
    124          fsdepw_b(:,:,:) = fsdepw_n(:,:,:) 
    125          hu_b (:,:) = hu(:,:) 
    126          hv_b (:,:) = hv(:,:) 
    127          hur_b(:,:) = hur(:,:) 
    128          hvr_b(:,:) = hvr(:,:) 
     121         e3uw_b(:,:,:)  = e3uw_n(:,:,:) 
     122         e3vw_b(:,:,:)  = e3vw_n(:,:,:) 
     123         dept_b(:,:,:)  = gdept_n(:,:,:) 
     124         gdepw_b(:,:,:) = gdepw_n(:,:,:) 
     125         hu_b (:,:) = hu_n(:,:) 
     126         hv_b (:,:) = hv_n(:,:) 
     127         r1_hu_b(:,:) = r1_hu_n(:,:) 
     128         r1_hv_b(:,:) = r1_hv_n(:,:) 
    129129      END IF 
    130130      ! 
    … …  
    220220               DO ji=1,jpi 
    221221                  IF (tmask(ji,jj,1) == 0._wp .OR. ptmask_b(ji,jj,1) == 0._wp) THEN 
    222                      fse3t_n(ji,jj,jk) = e3t_0(ji,jj,jk) * ( 1._wp + sshn(ji,jj) / ( ht_0(ji,jj) + 1._wp - ssmask(ji,jj) ) * tmask(ji,jj,jk) ) 
     222                     e3t_n(ji,jj,jk) = e3t_0(ji,jj,jk) * ( 1._wp + sshn(ji,jj) / ( ht_0(ji,jj) + 1._wp - ssmask(ji,jj) ) * tmask(ji,jj,jk) ) 
    223223                  ENDIF 
    224224               END DO 
    … …  
    226226         END DO 
    227227 
    228          CALL dom_vvl_interpol( fse3t_n(:,:,:), fse3u_n(:,:,:), 'U' ) 
    229          CALL dom_vvl_interpol( fse3t_n(:,:,:), fse3v_n(:,:,:), 'V' ) 
    230          CALL dom_vvl_interpol( fse3u_n(:,:,:), fse3f_n(:,:,:), 'F' ) 
     228         CALL dom_vvl_interpol( e3t_n(:,:,:), e3u_n(:,:,:), 'U' ) 
     229         CALL dom_vvl_interpol( e3t_n(:,:,:), e3v_n(:,:,:), 'V' ) 
     230         CALL dom_vvl_interpol( e3u_n(:,:,:), e3f_n(:,:,:), 'F' ) 
    231231 
    232232      ! Vertical scale factor interpolations 
    233233      ! ------------------------------------ 
    234          CALL dom_vvl_interpol( fse3t_n(:,:,:), fse3w_n (:,:,:), 'W'  ) 
    235          CALL dom_vvl_interpol( fse3u_n(:,:,:), fse3uw_n(:,:,:), 'UW' ) 
    236          CALL dom_vvl_interpol( fse3v_n(:,:,:), fse3vw_n(:,:,:), 'VW' ) 
     234         CALL dom_vvl_interpol( e3t_n(:,:,:), e3w_n (:,:,:), 'W'  ) 
     235         CALL dom_vvl_interpol( e3u_n(:,:,:), e3uw_n(:,:,:), 'UW' ) 
     236         CALL dom_vvl_interpol( e3v_n(:,:,:), e3vw_n(:,:,:), 'VW' ) 
    237237 
    238238      ! t- and w- points depth 
    239239      ! ---------------------- 
    240          fsdept_n(:,:,1) = 0.5_wp * fse3w_n(:,:,1) 
    241          fsdepw_n(:,:,1) = 0.0_wp 
    242          fsde3w_n(:,:,1) = fsdept_n(:,:,1) - sshn(:,:) 
     240         gdept_n(:,:,1) = 0.5_wp * e3w_n(:,:,1) 
     241         gdepw_n(:,:,1) = 0.0_wp 
     242         gde3w_n(:,:,1) = gdept_n(:,:,1) - sshn(:,:) 
    243243         DO jj = 1,jpj 
    244244            DO ji = 1,jpi 
    245245               DO jk = 2,mikt(ji,jj)-1 
    246                   fsdept_n(ji,jj,jk) = gdept_0(ji,jj,jk) 
    247                   fsdepw_n(ji,jj,jk) = gdepw_0(ji,jj,jk) 
    248                   fsde3w_n(ji,jj,jk) = gdept_0(ji,jj,jk) - sshn(ji,jj) 
     246                  gdept_n(ji,jj,jk) = gdept_0(ji,jj,jk) 
     247                  gdepw_n(ji,jj,jk) = gdepw_0(ji,jj,jk) 
     248                  gde3w_n(ji,jj,jk) = gdept_0(ji,jj,jk) - sshn(ji,jj) 
    249249               END DO 
    250250               IF (mikt(ji,jj) > 1) THEN 
    251251                  jk = mikt(ji,jj) 
    252                   fsdept_n(ji,jj,jk) = gdepw_0(ji,jj,jk) + 0.5_wp * fse3w_n(ji,jj,jk) 
    253                   fsdepw_n(ji,jj,jk) = gdepw_0(ji,jj,jk) 
    254                   fsde3w_n(ji,jj,jk) = fsdept_n(ji,jj,jk  ) - sshn   (ji,jj) 
     252                  gdept_n(ji,jj,jk) = gdepw_0(ji,jj,jk) + 0.5_wp * e3w_n(ji,jj,jk) 
     253                  gdepw_n(ji,jj,jk) = gdepw_0(ji,jj,jk) 
     254                  gde3w_n(ji,jj,jk) = gdept_n(ji,jj,jk  ) - sshn   (ji,jj) 
    255255               END IF 
    256256               DO jk = mikt(ji,jj)+1, jpk 
    257                   fsdept_n(ji,jj,jk) = fsdept_n(ji,jj,jk-1) + fse3w_n(ji,jj,jk) 
    258                   fsdepw_n(ji,jj,jk) = fsdepw_n(ji,jj,jk-1) + fse3t_n(ji,jj,jk-1) 
    259                   fsde3w_n(ji,jj,jk) = fsdept_n(ji,jj,jk  ) - sshn  (ji,jj) 
     257                  gdept_n(ji,jj,jk) = gdept_n(ji,jj,jk-1) + e3w_n(ji,jj,jk) 
     258                  gdepw_n(ji,jj,jk) = gdepw_n(ji,jj,jk-1) + e3t_n(ji,jj,jk-1) 
     259                  gde3w_n(ji,jj,jk) = gdept_n(ji,jj,jk  ) - sshn (ji,jj) 
    260260               END DO 
    261261            END DO 
    … …  
    266266         ht(:,:) = 0._wp ; hu(:,:) = 0._wp ; hv(:,:) = 0._wp 
    267267         DO jk = 1, jpkm1 
    268             hu(:,:) = hu(:,:) + fse3u_n(:,:,jk) * umask(:,:,jk) 
    269             hv(:,:) = hv(:,:) + fse3v_n(:,:,jk) * vmask(:,:,jk) 
    270             ht(:,:) = ht(:,:) + fse3t_n(:,:,jk) * tmask(:,:,jk) 
     268            hu_n(:,:) = hu_n(:,:) + e3u_n(:,:,jk) * umask(:,:,jk) 
     269            hv_n(:,:) = hv_n(:,:) + e3v_n(:,:,jk) * vmask(:,:,jk) 
     270            ht_n(:,:) = ht_n(:,:) + e3t_n(:,:,jk) * tmask(:,:,jk) 
    271271         END DO 
    272272         !                                        ! Inverse of the local depth 
    273          hur(:,:) = 1._wp / ( hu(:,:) + 1._wp - ssumask(:,:) ) * ssumask(:,:) 
    274          hvr(:,:) = 1._wp / ( hv(:,:) + 1._wp - ssvmask(:,:) ) * ssvmask(:,:) 
     273         r1_hu_n(:,:) = 1._wp / ( hu_n(:,:) + 1._wp - ssumask(:,:) ) * ssumask(:,:) 
     274         r1_hv_n(:,:) = 1._wp / ( hv_n(:,:) + 1._wp - ssvmask(:,:) ) * ssvmask(:,:) 
    275275 
    276276      END IF 
    … …  
    284284         DO jj = 1,jpj 
    285285            DO ji = 1,jpi 
    286                un(ji,jj,jk) = ub(ji,jj,jk)*pe3u_b(ji,jj,jk)*pumask_b(ji,jj,jk)/fse3u_n(ji,jj,jk)*umask(ji,jj,jk); 
    287                vn(ji,jj,jk) = vb(ji,jj,jk)*pe3v_b(ji,jj,jk)*pvmask_b(ji,jj,jk)/fse3v_n(ji,jj,jk)*vmask(ji,jj,jk); 
     286               un(ji,jj,jk) = ub(ji,jj,jk)*pe3u_b(ji,jj,jk)*pumask_b(ji,jj,jk)/e3u_n(ji,jj,jk)*umask(ji,jj,jk); 
     287               vn(ji,jj,jk) = vb(ji,jj,jk)*pe3v_b(ji,jj,jk)*pvmask_b(ji,jj,jk)/e3v_n(ji,jj,jk)*vmask(ji,jj,jk); 
    288288            END DO 
    289289         END DO 
    … …  
    296296      ztrp(:,:,:) = vb(:,:,:)*pe3v_b(:,:,:);  
    297297      zbvb(:,:)   = SUM(ztrp,DIM=3) 
    298       ztrp(:,:,:) = un(:,:,:)*fse3u_n(:,:,:);  
     298      ztrp(:,:,:) = un(:,:,:)*e3u_n(:,:,:);  
    299299      zbun(:,:)   = SUM(ztrp,DIM=3) 
    300       ztrp(:,:,:) = vn(:,:,:)*fse3v_n(:,:,:);  
     300      ztrp(:,:,:) = vn(:,:,:)*e3v_n(:,:,:);  
    301301      zbvn(:,:)   = SUM(ztrp,DIM=3) 
    302302 
    … …  
    305305      zhv1=0.0_wp ; 
    306306      DO jk  = 1,jpk 
    307         zhu1(:,:) = zhu1(:,:) + fse3u_n(:,:,jk) * umask(:,:,jk) 
    308         zhv1(:,:) = zhv1(:,:) + fse3v_n(:,:,jk) * vmask(:,:,jk) 
     307        zhu1(:,:) = zhu1(:,:) + e3u_n(:,:,jk) * umask(:,:,jk) 
     308        zhv1(:,:) = zhv1(:,:) + e3v_n(:,:,jk) * vmask(:,:,jk) 
    309309      END DO 
    310310       
    … …  
    392392                  IF (zwmaskn(ji,jj,jk) * zwmaskb(ji,jj,jk) == 1._wp .AND. (tmask(ji,jj,1)==0._wp .OR. ptmask_b(ji,jj,1)==0._wp) ) THEN 
    393393                     !compute weight 
    394                      zdzp1 = MAX(0._wp,fsdepw_n(ji,jj,jk+1) - pdepw_b(ji,jj,jk+1)) 
    395                      zdz   =           fsdepw_n(ji,jj,jk+1) - pdepw_b(ji,jj,jk  )  
    396                      zdzm1 = MAX(0._wp,pdepw_b(ji,jj,jk  )  - fsdepw_n(ji,jj,jk  )) 
     394                     zdzp1 = MAX(0._wp,gdepw_n(ji,jj,jk+1) - pdepw_b(ji,jj,jk+1)) 
     395                     zdz   =           gdepw_n(ji,jj,jk+1) - pdepw_b(ji,jj,jk  )  
     396                     zdzm1 = MAX(0._wp,pdepw_b(ji,jj,jk  ) - gdepw_n(ji,jj,jk  )) 
    397397                     IF (zdz .LT. 0._wp) THEN  
    398398                        CALL ctl_stop( 'STOP', 'rst_iscpl : unable to compute the interpolation' ) 
    … …  
    400400                     tsn(ji,jj,jk,jp_tem) = ( zdzp1*tsb(ji,jj,jk+1,jp_tem)      & 
    401401                        &                   + zdz  *tsb(ji,jj,jk  ,jp_tem)      & 
    402                         &                   + zdzm1*tsb(ji,jj,jk-1,jp_tem)      )/fse3t_n(ji,jj,jk) 
     402                        &                   + zdzm1*tsb(ji,jj,jk-1,jp_tem)      )/e3t_n(ji,jj,jk) 
    403403                     tsn(ji,jj,jk,jp_sal) = ( zdzp1*tsb(ji,jj,jk+1,jp_sal)      & 
    404404                        &                   + zdz  *tsb(ji,jj,jk  ,jp_sal)      & 
    405                         &                   + zdzm1*tsb(ji,jj,jk-1,jp_sal)      )/fse3t_n(ji,jj,jk) 
     405                        &                   + zdzm1*tsb(ji,jj,jk-1,jp_sal)      )/e3t_n(ji,jj,jk) 
    406406                  END IF 
    407407               END DO 

  • branches/2015/dev_merge_2015/NEMOGCM/NEMO/OPA_SRC/DYN/dynhpg.F90

    r6069 r6074  
    215215         ! compute density of the water displaced by the ice shelf  
    216216         DO jk = 1, jpk 
    217             CALL eos(ztstop(:,:,:),fsdept_n(:,:,jk),zrhd(:,:,jk)) 
     217            CALL eos(ztstop(:,:,:),gdept_n(:,:,jk),zrhd(:,:,jk)) 
    218218         END DO 
    219219       
    … …  
    228228            DO ji = 1, jpi 
    229229               ikt=mikt(ji,jj) 
    230                ziceload(ji,jj) = ziceload(ji,jj) + (znad + zrhd(ji,jj,1) ) * fse3w(ji,jj,1) * (1._wp - tmask(ji,jj,1)) 
     230               ziceload(ji,jj) = ziceload(ji,jj) + (znad + zrhd(ji,jj,1) ) * e3w(ji,jj,1) * (1._wp - tmask(ji,jj,1)) 
    231231               DO jk=2,ikt-1 
    232                   ziceload(ji,jj) = ziceload(ji,jj) + (2._wp * znad + zrhd(ji,jj,jk-1) + zrhd(ji,jj,jk)) * fse3w(ji,jj,jk) & 
     232                  ziceload(ji,jj) = ziceload(ji,jj) + (2._wp * znad + zrhd(ji,jj,jk-1) + zrhd(ji,jj,jk)) * e3w(ji,jj,jk) & 
    233233                     &                              * (1._wp - tmask(ji,jj,jk)) 
    234234               END DO 

  • branches/2015/dev_merge_2015/NEMOGCM/NEMO/OPA_SRC/TRA/trasbc.F90

    r6069 r6074  
    216216            DO jj = 2, jpj  
    217217               DO ji = fs_2, fs_jpim1 
    218                   zdep = 1._wp / fse3t_n(ji,jj,jk)  
     218                  zdep = 1._wp / e3t_n(ji,jj,jk)  
    219219                  tsa(ji,jj,jk,jp_tem) = tsa(ji,jj,jk,jp_tem) - htsc_iscpl(ji,jj,jk,jp_tem)                       & 
    220220                      &                                         * zdep 
Note: See TracChangeset for help on using the changeset viewer.