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 1472 for trunk/NEMO – NEMO

Changeset 1472 for trunk/NEMO


Ignore:
Timestamp:
2009-06-11T16:44:14+02:00 (15 years ago)
Author:
smasson
Message:

small bugfix in the coupling interface when sending surface current, see ticket:455

File:
1 edited

Legend:

Unmodified
Added
Removed

  • trunk/NEMO/OPA_SRC/SBC/sbccpl.F90

    r1469 r1472  
    10941094               DO jj = 2, jpjm1 
    10951095                  DO ji = fs_2, fs_jpim1   ! vector opt. 
    1096                      zotx1(ji,jj) = 0.5 * ( un      (ji,jj,1) + un      (ji-1,jj  ,1) ) * zfr_l(ji,jj)   
    1097                      zoty1(ji,jj) = 0.5 * ( vn      (ji,jj,1) + vn      (ji  ,jj-1,1) ) * zfr_l(ji,jj) 
    1098                      zitx1(ji,jj) = 0.5 * ( utau_ice(ji,jj  ) + utau_ice(ji-1,jj    ) ) *  fr_i(ji,jj) 
    1099                      zity1(ji,jj) = 0.5 * ( vtau_ice(ji,jj  ) + vtau_ice(ji  ,jj-1  ) ) *  fr_i(ji,jj) 
     1096                     zotx1(ji,jj) = 0.5 * ( un   (ji,jj,1) + un   (ji-1,jj  ,1) ) * zfr_l(ji,jj)   
     1097                     zoty1(ji,jj) = 0.5 * ( vn   (ji,jj,1) + vn   (ji  ,jj-1,1) ) * zfr_l(ji,jj) 
     1098                     zitx1(ji,jj) = 0.5 * ( u_ice(ji,jj  ) + u_ice(ji-1,jj    ) ) *  fr_i(ji,jj) 
     1099                     zity1(ji,jj) = 0.5 * ( v_ice(ji,jj  ) + v_ice(ji  ,jj-1  ) ) *  fr_i(ji,jj) 
    11001100                  END DO 
    11011101               END DO 
    … …  
    11031103               DO jj = 2, jpjm1 
    11041104                  DO ji = fs_2, fs_jpim1   ! vector opt. 
    1105                      zotx1(ji,jj) = 0.5  * ( un(ji,jj,1)         + un(ji-1,jj-1,1)  ) * zfr_l(ji,jj)   
    1106                      zoty1(ji,jj) = 0.5  * ( vn(ji,jj,1)         + vn(ji  ,jj-1,1)  ) * zfr_l(ji,jj)   
    1107                      zitx1(ji,jj) = 0.25 * ( utau_ice(ji+1,jj+1) + utau_ice(ji,jj+1)   & 
    1108                         &                  + utau_ice(ji+1,jj  ) + utau_ice(ji,jj  ) ) * fr_i(ji,jj) 
    1109                      zity1(ji,jj) = 0.25 * ( vtau_ice(ji+1,jj+1) + vtau_ice(ji,jj+1)   & 
    1110                         &                  + vtau_ice(ji+1,jj  ) + vtau_ice(ji,jj  ) ) * fr_i(ji,jj) 
     1105                     zotx1(ji,jj) = 0.5  * ( un(ji,jj,1)      + un(ji-1,jj-1,1) ) * zfr_l(ji,jj)   
     1106                     zoty1(ji,jj) = 0.5  * ( vn(ji,jj,1)      + vn(ji  ,jj-1,1) ) * zfr_l(ji,jj)   
     1107                     zitx1(ji,jj) = 0.25 * ( u_ice(ji+1,jj+1) + u_ice(ji,jj+1)                     & 
     1108                        &                  + u_ice(ji+1,jj  ) + u_ice(ji,jj  )  ) * fr_i(ji,jj) 
     1109                     zity1(ji,jj) = 0.25 * ( v_ice(ji+1,jj+1) + v_ice(ji,jj+1)                     & 
     1110                        &                  + v_ice(ji+1,jj  ) + v_ice(ji,jj  )  ) * fr_i(ji,jj) 
    11111111                  END DO 
    11121112               END DO 
    … …  
    11141114               DO jj = 2, jpjm1 
    11151115                  DO ji = fs_2, fs_jpim1   ! vector opt. 
    1116                      zotx1(ji,jj) = 0.5  * ( un(ji,jj,1)         + un(ji-1,jj-1,1)  ) * zfr_l(ji,jj)   
    1117                      zoty1(ji,jj) = 0.5  * ( vn(ji,jj,1)         + vn(ji  ,jj-1,1)  ) * zfr_l(ji,jj)   
    1118                      zitx1(ji,jj) = 0.25 * ( utau_ice(ji-1,jj-1) + utau_ice(ji,jj-1)   & 
    1119                         &                  + utau_ice(ji-1,jj  ) + utau_ice(ji,jj  ) ) * fr_i(ji,jj) 
    1120                      zity1(ji,jj) = 0.25 * ( vtau_ice(ji-1,jj-1) + vtau_ice(ji,jj-1)   & 
    1121                         &                  + vtau_ice(ji-1,jj  ) + vtau_ice(ji,jj  ) ) * fr_i(ji,jj) 
     1116                     zotx1(ji,jj) = 0.5  * ( un(ji,jj,1)      + un(ji-1,jj-1,1) ) * zfr_l(ji,jj)   
     1117                     zoty1(ji,jj) = 0.5  * ( vn(ji,jj,1)      + vn(ji  ,jj-1,1) ) * zfr_l(ji,jj)   
     1118                     zitx1(ji,jj) = 0.25 * ( u_ice(ji-1,jj-1) + u_ice(ji,jj-1)                     & 
     1119                        &                  + u_ice(ji-1,jj  ) + u_ice(ji,jj  )  ) * fr_i(ji,jj) 
     1120                     zity1(ji,jj) = 0.25 * ( v_ice(ji-1,jj-1) + v_ice(ji,jj-1)                     & 
     1121                        &                  + v_ice(ji-1,jj  ) + v_ice(ji,jj  )  ) * fr_i(ji,jj) 
    11221122                  END DO 
    11231123               END DO 
    … …  
    11291129               DO jj = 2, jpjm1 
    11301130                  DO ji = fs_2, fs_jpim1   ! vector opt. 
    1131                      zotx1(ji,jj) = 0.5 * ( un      (ji,jj,1) + un      (ji-1,jj  ,1) ) * zfr_l(ji,jj) & 
    1132                         &         + 0.5 * ( utau_ice(ji,jj  ) + utau_ice(ji-1,jj    ) ) *  fr_i(ji,jj) 
    1133                      zoty1(ji,jj) = 0.5 * ( vn      (ji,jj,1) + vn      (ji  ,jj-1,1) ) * zfr_l(ji,jj) & 
    1134                         &         + 0.5 * ( vtau_ice(ji,jj  ) + vtau_ice(ji  ,jj-1  ) ) *  fr_i(ji,jj) 
     1131                     zotx1(ji,jj) = 0.5 * ( un   (ji,jj,1) + un   (ji-1,jj  ,1) ) * zfr_l(ji,jj)  & 
     1132                        &         + 0.5 * ( u_ice(ji,jj  ) + u_ice(ji-1,jj    ) ) *  fr_i(ji,jj) 
     1133                     zoty1(ji,jj) = 0.5 * ( vn   (ji,jj,1) + vn   (ji  ,jj-1,1) ) * zfr_l(ji,jj)  & 
     1134                        &         + 0.5 * ( v_ice(ji,jj  ) + v_ice(ji  ,jj-1  ) ) *  fr_i(ji,jj) 
    11351135                  END DO 
    11361136               END DO 
    … …  
    11381138               DO jj = 2, jpjm1 
    11391139                  DO ji = fs_2, fs_jpim1   ! vector opt. 
    1140                      zotx1(ji,jj) = 0.5  * ( un(ji,jj,1)         + un(ji-1,jj-1,1)   ) * zfr_l(ji,jj) &    
    1141                         &         + 0.25 * ( utau_ice(ji+1,jj+1) + utau_ice(ji,jj+1)   & 
    1142                         &                  + utau_ice(ji+1,jj  ) + utau_ice(ji,jj  ) ) *  fr_i(ji,jj) 
    1143                      zoty1(ji,jj) = 0.5  * ( vn(ji,jj,1)         + vn(ji  ,jj-1,1)   ) * zfr_l(ji,jj) &  
    1144                         &         + 0.25 * ( vtau_ice(ji+1,jj+1) + vtau_ice(ji,jj+1)   & 
    1145                         &                  + vtau_ice(ji+1,jj  ) + vtau_ice(ji,jj  ) ) *  fr_i(ji,jj) 
     1140                     zotx1(ji,jj) = 0.5  * ( un(ji,jj,1)      + un(ji-1,jj-1,1) ) * zfr_l(ji,jj)  &    
     1141                        &         + 0.25 * ( u_ice(ji+1,jj+1) + u_ice(ji,jj+1)                     & 
     1142                        &                  + u_ice(ji+1,jj  ) + u_ice(ji,jj  ) ) *  fr_i(ji,jj) 
     1143                     zoty1(ji,jj) = 0.5  * ( vn(ji,jj,1)      + vn(ji  ,jj-1,1) ) * zfr_l(ji,jj)  &  
     1144                        &         + 0.25 * ( v_ice(ji+1,jj+1) + v_ice(ji,jj+1)                     & 
     1145                        &                  + v_ice(ji+1,jj  ) + v_ice(ji,jj  ) ) *  fr_i(ji,jj) 
    11461146                  END DO 
    11471147               END DO 
    … …  
    11491149               DO jj = 2, jpjm1 
    11501150                  DO ji = fs_2, fs_jpim1   ! vector opt. 
    1151                      zotx1(ji,jj) = 0.5  * ( un(ji,jj,1)         + un(ji-1,jj-1,1)   ) * zfr_l(ji,jj) &    
    1152                         &         + 0.25 * ( utau_ice(ji-1,jj-1) + utau_ice(ji,jj-1)   & 
    1153                         &                  + utau_ice(ji-1,jj  ) + utau_ice(ji,jj  ) ) *  fr_i(ji,jj) 
    1154                      zoty1(ji,jj) = 0.5  * ( vn(ji,jj,1)         + vn(ji  ,jj-1,1)   ) * zfr_l(ji,jj) &  
    1155                         &         + 0.25 * ( vtau_ice(ji-1,jj-1) + vtau_ice(ji,jj-1)   & 
    1156                         &                  + vtau_ice(ji-1,jj  ) + vtau_ice(ji,jj  ) ) *  fr_i(ji,jj) 
     1151                     zotx1(ji,jj) = 0.5  * ( un(ji,jj,1)      + un(ji-1,jj-1,1) ) * zfr_l(ji,jj)  &    
     1152                        &         + 0.25 * ( u_ice(ji-1,jj-1) + u_ice(ji,jj-1)                     & 
     1153                        &                  + u_ice(ji-1,jj  ) + u_ice(ji,jj  ) ) *  fr_i(ji,jj) 
     1154                     zoty1(ji,jj) = 0.5  * ( vn(ji,jj,1)      + vn(ji  ,jj-1,1) ) * zfr_l(ji,jj)  &  
     1155                        &         + 0.25 * ( v_ice(ji-1,jj-1) + v_ice(ji,jj-1)                     & 
     1156                        &                  + v_ice(ji-1,jj  ) + v_ice(ji,jj  ) ) *  fr_i(ji,jj) 
    11571157                  END DO 
    11581158               END DO 
Note: See TracChangeset for help on using the changeset viewer.