- Timestamp:
- 12/20/17 15:47:20 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
codes/icosagcm/devel/src/unstructured/timestep_unstructured.F90
r642 r645 41 41 dPhi_fast, dPhi_slow, dW_fast, dW_slow) BINDC(caldyn_unstructured) ! OUT : tendencies 42 42 DBL, VALUE :: tau 43 FIELD_MASS :: rhodz, drhodz, pk, berni ! IN, OUT, DIAG44 FIELD_THETA :: theta_rhodz, dtheta_rhodz, theta ! IN, OUT, DIAG45 FIELD_GEOPOT :: wflux, w, geopot, & ! DIAG, INOUT46 dPhi_fast, dPhi_slow, dW_fast, dW_slow ! OUT47 FIELD_U :: u,du_fast,du_slow,hflux, qu! INOUT,OUT,OUT,DIAG48 FIELD_Z :: qv ! DIAG49 FIELD_PS :: ps,dmass_col,mass_col ! OUT,OUT,IN (if eta_mass) or OUT,UNUSED,UNUSED (if eta_lag)43 FIELD_MASS :: rhodz, drhodz, pk, berni ! IN, OUT, DIAG 44 FIELD_THETA :: theta_rhodz, dtheta_rhodz, theta ! IN, OUT, DIAG 45 FIELD_GEOPOT :: wflux, w, geopot, & ! DIAG, INOUT 46 dPhi_fast, dPhi_slow, dW_fast, dW_slow ! OUT 47 FIELD_U :: u,du_fast,du_slow,hflux,Ftheta,qu ! INOUT,OUT,OUT,DIAG 48 FIELD_Z :: qv ! DIAG 49 FIELD_PS :: ps,dmass_col,mass_col ! OUT,OUT,IN (if eta_mass) or OUT,UNUSED,UNUSED (if eta_lag) 50 50 FIELD_UL :: wwuu 51 51 DBL :: time1,time2 … … 71 71 CALL compute_pvort_only(rhodz,u,qv,qu) 72 72 CALL compute_caldyn_slow_hydro(rhodz,theta,u, berni,hflux,du_slow) 73 CALL compute_coriolis(hflux,theta,qu, drhodz,dtheta_rhodz,du_slow)73 CALL compute_coriolis(hflux,theta,qu,Ftheta, drhodz,dtheta_rhodz,du_slow) 74 74 IF(caldyn_eta == eta_mass) THEN 75 75 CALL caldyn_vert(drhodz,rhodz,theta,u, dmass_col,wflux,dtheta_rhodz,du_slow,wwuu) … … 91 91 CALL compute_pvort_only(rhodz,u,qv,qu) 92 92 CALL compute_caldyn_slow_NH(u,rhodz,geopot,W, hflux,du_slow,dPhi_slow,dW_slow) 93 CALL compute_coriolis(hflux,theta,qu, drhodz,dtheta_rhodz,du_slow)93 CALL compute_coriolis(hflux,theta,qu,Ftheta, drhodz,dtheta_rhodz,du_slow) 94 94 IF(caldyn_eta == eta_mass) THEN 95 95 CALL caldyn_vert(drhodz,rhodz,theta,u, dmass_col,wflux,dtheta_rhodz,du_slow,wwuu) … … 114 114 FIELD_MASS :: rhodz, pk, berni ! IN, DIAG 115 115 FIELD_THETA :: theta_rhodz, theta ! IN, DIAG 116 FIELD_U :: u,hflux,qu ! INOUT,DIAG116 FIELD_U :: u,hflux,qu,Ftheta ! INOUT,DIAG*3 117 117 FIELD_GEOPOT :: geopot, w, wflux, w_il ! IN, INOUT, DIAG, DIAG 118 118 FIELD_Z :: qv ! DIAG … … 150 150 CALL compute_pvort_only(rhodz,u,qv,qu) 151 151 CALL compute_caldyn_slow_hydro(rhodz,theta,u, berni,hflux,du_slow(:,:,stage)) 152 CALL compute_coriolis(hflux,theta,qu, drhodz(:,:,stage), dtheta_rhodz(:,:,:,stage),du_slow(:,:,stage))152 CALL compute_coriolis(hflux,theta,qu,Ftheta, drhodz(:,:,stage), dtheta_rhodz(:,:,:,stage),du_slow(:,:,stage)) 153 153 IF(caldyn_eta == eta_mass) THEN 154 154 CALL caldyn_vert(drhodz(:,:,stage),rhodz,theta,u, & … … 169 169 CALL compute_pvort_only(rhodz,u,qv,qu) 170 170 CALL compute_caldyn_slow_NH(u,rhodz,geopot,W, hflux, du_slow(:,:,stage), dPhi_slow(:,:,stage), dW_slow(:,:,stage)) 171 CALL compute_coriolis(hflux,theta,qu, drhodz(:,:,stage),dtheta_rhodz(:,:,:,stage), du_slow(:,:,stage))171 CALL compute_coriolis(hflux,theta,qu,Ftheta, drhodz(:,:,stage),dtheta_rhodz(:,:,:,stage), du_slow(:,:,stage)) 172 172 IF(caldyn_eta == eta_mass) THEN 173 173 CALL caldyn_vert(drhodz(:,:,stage),rhodz,theta,u, &
Note: See TracChangeset
for help on using the changeset viewer.