- Timestamp:
- 2019-07-29T16:47:01+02:00 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/dev_r11265_ASINTER-01_Guillaume_ABL1D/src/OCE/DIA/diawri.F90
r11360 r11363 654 654 CALL histdef( nid_T, "sowindsp", "wind speed at 10m" , "m/s" , & ! wndm 655 655 & jpi, jpj, nh_T, 1 , 1, 1 , -99 , 32, clop, zsto, zout ) 656 !656 ! 657 657 IF( ln_abl ) THEN 658 ! !!! nid_A : 3D 659 CALL histdef( nid_A, "t_abl", "Potential Temperature" , "K" , & ! t_abl 658 CALL histdef( nid_A, "t_abl", "Potential Temperature" , "K" , & ! t_abl 660 659 & jpi, jpj, nh_A, ipka, 1, ipka, nz_A, 32, clop, zsto, zout ) 661 660 CALL histdef( nid_A, "q_abl", "Humidity" , "kg/kg" , & ! q_abl … … 677 676 & jpi, jpj, nh_A, 1 , 1, 1 , -99 , 32, clop, zsto, zout ) 678 677 #endif 679 CALL histend( nid_A, snc4chunks=snc4set ) 680 ! 681 ENDIF 682 ! 678 CALL histend( nid_A, snc4chunks=snc4set ) 679 ENDIF 680 ! 683 681 IF( ln_icebergs ) THEN 684 682 CALL histdef( nid_T, "calving" , "calving mass input" , "kg/s" , & … … 838 836 CALL histwrite( nid_T, "soicecov", it, fr_i , ndim_hT, ndex_hT ) ! ice fraction 839 837 CALL histwrite( nid_T, "sowindsp", it, wndm , ndim_hT, ndex_hT ) ! wind speed 840 !838 ! 841 839 IF( ln_abl ) THEN 842 843 844 845 846 END DO 847 840 ALLOCATE( zw3d_abl(jpi,jpj,jpka) ) 841 IF( ln_mskland ) THEN 842 DO jk=1,jpka 843 zw3d_abl(:,:,jk) = tmask(:,:,1) 844 END DO 845 ELSE 848 846 zw3d_abl(:,:,:) = 1._wp 849 847 ENDIF 850 851 852 853 854 855 856 857 848 CALL histwrite( nid_A, "pblh" , it, pblh(:,:) *zw3d_abl(:,:,1 ), ndim_hA, ndex_hA ) ! pblh 849 CALL histwrite( nid_A, "u_abl" , it, u_abl (:,:,2:jpka,nt_n )*zw3d_abl(:,:,2:jpka), ndim_A , ndex_A ) ! u_abl 850 CALL histwrite( nid_A, "v_abl" , it, v_abl (:,:,2:jpka,nt_n )*zw3d_abl(:,:,2:jpka), ndim_A , ndex_A ) ! v_abl 851 CALL histwrite( nid_A, "t_abl" , it, tq_abl (:,:,2:jpka,nt_n,1)*zw3d_abl(:,:,2:jpka), ndim_A , ndex_A ) ! t_abl 852 CALL histwrite( nid_A, "q_abl" , it, tq_abl (:,:,2:jpka,nt_n,2)*zw3d_abl(:,:,2:jpka), ndim_A , ndex_A ) ! q_abl 853 CALL histwrite( nid_A, "tke_abl", it, tke_abl (:,:,2:jpka,nt_n )*zw3d_abl(:,:,2:jpka), ndim_A , ndex_A ) ! tke_abl 854 CALL histwrite( nid_A, "avm_abl", it, avm_abl (:,:,2:jpka )*zw3d_abl(:,:,2:jpka), ndim_A , ndex_A ) ! avm_abl 855 CALL histwrite( nid_A, "avt_abl", it, avt_abl (:,:,2:jpka )*zw3d_abl(:,:,2:jpka), ndim_A , ndex_A ) ! avt_abl 858 856 #if defined key_si3 859 857 CALL histwrite( nid_A, "oce_frac" , it, ato_i(:,:) , ndim_hA, ndex_hA ) ! ato_i 860 858 #endif 861 862 863 !859 DEALLOCATE(zw3d_abl) 860 ENDIF 861 ! 864 862 IF( ln_icebergs ) THEN 865 863 !
Note: See TracChangeset
for help on using the changeset viewer.