Changeset 10043
- Timestamp:
- 2018-08-07T12:25:11+02:00 (7 years ago)
- Location:
- branches/UKMO/dev_r5518_cleanup_1d_cpl/NEMOGCM/NEMO/OPA_SRC/SBC
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/dev_r5518_cleanup_1d_cpl/NEMOGCM/NEMO/OPA_SRC/SBC/cpl_oasis3.F90
r10041 r10043 329 329 ! Define 0D (Greenland or Antarctic ice mass) or 1D (river outflow) coupling fields 330 330 IF (srcv(ji)%dimensions <= 1) THEN 331 WRITE(numout,*) "RSRH 0d or 1d define field ",zclname; flush(numout)332 331 IF (nproc == 0) THEN 333 332 … … 352 351 OASIS_In , (/ 0, 0 /) , OASIS_REAL, nerror ) 353 352 END IF 354 WRITE(numout,*) "RSRH 0d field done ",zclname,nerror; flush(numout)355 353 ELSE 356 WRITE(numout,*) "RSRH 2d define field ",zclname; flush(numout)357 354 ! It's a "normal" 2D (or pseudo 3D) coupling field. 358 355 CALL oasis_def_var (srcv(ji)%nid(jc,jm), zclname, id_part , (/ 2, 0 /), & 359 356 OASIS_In , ishape , OASIS_REAL, nerror ) 360 WRITE(numout,*) "RSRH 2d field done ",zclname,nerror; flush(numout)361 357 ENDIF 362 358 … … 376 372 ! End of definition phase 377 373 !------------------------------------------------------------------ 378 WRITE(numout,*) "RSRH NEMO calling enddef";flush(numout)379 374 CALL oasis_enddef(nerror) 380 WRITE(numout,*) "RSRH NEMO finished enddef", nerror;flush(numout)381 375 IF( nerror /= OASIS_Ok ) CALL oasis_abort ( ncomp_id, 'cpl_define', 'Failure in oasis_enddef') 382 376 ! -
branches/UKMO/dev_r5518_cleanup_1d_cpl/NEMOGCM/NEMO/OPA_SRC/SBC/sbccpl.F90
r10041 r10043 485 485 IF( TRIM( sn_rcv_cal%cldes ) == 'coupled' ) srcv(jpr_cal)%laction = .TRUE. 486 486 srcv(jpr_grnm )%clname = 'OGrnmass' 487 srcv(jpr_grnm )%dimensions = 0 ! Scalar field 488 write(numout,*) "RSRH set up grnmss dimension:" 489 487 srcv(jpr_grnm )%dimensions = 0 ! Scalar field 490 488 491 489 IF( TRIM( sn_rcv_grnm%cldes ) == 'coupled' ) srcv(jpr_grnm)%laction = .TRUE. … … 959 957 ! 960 958 CASE( 0 ) ! Scalar field 961 WRITE(numout,*) "RSRH allocate zero dim field z3",jn ; flush(numout)962 963 ! We have a scalar field964 959 ALLOCATE( frcv(jn)%z3(1,1,1) ) 965 960 966 961 CASE( 1 ) ! 1D field 967 968 ! In the special case of 1D fields we can't allocate the array yet as969 ! we don't know what nn_cpl_river is.970 WRITE(numout,*) "RSRH allocate 1 dim field z3",jn, nn_cpl_river ; flush(numout)971 972 962 ALLOCATE( frcv(jn)%z3(nn_cpl_river,1,1) ) 973 963 974 CASE DEFAULT 975 976 WRITE(numout,*) "RSRH allocate 2 dim field z3",jn,srcv(jn)%nct ; flush(numout) 977 ! We have a "normal" 2D (or pseudo 3D) field. 964 CASE DEFAULT ! 2D (or pseudo 3D) field. 978 965 ALLOCATE( frcv(jn)%z3(jpi,jpj,srcv(jn)%nct) ) 966 979 967 END SELECT 980 968 END IF … … 1118 1106 ! 1119 1107 IF( ln_mixcpl ) zmsk(:,:) = 1. - xcplmask(:,:,0) 1120 write(numout,*) "RSRH in cpl_rcv",kt ; flush(numout)1121 1108 ! 1122 1109 ! ! ======================================================= ! … … 1126 1113 DO jn = 1, jprcv ! received fields sent by the atmosphere 1127 1114 IF( srcv(jn)%laction ) THEN 1128 write(numout,*) "RSRH recieving field via call to cpl_rcv",kt,jn ; flush(numout)1129 1115 1130 1116 IF ( srcv(jn)%dimensions == 0 ) THEN 1131 write(numout,*) "RSRH recieving 0d or 1d field",kt,jn ; flush(numout)1132 1117 CALL cpl_rcv_1d( jn, isec, frcv(jn)%z3, SIZE(frcv(jn)%z3), nrcvinfo(jn) ) 1133 1118 ELSE 1134 write(numout,*) "RSRH recieving 2d field",kt,jn ; flush(numout)1135 1119 CALL cpl_rcv( jn, isec, frcv(jn)%z3, xcplmask(:,:,1:nn_cplmodel), nrcvinfo(jn) ) 1136 1120 END IF 1137 write(numout,*) "RSRH completed recieve of field via call to cpl_rcv",kt,jn ; flush(numout)1138 1121 1139 1122 END IF 1140 1123 END DO 1141 write(numout,*) "RSRH still in cpl_rcv",kt ; flush(numout)1142 1124 ! ! ========================= ! 1143 1125 IF( srcv(jpr_otx1)%laction ) THEN ! ocean stress components ! … … 1440 1422 zmask_sum = glob_sum( tmask(:,:,1) ) 1441 1423 ENDIF 1442 write(numout,*) "RSRH still in cpl_rcv at teswt for grnm",kt ; flush(numout)1443 1424 1444 1425 IF( srcv(jpr_grnm)%laction .AND. nn_coupled_iceshelf_fluxes == 1 ) THEN 1445 1426 ! This is a zero dimensional, single value field. 1446 1427 zgreenland_icesheet_mass_in = frcv(jpr_grnm)%z3(1,1,1) 1447 write(numout,*) "RSRH still in cpl_rcv inside teswt for grnm",kt ; flush(numout)1448 1428 1449 1429 greenland_icesheet_timelapsed = greenland_icesheet_timelapsed + rdt
Note: See TracChangeset
for help on using the changeset viewer.