Changeset 11605
- Timestamp:
- 2019-09-26T21:23:26+02:00 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/dev_r11233_AGRIF-05_jchanut_vert_coord_interp/src/NST/vremap.F90
r11603 r11605 73 73 74 74 DO jk = 2, kjpk_in 75 zwork(jk,1)=1. /(phin(jk-1)+phin(jk))75 zwork(jk,1)=1._wp/(phin(jk-1)+phin(jk)) 76 76 END DO 77 77 … … 93 93 q001 = phin(jk)*zwork2(jk+1,1) 94 94 q002 = phin(jk)*zwork2(jk,1) 95 IF (q001*q002 < 0 ) then95 IF (q001*q002 < 0._wp) then 96 96 q001 = 0._wp 97 97 q002 = 0._wp … … 129 129 END DO 130 130 131 zwork2(1,3) = 2 *coeffremap(1,1)-zwork2(2,3)132 zwork2(kjpk_in+1,3)=2 *coeffremap(kjpk_in,1)-zwork2(kjpk_in,3)131 zwork2(1,3) = 2._wp*coeffremap(1,1)-zwork2(2,3) 132 zwork2(kjpk_in+1,3)=2._wp*coeffremap(kjpk_in,1)-zwork2(kjpk_in,3) 133 133 134 134 DO jk = 1, kjpk_in 135 135 q01=zwork2(jk+1,3)-coeffremap(jk,1) 136 136 q02=coeffremap(jk,1)-zwork2(jk,3) 137 q001=2. *q01138 q002=2. *q02139 IF (q01*q02<0 ) then137 q001=2._wp*q01 138 q002=2._wp*q02 139 IF (q01*q02<0._wp) then 140 140 q01=0._wp 141 141 q02=0._wp … … 190 190 q001 = (zbox-z_win(kbox))/phin(kbox) 191 191 q002 = (ztop-z_win(kbox))/phin(kbox) 192 q01=q001**2+q002**2+q001*q002+1. -2.*(q001+q002)193 q02=q01-1. +(q001+q002)192 q01=q001**2+q002**2+q001*q002+1._wp-2._wp*(q001+q002) 193 q02=q01-1._wp+(q001+q002) 194 194 q0=1._wp-q01-q02 195 195 ELSE
Note: See TracChangeset
for help on using the changeset viewer.