[3] | 1 | ! define time axis |
---|
| 2 | it= kt-nit000+1 |
---|
| 3 | IF( lwp .AND. MOD( kt, nwrite ) == 0 ) THEN |
---|
| 4 | WRITE(numout,*) ' diamld : write NetCDF fields' |
---|
| 5 | ENDIF |
---|
| 6 | |
---|
| 7 | CALL histwrite( nidtrd,"somlttml",it,rmld ,ndimtrd1,ndextrd1) ! Mixed-layer depth |
---|
| 8 | CALL histwrite( nidtrd,"somltemp",it,tml ,ndimtrd1,ndextrd1) ! Mixed-layer temperature |
---|
| 9 | |
---|
| 10 | |
---|
| 11 | ! Temperature trends |
---|
| 12 | ! ------------------ |
---|
| 13 | CALL histwrite( nidtrd,"somlttto",it,ztmltot ,ndimtrd1,ndextrd1) ! total |
---|
| 14 | CALL histwrite( nidtrd,"somlttax",it,tmltrd(:,:, 1),ndimtrd1,ndextrd1) ! i- adv. |
---|
| 15 | CALL histwrite( nidtrd,"somlttay",it,tmltrd(:,:, 2),ndimtrd1,ndextrd1) ! j- adv. |
---|
| 16 | CALL histwrite( nidtrd,"somlttaz",it,tmltrd(:,:, 3),ndimtrd1,ndextrd1) ! vertical adv. |
---|
| 17 | CALL histwrite( nidtrd,"somlttdh",it,tmltrd(:,:, 4),ndimtrd1,ndextrd1) ! hor. lateral diff. |
---|
| 18 | CALL histwrite( nidtrd,"somlttfo",it,tmltrd(:,:, 5),ndimtrd1,ndextrd1) ! forcing |
---|
| 19 | |
---|
| 20 | CALL histwrite( nidtrd,"somlbtdz",it,tmltrd(:,:, 6),ndimtrd1,ndextrd1) ! vert. diffusion |
---|
| 21 | CALL histwrite( nidtrd,"somlbtdt",it,ztmlres ,ndimtrd1,ndextrd1) ! dh/dt entrainment (residual) |
---|
| 22 | IF( l_traldf_iso ) THEN |
---|
| 23 | CALL histwrite( nidtrd,"somlbtdv",it,tmltrd(:,:, 7),ndimtrd1,ndextrd1) ! vert. lateral diff. |
---|
| 24 | ENDIF |
---|
| 25 | #if defined key_traldf_eiv |
---|
| 26 | CALL histwrite( nidtrd,"somlgtax",it,tmltrd(:,:, 8),ndimtrd1,ndextrd1) ! i- adv. (eiv) |
---|
| 27 | CALL histwrite( nidtrd,"somlgtay",it,tmltrd(:,:, 9),ndimtrd1,ndextrd1) ! j- adv. (eiv) |
---|
| 28 | CALL histwrite( nidtrd,"somlgtaz",it,tmltrd(:,:,10),ndimtrd1,ndextrd1) ! vert. adv. (eiv) |
---|
| 29 | z2d(:,:) = tmltrd(:,:,8) + tmltrd(:,:,9) + tmltrd(:,:,10) |
---|
| 30 | CALL histwrite( nidtrd,"somlgtat",it,z2d ,ndimtrd1,ndextrd1) ! total adv. (eiv) |
---|
| 31 | #endif |
---|
| 32 | CALL histwrite( nidtrd,"somlsalt",it,sml ,ndimtrd1,ndextrd1) ! Mixed-layer salinity |
---|
| 33 | |
---|
| 34 | ! Salinity trends |
---|
| 35 | ! --------------- |
---|
| 36 | CALL histwrite( nidtrd,"somltsto",it,zsmltot ,ndimtrd1,ndextrd1) ! total |
---|
| 37 | CALL histwrite( nidtrd,"somltsax",it,smltrd(:,:, 1),ndimtrd1,ndextrd1) ! i- adv. |
---|
| 38 | CALL histwrite( nidtrd,"somltsay",it,smltrd(:,:, 2),ndimtrd1,ndextrd1) ! j- adv. |
---|
| 39 | CALL histwrite( nidtrd,"somltsaz",it,smltrd(:,:, 3),ndimtrd1,ndextrd1) ! vert. adv. |
---|
| 40 | CALL histwrite( nidtrd,"somltsdh",it,smltrd(:,:, 4),ndimtrd1,ndextrd1) ! hor. lateral diff. |
---|
| 41 | CALL histwrite( nidtrd,"somltsfo",it,smltrd(:,:, 5),ndimtrd1,ndextrd1) ! forcing |
---|
| 42 | CALL histwrite( nidtrd,"somlbsdz",it,smltrd(:,:, 6),ndimtrd1,ndextrd1) ! vert. diff. |
---|
| 43 | CALL histwrite( nidtrd,"somlbsdt",it,zsmlres ,ndimtrd1,ndextrd1) ! dh/dt entrainment (residual) |
---|
| 44 | IF( l_traldf_iso ) THEN |
---|
| 45 | CALL histwrite( nidtrd,"somlbsdv",it,smltrd(:,:, 7),ndimtrd1,ndextrd1) ! vert. lateral diff; |
---|
| 46 | ENDIF |
---|
| 47 | #if defined key_traldf_eiv |
---|
| 48 | CALL histwrite( nidtrd,"somlgsax",it,smltrd(:,:, 8),ndimtrd1,ndextrd1) ! i-adv. (eiv) |
---|
| 49 | CALL histwrite( nidtrd,"somlgsay",it,smltrd(:,:, 9),ndimtrd1,ndextrd1) ! j-adv. (eiv) |
---|
| 50 | CALL histwrite( nidtrd,"somlgsaz",it,smltrd(:,:,10),ndimtrd1,ndextrd1) ! vert. adv. (eiv) |
---|
| 51 | z2d(:,:) = smltrd(:,:,8) + smltrd(:,:,9) + smltrd(:,:,10) |
---|
| 52 | CALL histwrite( nidtrd,"somlgsat",it,z2d ,ndimtrd1,ndextrd1) ! total adv. (eiv) |
---|
| 53 | #endif |
---|
| 54 | |
---|