94 |
end do |
end do |
95 |
end do |
end do |
96 |
end do |
end do |
97 |
elseif ( buoyancyRelation .eq. 'ATMOSPHERIC' |
elseif ( buoyancyRelation .eq. 'OCEANICP' ) then |
|
& .or. buoyancyRelation .eq. 'OCEANICP' ) then |
|
98 |
C |
C |
99 |
C in pressure coordinates the pressure is just the coordinate of |
C in pressure coordinates the pressure is just the coordinate of |
100 |
C the tracer point, that is, its constant in time |
C the tracer point, that is, its constant in time |
110 |
end do |
end do |
111 |
end do |
end do |
112 |
end do |
end do |
|
else |
|
|
STOP 'INI_EOS: We should never reach this point!' |
|
113 |
endif |
endif |
114 |
|
|
115 |
if ( equationOfState .eq. 'LINEAR' ) then |
if ( equationOfState .eq. 'LINEAR' ) then |
303 |
end if |
end if |
304 |
|
|
305 |
_BEGIN_MASTER( myThid ) |
_BEGIN_MASTER( myThid ) |
306 |
|
C-- Check EOS initialisation: |
307 |
|
|
308 |
call check_eos( myThid ) |
call check_eos( myThid ) |
309 |
|
|
310 |
|
IF ( buoyancyRelation .EQ. 'OCEANIC' .AND. |
311 |
|
& ( equationOfState .EQ. 'JMD95P'.OR. |
312 |
|
& equationOfState .EQ. 'MDJWF' .OR. |
313 |
|
& equationOfState .EQ. 'UNESCO' ) ) THEN |
314 |
|
WRITE(msgBuf,'(A)') |
315 |
|
& 'S/R INI_EOS: WARNING >>> Wrong Restart !!!' |
316 |
|
CALL PRINT_MESSAGE(msgBuf, errorMessageUnit, |
317 |
|
& SQUEEZE_RIGHT,myThid) |
318 |
|
WRITE(msgBuf,'(2A)') 'S/R INI_EOS: ', |
319 |
|
& 'not correct using Z-coord with EOS function of P' |
320 |
|
CALL PRINT_MESSAGE(msgBuf, errorMessageUnit, |
321 |
|
& SQUEEZE_RIGHT,myThid) |
322 |
|
WRITE(msgBuf,'(A)') |
323 |
|
& 'S/R INI_EOS: WARNING <<< Wrong Restart !!!' |
324 |
|
CALL PRINT_MESSAGE(msgBuf, errorMessageUnit, |
325 |
|
& SQUEEZE_RIGHT,myThid) |
326 |
|
ENDIF |
327 |
|
|
328 |
_END_MASTER( myThid ) |
_END_MASTER( myThid ) |
329 |
|
|
330 |
return |
return |