C $Header: /home/ubuntu/mnt/e9_copy/MITgcm_contrib/dgoldberg/streamice/streamice_init_varia.F,v 1.2 2012/07/19 18:46:56 dgoldberg Exp $ C $Name: $ #include "STREAMICE_OPTIONS.h" C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----| CBOP SUBROUTINE STREAMICE_INIT_VARIA( myThid ) C /============================================================\ C | SUBROUTINE STREAMICE_INIT_VARIA | C | o Routine to initialize STREAMICE variables. | C |============================================================| C | Initialize STREAMICE parameters and variables. | C \============================================================/ IMPLICIT NONE C === Global variables === #include "SIZE.h" #include "GRID.h" #include "EEPARAMS.h" #include "PARAMS.h" #include "STREAMICE.h" #include "STREAMICE_CG.h" #include "STREAMICE_ADV.h" C === Routine arguments === C myThid - Number of this instance of STREAMICE_INIT_VARIA INTEGER myThid CEndOfInterface #ifdef ALLOW_STREAMICE C === Local variables === C I,J,bi,bj - Loop counters INTEGER i, j, k, bi, bj, Gi, Gj INTEGER col_y, col_x _RL slope_pos, c1 CHARACTER*(MAX_LEN_MBUF) msgBuf CEOP C ZERO OUT FLOATING POINT ARRAYS DO bj = myByLo(myThid), myByHi(myThid) DO bi = myBxLo(myThid), myBxHi(myThid) DO j=1-Oly,sNy+Oly DO i=1-Olx,sNx+Olx H_streamIce(i,j,bi,bj) = 0. _d 0 U_streamice(i,j,bi,bj) = 0. _d 0 V_streamice(i,j,bi,bj) = 0. _d 0 visc_streamice(i,j,bi,bj) = 0. _d 0 tau_beta_eff_streamice(i,j,bi,bj) = 0. _d 0 float_frac_streamice(i,j,bi,bj) = 0. _d 0 base_el_streamice(i,j,bi,bj) = 0. _d 0 surf_el_streamice(i,j,bi,bj) = 0. _d 0 area_shelf_streamice(i,j,bi,bj) = 0. _d 0 mass_ice_streamice(i,j,bi,bj) = 0. _d 0 C_basal_friction(i,j,bi,bj) = C_basal_fric_const #ifdef ALLOW_AUTODIFF_TAMC ru_old_si(i,j,bi,bj) = 0. _d 0 rv_old_si(i,j,bi,bj) = 0. _d 0 zu_old_si(i,j,bi,bj) = 0. _d 0 zv_old_si(i,j,bi,bj) = 0. _d 0 h_after_uflux_SI(i,j,bi,bj) = 0. _d 0 #endif ENDDO ENDDO ENDDO ENDDO DO j = 1-oly, sNy+oly DO i = 1-olx, sNx+olx DO bj = myByLo(myThid), myByHi(myThid) DO bi = myBxLo(myThid), myBxHi(myThid) cc DO k=1,4 DO col_x=-1,1 DO col_y=-1,1 streamice_cg_A1(i,j,bi,bj,col_x,col_y)=0.0 streamice_cg_A2(i,j,bi,bj,col_x,col_y)=0.0 streamice_cg_A3(i,j,bi,bj,col_x,col_y)=0.0 streamice_cg_A4(i,j,bi,bj,col_x,col_y)=0.0 ENDDO ENDDO cc ENDDO ENDDO ENDDO ENDDO ENDDO C INIT. INTEGER ARRAYS DO bj = myByLo(myThid), myByHi(myThid) DO bi = myBxLo(myThid), myBxHi(myThid) DO j=1-Oly,sNy+Oly DO i=1-Olx,sNx+Olx STREAMICE_hmask(i,j,bi,bj) = -1.0 STREAMICE_umask(i,j,bi,bj) = 0.0 STREAMICE_vmask(i,j,bi,bj) = 0.0 STREAMICE_ufacemask(i,j,bi,bj) = 0.0 STREAMICE_vfacemask(i,j,bi,bj) = 0.0 STREAMICE_float_cond(i,j,bi,bj) = 0.0 ENDDO ENDDO ENDDO ENDDO !ph SELECT CASE (TRIM(STREAMICEthickInit)) IF ( STREAMICEthickInit.EQ.'PARAM' ) THEN WRITE(msgBuf,'(A)') 'initializing analytic thickness' CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, & SQUEEZE_RIGHT , 1) slope_pos = shelf_edge_pos - shelf_flat_width c1 = 0.0 IF (shelf_slope_scale .GT. 0.0) THEN c1 = 1.0 / shelf_slope_scale ENDIF DO bj = myByLo(myThid), myByHi(myThid) DO bi = myBxLo(myThid), myBxHi(myThid) DO j=1,sNy DO i=1,sNx Gi = (myXGlobalLo-1)+(bi-1)*sNx+i Gj = (myYGlobalLo-1)+(bj-1)*sNy+j IF ((Gi.lt.Nx).and.(Gj.lt.Ny)) THEN C IF (flow_dir .EQ. 2.0) THEN IF (.TRUE.) THEN IF (xC(i-1,j,bi,bj).GE.shelf_edge_pos) THEN area_shelf_streamice(i,j,bi,bj) = 0. _d 0 STREAMICE_hmask(i,j,bi,bj) = 0. _d 0 ELSE IF (xC(i,j,bi,bj).GT.slope_pos) THEN H_streamice (i,j,bi,bj) = shelf_min_draft ELSE H_streamice (i,j,bi,bj) = (shelf_min_draft + & (shelf_max_draft - shelf_min_draft) * & min (1.0, (c1*(slope_pos-xC(i,j,bi,bj)))**2)) ENDIF IF (xC(i,j,bi,bj).GT.shelf_edge_pos) THEN area_shelf_streamice(i,j,bi,bj) = rA(i,j,bi,bj) * & (shelf_edge_pos-xG(i,j,bi,bj)) / & (xG(i+1,j,bi,bj)-xG(i,j,bi,bj)) IF (area_shelf_streamice(i,j,bi,bj).gt. 0._d 0) THEN STREAMICE_hmask(i,j,bi,bj) = 2.0 ELSE STREAMICE_hmask(i,j,bi,bj) = 0.0 H_streamice(i,j,bi,bj) = 0.0 ENDIF ELSE area_shelf_streamice(i,j,bi,bj) = rA(i,j,bi,bj) STREAMICE_hmask(i,j,bi,bj) = 1.0 ENDIF c WRITE(msgBuf,'(A,I3,I3,E9.2)') 'thickness', c & i,j,shelf_edge_pos c CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, c & SQUEEZE_RIGHT , 1) ENDIF ENDIF ENDIF ENDDO ENDDO ENDDO ENDDO ! DO bj = myByLo(myThid), myByHi(myThid) ! DO bi = myBxLo(myThid), myBxHi(myThid) ! DO j=1,sNy ! DO i=1,sNx ! STREAMICE_dummy_array(i,j,bi,bj) = ! & REAL(STREAMICE_hmask(i,j,bi,bj)) ! WRITE(msgBuf,'(F3.0)') STREAMICE_dummy_array(i,j,bi,bj) ! CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, ! & SQUEEZE_RIGHT , 1) ! ENDDO ! ENDDO ! ENDDO ! ENDDO ! ! CALL WRITE_FLD_XY_RL("STREAMICE_hmask"," ", ! & STREAMICE_dummy_array,0,myThid) ! CALL WRITE_FLD_XY_RL("h_streamice"," ", ! & H_streamice,0,myThid) ELSE IF ( STREAMICEthickInit.EQ.'FILE' ) THEN IF ( STREAMICEthickFile .NE. ' ' ) THEN _BARRIER C The 0 is the "iteration" argument. The ' ' is an empty suffix CALL READ_FLD_XY_RS( STREAMICEthickFile, ' ', H_streamice, & 0, myThid ) DO bj = myByLo(myThid), myByHi(myThid) DO bi = myBxLo(myThid), myBxHi(myThid) DO j=1,sNy DO i=1,sNx Gi = (myXGlobalLo-1)+(bi-1)*sNx+i Gj = (myYGlobalLo-1)+(bj-1)*sNy+j IF ((Gi.lt.Nx).and.(Gj.lt.Ny)) THEN IF (H_streamice(i,j,bi,bj).GT.0. _d 0) THEN area_shelf_streamice(i,j,bi,bj) = rA(i,j,bi,bj) STREAMICE_hmask(i,j,bi,bj) = 1.0 ELSE area_shelf_streamice(i,j,bi,bj) = 0. _d 0 STREAMICE_hmask(i,j,bi,bj) = 0. _d 0 ENDIF ENDIF ENDDO ENDDO ENDDO ENDDO ELSE WRITE(msgBuf,'(A)') 'INIT THICKNESS - FILENAME MISSING' CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, & SQUEEZE_RIGHT , 1) ENDIF ELSE WRITE(msgBuf,'(A)') 'INIT THICKNESS - NOT IMPLENTED' CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, & SQUEEZE_RIGHT , 1) ENDIF CALL STREAMICE_UPD_FFRAC_UNCOUPLED ( myThid ) _EXCH_XY_RL(H_streamice, myThid ) _EXCH_XY_RL(STREAMICE_hmask, myThid ) _EXCH_XY_RL(area_shelf_streamice, myThid ) CALL WRITE_FLD_XY_RL ( "H_streamice", "init", & H_streamIce, 0, myThid ) CALL WRITE_FLD_XY_RL ( "area_shelf_streamice", "init", & area_shelf_streamice, 0, myThid ) CALL WRITE_FLD_XY_RL ( "STREAMICE_hmask", "init", & STREAMICE_hmask, 0, myThid ) ! CALL STREAMICE_VELMASK_UPD (myThid) ! CALL STREAMICE_UPD_FFRAC_UNCOUPLED ( myThid ) ! CALL STREAMICE_VEL_SOLVE( myThid ) CALL WRITE_FLD_XY_RL ( "U_init", "", & U_streamice, 0, myThid ) CALL WRITE_FLD_XY_RL ( "V_init", "", & V_streamice, 0, myThid ) ! CALL WRITE_FULLARRAY_RL ("H",H_streamice,1,0,0,1,0,myThid) ! CALL WRITE_FULLARRAY_RL ("hmask",STREAMICE_hmask,1,0,0,1,0,myThid) ! CALL WRITE_FULLARRAY_RL ("umask",STREAMICE_umask,1,0,0,1,0,myThid) #endif /* ALLOW_STREAMICE */ RETURN END